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Preface 


The advent of coherent sources at optical frequencies opened new and important 
possibilities for increased capacity and efficiency in communication systems. These 
increases have yet to be realized. The delay is due in part to the time required to 
develop an adequate technological base and in part to the emergence of fundamental 
problems that do not arise at lower frequencies. To assess some of these problems and 
to delineate a path toward their solutions, NASA/ERC and MIT/RLE invited some of the 
leaders in the field of optical communication to a two-week Workshop that was held at 
Williams College. 

Broadly stated, the Workshop was concerned with the knowledge that must be 
acquired to determine the structure and performance of near-optimum optical systems. 
The feasibility of implementing these systems with existing devices and the character- 
istics of necessary, but unavailable, devices were also considered. We hope that the 
Workshop, and this report of its deliberations, will contribute to the further develop- 
ment of optical communication systems. 

The success of the Workshop was ensured by the expertise and energy that all of the 
participants brought to it, but it was achieved more easily because of the guidance pro- 
vided by the Steering Committee. Paul L. Smith's efficiency as business manager and 
general administrator and the efforts of Charles Janke, Mary Owens, Marilyn Pierce, 
and Bunny Frank also contributed greatly to the smooth operation and enjoyment of the 
Workshop. 

We particularly thank Joseph Hull, whose early commitment to our suggestion for 
the Workshop encouraged us to pursue it, and Estil V. Hoversten, whose contribution 
to its preparation, operation, and summation was a principal factor in its success. 

Sherman Karp Robert S. Kennedy 

NASA/Electronics Research Center M. I. T. /Research Laboratory of Electronics 
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I. Introduction 


Optical communication systems for space applications are potentially of great impor- 
tance to the National Aeronautics and Space Administration. The high data rate require- 
ments for space-to-space, space -to -ground, and computer-to-computer communication, 
as well as the requirements for deep-space applications, constitute a strong justification 
for the pursuit of a research program in optical communication. 

Many system studies for deep-space applications have been performed by 
contractors for the National Aeronautics and Space Administration. The coordination 
and effective utilization of this effort has been complicated by the multidisciplinary 
nature of the work. For example, the technical language barrier between the optical 
communication theorist and the quantum physicist, or between the solid-state device 
researcher and the physicist working on the statistical properties of the atmosphere, 
is a significant problem in the advancement of the program. 

In an attempt to improve this coordination, the Massachusetts Institute of Tech- 
nology and NASA jointly organized a Workshop, which was held at Williams College, 
Williamstown, Massachusetts, from August 4 to 17, 1968. A selected group of partici- 
pants from the areas of the physical sciences having a direct bearing on optical com- 
munication were invited to participate under the guidance of communication engineers. 
Group interactions were emphasized rather than individual work. The partici- 
pants were selected from industry, universities, and NASA centers. The object 
of the Workshop was to identify some of the problems that represent the most 
significant obstacles in the development of optical communications systems to 
their full potential. 

To provide a common basis for work, the first two and one -half days were devoted 
to the presentation of status papers on various aspects of optical communication by mem- 
bers of the group. Abstracts of these papers are presented in Appendix A. After these 
presentations, the participants formed working groups pertaining to three classes of 
channels: atmospheric, cloud, and free-space. In each of these groups, the objective 
was to assess the present level of understanding and to identify problem areas in which 
research should be concentrated to enhance this understanding. This report is intended 
to present the conclusions of the working groups. 

The emphasis of the Workshop was directed toward digital, rather than analog, sig- 
naling schemes, partly because digital systems provide greater flexibility in signal 
selection and partly because they are more easily understood. Many of the conclusions 
are applicable to both digital and analog systems. The approach to the problems reflects 
the perspective of contemporary communication engineering. The following questions 
are typical of those that were posed. 

1 . What is and must be known about the channel in order to determine the optimum 
receiver and transmitter? How can the missing information be acquired? 

E. Given the channel description, is the (mathematically) optimum receiver already 
known? If not, how can it be determined? 
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3. What new devices are needed if near-optimum systems are to be realized? What 
performance can be realized with existing devices? What signal sets should be used? 

4. Given the (mathematically) optimum receiver, can it be implemented easily? If 
not, what are some possible avenues along which "good" practical systems might be 
developed? 

The principal conclusions and recommendations of the Workshop are summarized 
in Section II. In Section III system problems associated with free-space channels are 
discussed. Many of the problems encountered in the free-space system are common 
also to the atmospheric and cloud channels. Section IV represents the largest group 
effort of the Workshop; the atmospheric channel considered there poses many challenging 
problems, as indicated by the suggested work areas. In Section V the problem of com- 
municating through clouds is considered. 



II. Summary and Conclusions 


Several general conclusions, which are applicable to all of the channels discussed, 
emerged from the workshop. We present them first and then consider the atmospheric 
and cloud channels. The free-space channel is not discussed separately, since the con- 
clusions reached by that group applied to all three classes of channels. 

For emphasis, sentences suggesting important areas for further investigation have 
been underscored. These suggestions reflect the bias of the participants, as well as the 
objectives of the workshop. Omission of an area in no way implies that it is unimportant. 

2. 1 GENERAL CONCLUSIONS 

2. 1. 1 Large Collecting Optics 

With proper design, the performance of a system operating through free space, the 
atmosphere, or a cloud will continue to improve as the collecting area of the receiver 
is increased (for a given field of view). The investigation and development of large- 
aperture systems should be intensified. 

The improvement in performance that can be gained by proper design (diversity com- 
bining) is substantial and does not require that the collecting optics be diffraction-limited. 
As an example, for a given aperture area, A^, and a given performance, a heterodyne 
system employing proper combining will require on the order of A c /A r times the trans- 
mitter power required for a (single) heterodyning system, where A q is the area over 
which the field is coherent in the receiver aperture. Comparable performance can be 
realized with direct-detection systems. Difficulties in implementation occur because 
either a large number of parallel heterodyne receivers or a very narrow -band optical 
filter and large detector area is required. These difficulties will be discussed further 
in connection with device needs, 

2. 1. 2 Device Needs 

A continued development of space qualified optical components such as sources, mod- 
ulators, and detectors is vital to the progress of optical communication systems. It is 
important, however, to recognize the limitations of the components which will emerge 
from this development and the importance of some devices not now under development. 

Signal Generation: A wideband modulator is certainly a necessity for wideband FM 
and is very useful in many other communication systems, but it is not a necessity for 
high-rate digital systems. The required set of signals can often be generated by 
exploiting the pulse properties of lasers (mode locking, Q-switching, etc.). In fact, such 
direct generation has some distinct advantages, when it can be employed, because the 
electronic bandwidth of the drivers need only be the information bandwidth rather than 
the channel bandwidth. The point here is that the extreme versatility of a wideband 
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optical modulator is not always required. When it is not required substantial system 
simplifications often can be realized by using a signal-generation technique that is tai- 
lored to the application; for example, using a fast optical switch to select one of M pos- 
sible waveforms. The development of specialized signal generation techniques should 
be pursued. 

Signal Processing: Some promising systems for communication through the atmo- 

sphere would employ devices that have yet to be developed (see section 4. 3). These same 
systems would be useful through clouds, or through free space with large receiving 
antennas that are not diffraction- limited (see sections 3.3 and 5.2). For example, one 
near-optimum system uses frequency and spatial filtering to reject as much background 
noise as possible, then independently detects the signal energy arriving from different 
angles, and finally combines the detector outputs. Such systems can be implemented 

as a collection of parallel direct detection or heterodyning systems, but the number of 

3 8 

subsystems required is often enormous (10 -10 ). An alternative implementation would 
consist of large focusing optics, an aperture stop, a single very narrow-band optical 
filter, and an ideal photon counter. 

At present, neither of these implementations is very appealing. If the paralleling 
approach is to be practical, techniques for the fabrication of arrays containing a large 
number of parallel systems must be developed. In principle, this might be done by a 
microelectronic array in which each element consists of an optical detector followed by 
an IF filter and a square-law detector. The processed outputs would then be added 
together. No such cascades have been developed for parallel operation. 

The second implementation avoids the array problem by employing a large area 
photodetector preceded by an optical filter. For the atmospheric and cloud channels, 
the product of the filter's area and field of view should be large relative to that of a 
diffraction-limited system, while its bandwidth should be very "small." Such devices 
are possible, in principle, and have been realized for some isolated applications by 
exploiting atomic effects. The optical filters now available fall far short of the ideal 
requirements of communication systems. The possibility of exploiting novel physical 
phenomena to gain improved filter performance should be explored. 

Another obstacle to the development of systems along this second avenue is that rela- 
tively large-area, semi-ideal, photon counters are required. At present, these are only 
available at wavelengths shorter than ~0. 9 pm. Thus an infrared photon counter would 
be very desirable. 

Receivers and/or transmitters that measure and compensate for the phase distortion 
introduced by atmospheric turbulence, clouds, and antenna imperfections are also prom- 
ising (see sections 3. 3 and 4. 3. 3). The implementation of such systems requires, in 
essence, a lens, or phase grating, whose phase delay can be varied in time and space. 

As in the system just discussed, this variation can be achieved by paralleling a large 
number of simple systems, for example, the active mirror concept. It may be possible, 
however, to obtain this variation more simply and efficiently at the focal plane of a lens 
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with a distributed electro- optic system. The feasibility of a temporally and spatially 
adjustable phase grating should be explored. 

Another promising system employs a transmitted reference, or "pilot tone," which 
is filtered out of the signal at the receiver and "spatially amplified" with the spatial 
phase kept intact. In some situations, the amplified reference can be used as a local 
oscillator which is in phase with the signal at every point on the detector surface (see 
sections 3. 3 and 4. 3. 3). Consequently, the need for a large number of parallel systems 
or for a narrow-band optical filter for, say, the atmospheric channel can be avoided. 

The characteristics of wide-angle optical spatial amplifiers should be explored. 

Multiplexing: The electronic bandwidth of existing and contemplated modulators and 

detectors significantly limits the information rate that can be realized with these devices. 
When higher information rates are required, multiplexing of the optical waveforms will 
be necessary. The development of optical multiplexing techniques and devices is prob- 
ably a prerequisite to operation at data rates exceeding 10 bits/sec. 


2.1.3 Systems Analysis 

Many valuable studies of optical communication systems, particularly free-space 
systems, have been undertaken. The quantitative understanding of such systems is still 
meager. Thus, just as a continued improvement of optical devices is important to opti- 
cal communication, so also is a mor e complete understanding of system performance 
characteristics . 

The inadequacy of our present understanding is due, in part, to the novel charac- 
teristics of some optical channels, for example, the atmosphere; but the dominance of 
quantum effects at optical frequencies has been a more fundamental obstacle to quanti- 
tative results. Thus, although their output signal-to-noise ratios are comparable, 
direct detection and heterodyning systems yield different performance levels. We shall 
very briefly summarize the present understanding of these two systems and then consider 
the problem of optimum quantum systems. 

Heterodyning Systems: Broadly stated, the analysis of heterodyning systems poses 

no fundamentally new difficulties. This is not to say that all systems of interest have 
been, or can easily be, analyzed. The required analysis techniques are well under- 
stood, and results are available (see sections 3. 2. 1 and 4. 3. 2) for some important sit- 
uations. These analyses have made it possible to compare different signaling and post 
heterodyning processing schemes. 

Direct Detection Systems: Unlike heterodyne detection, the analysis of direct detec- 

tion systems involves a class of problems that have not been investigated previously. 
Consequently, the performance of only a few direct-detection systems is now known (see 
sections 3. 2. 1 and 4. 3. 1). The lack of quantitative results has made it difficult to deter- 
mine optimum waveforms and post detection processing. It has also impeded perfor- 
mance comparisons between heterodyning and direct-detection systems. 
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Optimum Quantum Systems: It appears, at least in the binary case, that the optimum 

receiver for combating quantum noise is sometimes neither a direct-detection nor a het- 
erodyning system. On the other hand, in all of the available results, either the direct- 
detection or heterodyning system is essentially as good as the optimum system. Thus, 
the question arises: Is the optimum quantum system ever significantly superior to both 

the heterodyning and direct-detection system ? 

If it can be shown that either a direct -detection or heterodyning system is "essen- 
tially" optimum for all channels of interest, the question becomes relatively unimportant 
to practical space communication. If, however, any example can be found for which the 
optimum quantum receiver is far superior to heterodyning and direct detection, the ques- 
tion becomes extremely significant. Even if essentially optimum systems can be real- 
ized by heterodyning or direct detection, the understanding of such systems will be 
enhanced and simplified when a unified quantum communication theory is available. 

In addition to the general conclusions just discussed, there were conclusions per- 
taining specifically to the atmospheric and cloud channels. 

2. 2 ATMOSPHERIC CHANNEL 

While a great deal of effort has been expended, for several years past, on the prob- 
lem of obtaining an adequate theoretical and experimental understanding of the effects 
of atmospheric turbulence on optical propagation, the present state of knowledge is still 
not entirely satisfactory from a communication theory viewpoint. In particular, at the 
present time there is no general consensus of opinion on an adequate theoretical solu- 
tion to the turbulent optical propagation problem. The workshop participants felt that 
the Rytov approximation (see section 4. 1. 1) was often adequate for communication sys- 
tem studies, particularly when an Earth-Space path was involved; this feeling was based 
largely on the reasonable agreement between this theory and existing measurements. On 
the other hand, there was no great confidence that this approximation technique ade- 
quately extrapolates to answer new questions where confirming measurements are not 
available. A better u nderstanding of the limitations, as far as answering questions of 
interest in communication system design, of the Rytov technique is required; also, solu- 
tion techniques that are bot h complete, in the sense of answering the communication 
questions of interest, and also more general, in the sense of being applicable to most 
physical situations, should be explored. Specific theoretical questions of interest in the 
channel characterization problem are discussed in section 4. 1. 5. 

At the present time, data on the microscale changes of the atmosphere's refractive 
index are scarce and their interrelationships with other atmospheric conditions are not 
well understood. There are almost no direct data on the altitude variation of the refrac- 
tive index process. The lack of confidence in the theoretical propagation results has 
also made indirect refractive index variation measurements suspect. Effort should be 
expended to understand the relationship of the refractive index variation process to other 
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atmospheric conditions and to better characterize this process, for example, examine 

the time behavior and ergodicity of the process. Measurements of the height distribu- 
2 

tion of C N , the refractive index structure constant, should be made; the measurements 
should include direct temperature measurements, and if feasible, the development of 
some remote sensing technique. 

Measurements of various parameters of the received optical field have been made. 
Many of these measurement programs have been concentrated on a single or small num- 
ber of parameters, however, and often have not been guided by a communication objec- 
tive. In particular, there are relatively few phase measurements available. There is 
a need f or a comprehensive measurement program with the objective of supplying the 
information required to design and analyze communication systems . The emphasis 
should be on characterizing as completely as possible the received field, including the 
effects of variations in various transmitter and receiver parameters. Thus such quan- 
tities as the statistics of the received intensity and phase processes and the intensity and 
phase coherence bandwidths should be measured. Furthermore, the variations in the 
statistical properties of the intensity and phase processes with such parameters as the 
transmitting aperture size and illumination, the receiver aperture size -and field of view, 
and the source wavelength are needed. Particular measurements of interest are listed 
in section 4. 1. 5. 

The measurements discussed above should be made over propagation paths that are 
typical of the communication links that are of interest. In particular, data are needed 
for ground-space paths if a ground-space link is the one of interest. With a better theo- 
retical understanding of turbulent propagation and measurements of the turbulence alti- 
tude profile, it may be possible to extrapolate careful ground measurements to ground- 
space paths. It will always be preferable to have measurements over the actual path, 
however. Thus careful consideration should be given to using stellar measurements 
when they are applicable. Such consideration should include performing additional stel- 
lar measurements if it is found that appropriate data are not now available, for example, 
statistical measurements and daylight measurements. Obviously as sources and detec- 
tors become available in space they should be used to make the desired measurements. 
Indeed the sources and detectors put into space in the near future should be so designed 
that they maximize the number of meaningful channel characterization measurements 
that can be made . 

The analysis and synthesis of optical communication systems for the atmospheric 
channel has been impeded by the lack of a tractable representation for the received field 
process. Often the appropriate process model for the received field is, conditionally, 
the sum of log-normal and Gaussian processes, which presents analysis difficulties, 
Work is needed on altern ative channel representations and appropriate approximations. 

In particular, the relationship between aperture and focal plane representations should 
be examined with emphasis on the form in which the degrees of freedom appear. An 
angular spectrum representation for the received field should be investigated both 
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theoretically and experimentally. (See section 4. 2 for more details.) It should be empha- 
sized that this channel representation work is not independent of the measurement pro- 
gram suggested above, but rather should both be guided by measurement results and 
used to suggest appropriate measurement and analysis questions. 

Communication systems for the atmospheric channel must process the received field 
in such a way as to properly treat the spatial diversity and combat the fading. For both 
direct-detection and most heterodyning systems, the best processing will require an 
array of detectors. At the present time, however, neither the general performance of 
optimized single -detector systems nor the performance of array systems is known for 
the case of direct detection; further, the optimum processing is also, in general, 
unknown. Thus the tradeoff between performance and complexity for direct-detection . 
systems is not known. 

For heterodyning systems, more results are available. Some performance results 
exist for both single -detector and array systems. These results indicate that proper 
spatial diversity combining is critical when heterodyne detection is used and that useful 
aperture sizes may be relatively small or the number of array elements large, particu- 
larly at visible frequencies. Phase-front tracking receivers can be used to increase the 
effective aperture area for a single heterodyne detector (by approximately a factor of 10) 
or reduce the number of array elements required for a given aperture. Heterodyning 
systems also have many unanswered questions. For example, the data processing is 
not known for signaling times that are long relative to the channel coherence time, or 
for analog modulation. As discussed in section 2. 1. 2, the important question of prac- 
tical implementation of heterodyning array systems is unsolved. Analysis work should 
be performed to provide a better understanding of both direct-detection and heterodyning 
systems for the atmospheric channel. As the design tradeoffs become better understood, 
experimental work should be performed to validate the approximations and assumptions 
that went into the system designs. 

The use of adaptive systems for the atmospheric channel was discussed at the Work- 
shop. Such systems include both channel-measurement receivers and adaptive trans- 
mitters (see section 4. 3). Channel-measurement techniques look promising, because 
of the channel's long coherence time. With respect to measurement receivers, the 
important issues are how well the channel variations can be measured, the performance 
improvement that is possible by using these noisy measurements, the system imple- 
mentations required, and the signal design problem. 

With respect to transmitter adaptation, the question is; To what degree can a trans- 
mitter compensate for distortions that are introduced by the channel? The answer 
involves the issues of channel measurement and reciprocity, as well as transmitter 
strategy. Specific examples of adaptive transmitter techniques are variable-rate sys- 
tems, adaptive signal design systems, and adaptive antenna (beam -shaping) techniques. 
The performance improvement that might be realized by both channel-measurement 
receivers and adaptive transmitter techniques and practical means of implementing them 
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should be investigated. The questions of channel "state" determination should be studied 
and measurements should be made to determine the validity of the concept of atmospheric 
reciprocity (see section 4. 3). 

The use of transmitted reference signals and 'imaging' preamplifiers, as discussed 
in sections 3. 3 and 4. 3, are techniques with promise for the atmospheric channel. Thus 
transmitted reference signaling schemes and imaging preamplifier receiving techniques 
should be studied further to determine their performance and feasibility. 


2.3 SCATTERING CHANNELS 

The feasibility of optical communication through strongly scattering media such as 
clouds, fog, and haze appears to have received very little attention. The available 
knowledge and experience with related channels suggests, however, that such communi- 
cation may be feasible — at least at modest information rates. 

The effect of a given nonabsorbing scattering medium depends upon its location rela- 
tive to the system terminals. If it is located "near" the transmitter, its major effect is 
to spread the transmitted beam and hence to reduce the energy that can be collected by 
the receiver. This reduction imposes a basic limitation upon the potential for communi- 
cation through scattering media near the transmitter. If the scattering medium is located 
"near" the receiver, its major effect is to destroy the spatial and temporal coherence 
properties of the received field; however, its presence does not lead to an enormous 
reduction in received energy if the receiver is properly designed. Consequently, a non- 
absorbing cloud near the receiver, for example, space-to-ground communication, does 
not limit the potential for communication as severely as does a cloud near the trans- 
mitter, for example, ground-to-space communication. Since all-weather operation of 
optical systems is sometimes important, the feasibility of communication through clouds, 
fog, and haze, should be investigated. 

The initial effort should be to obtain a description of the radiation emerging from 
these channels that is adequate for the analysis of communication systems. Some means 
of obtaining this description are discussed in Section V. The difficulties of an exact and 
complete analytical description dictate that approximate analysis, simulation, and experi- 
mental data be employed. Even in the absence of a complete description, valuable esti- 
mates of the potential for communication through clouds, fog, and haze, can and should 
be obtained from some simple experiments (section 5. 3). 

A question that emerged at the Workshop concerned the absorption coefficient of water 
and ice at 10. 6 pm. The tabulated values of absorption are larger than the values sug- 
gested by recent measurements on natural fog and rain, as well as on artificial clouds. 
This uncertainty is important because the tabulated data suggest that communication 
through clouds at 10. 6 pm is not competitive with visible wavelengths, while the recent 
measurements suggest that it may be. Experiments to resolve the uncertainty in the 


absorption of ice and water at 10. 6 pm should be performed. 
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A thorough investigation of the design and performance of systems for scattering 
channels must await a more complete description of the channel characteristics. Pre- 
liminary estimates of the performance characteristics and device requirements for opti- 
cal communication through clouds, fog, and haze can, and should, be performed. The 
general structure of these systems will probably be similar to atmospheric systems, 
but the device requirements and performance limitations will be more severe (sec- 
tion 5.1). 


III. Free-Space Systems 


Free- space systems, as typified by satellite-to- satellite communication, were not 
major topics of discussion at the Workshop, but many of the unresolved problems 
encountered in atmospheric and cloud systems are also encountered, in a simpler form, 
in free-space systems. Because of this, and since these systems are important in their 
own right, we shall first consider free-space systems. 

By a free-space system we mean one in which the propagating medium only attenu- 
ates the transmitted signal, for example, interplanetary space or a nonturbulent atmo- 
sphere. Such systems are in many ways the optical equivalent of line-of-sight 
microwave systems. The quantum and pointing issues, however, that are unimportant 
at microwave frequencies become fundamental at optical frequencies. 

At microwave frequencies the limitations upon antenna size and beamwidth are such 
that adequate pointing accuracy can be attained without limiting the communication sys- 
tem's performance; at optical frequencies pointing accuracy and antenna stability impose 
important constraints upon the usable antenna sizes. Although this suggests that antenna 
pointing should be considered as an integral part of the communication system, a com- 
prehensive discussion of the pointing problem was beyond the scope of this Workshop 
and the problem was not discussed in detail. 

Traditionally, it has been possible to view a "radio receiver" as a device that adds 
some noise to its input, or inputs, and then performs noise-free operations (mixing, fil- 
tering, detection) upon the resultant sum. The important point is that the nature 
and intensity of the added noise is independent of the operations performed. In particu- 
lar, it has often been possible to characterize a receiver by a temperature that is a 
measure of the quality of the receiver front end. Consequently, questions of optimum 
system design could, to a great extent, be separated from questions of receiver noise 
characteristics. At optical frequencies the situation is fundamentally different because 
of quantum effects. In particular, the character of the noise associated with a receiver 
is strongly dependent upon the operations that it perferms. Moreover, this situation 
is a consequence of fundamental (quantum) physical considerations; it cannot be relieved 
by improved component design. 

Three approaches to optical communication systems have been fruitful thus far. The 
first approach employs heterodyne detection and introduces a quantum noise that is a 
Gaussian random process. Broadly stated, these systems are well understood. The 
second approach uses energy-detection devices and introduces a "quantum noise" that 
is a shot- noise process. There is renewed interest in problems of signal design, detec- 
tion, and parameter estimation for such processes. The final approach, and probably 
the one least understood, is the reformulation of the communication problem in quantum- 
mechanical terms. As we shall discuss, there are fundamental problems in this 
approach, but it appears to be the only way of determining the ultimate performance 
limits of optical communication systems. Before considering this approach we shall 
review the present technology in components and receiver structures. 
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3.1 EXISTING RECEIVER COMPONENTS 

All optical receivers that have been constructed, or contemplated, have employed 
direct detection (energy detection), heterodyning, homodyning, or preamplification. Since 
heterodyning and homodyning have only been implemented efficiently with energy detec- 
tion and amplification, we only consider energy detectors and amplifiers. Components, 
such as dielectric filters, which are well understood, are not discussed nor are some 
of the sophisticated components, such as parametric converters and second- harmonic 
generators . 


3.1.1 Energy Detectors 

Many parameters of an energy detector are useful in the design of optical communi- 
cation systems: sensitivity, bandwidth, internal noise power, area, and so forth. A 
more complete statistical description of the detector is required, however, if its ultimate 
utility in a communication system is to be judged. One of the questions discussed at this 
Workshop was whether or not such a description exists, in principle, for the various 
classes of available detectors. The conclusion was that it does — at least in broad out- 
line. Since this general description has not been presented in a form that is accessible 
to many communication engineers, we present it here. 

Ideal Energy Detectors: It is generally agreed that the output of an ideal energy 

1-5 

detector of quantum efficiency p is a conditional Poisson process. That is, given a 
linearly polarized classical field , E(r, t) = Re [V(r, t) exp j2irvt], incident upon the detec- 
tor, the sequence of detections forms a time-varying Poisson process of events whose 
rate parameter \(t) is 


Mt) = 


P / |v(r, t )| 2 dr 
_ — 


( 1 ) 


where the integral is over the detector area. If the process V(r,t) is not narrow-band, 
one should also include the wavelength dependence in p and v. An equivalent descrip- 
tion of the detected process could be stated thus: Given the incident field, the numbers 

of detected photons in any two disjoint time intervals are independent of each other, and 
the probability of detecting k photons in the time t to t+T is 


K k — 

P(k) = |t exp -K, 


(2) 


where 


If a photon description of the field is used, an identical result can be obtained by con- 
sidering the detector as a binary branching process. 
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_ ot+T 

K=\ X(t) dt. (3) 

J t 

It is important to note that the classical field incident upon the detector is presumed 
to be known. If it is random, the conditional Poisson process must be averaged over the 
statistics of the incident field. Thus, for example, if the incident field is an extremely 
narrow-band, zero-mean, Gaussian random process, the probability of k counts in a 
time that is small compared with the inverse of the optical bandwidth will obey a Bose- 
Einstein distribution. 

Real Photon Detectors : Real detectors differ from ideal detectors in several ways. 
Although the details of these differences depend upon the specific detector that is con- 
sidered, it appears that the following characterization is generally accepted by 
physicists. 

Given the classical field incident upon the detector, the output of the detector and its 
associated amplifier can be expressed as 

f(t) = h(t-t i> G) + n T (t) + n g (t) (4) 

i 

and can be interpreted as illustrated in Fig. 1. In this expression n T (t) is a thermal 
noise (zero-mean, stationary, and Gaussian) that is typically white over the band of 
interest. The t. are the detection times for an ideal photon detector, that is, given \(t) 
in Fig. 1, they comprise a Poisson random process with rate parameter \(t). The func- 
tion h(t-G, G) is the shape of the detector output pulse which is caused by a photon 


POISSON 
PROCESS 
RATE — 
PARAMETER 
X(t) 


LINEAR, RANDOM, 
TIME-VARYING FILTER 
IMPULSE RESPONSE h(t-r,r) 


ijM r» s (t) 




f(t) 


r 


POISSON PROCESS 

Fig. 1. Photon detector model. 


detected at time G . In general, this pulse can itself be a random function. As is shown 
in Fig. 1, 2 h(t-G, t.) can be interpreted as the output of a random time -varying filter 

driven by a conditionally Poisson process. Finally, n (t) is a signal-independent "shot 

s 

noise" that can also be viewed as the output of a random time -varying filter driven by 
an (unconditional) Poisson process. This can represent detector dark current. 



14 


OPTICAL SPACE COMMUNICATION 


To illustrate the application of Eq. (4), we consider several examples in the absence 
of thermal noise and dark current. 

i. Ideal photon detector: h(t-h, t.) = 6(t-t.), where 6( • ) is a unit impulse. 

ii. Photodiodes: to first order h(t-C, t^) = g(t-T), where g( • ) is a fixed known func- 
tion whose Fourier transform is the detector frequency response. 

iii. Photomultiplier : to first order h(t-T, t^) = cr^gft-F), where the set of a\ are inde- 
pendent, identically distributed, random variables. 

iv. Photoconductor: to first order h(t-t., t.) = g[(t— t^l/cr^], where the cr^ are indepen- 
dent, identically distributed, random variables. The probability density of cn is 

1 tr - 
P(®’ i ) = ~ exp - — . 

v. Avalanche Photodiode: Similar to iii, except that the probability density of o\ differs. 
Typically, the variance of cr^ is of the form (<r) a with a < 2, where o r is the mean value 
of <r. . 

l 

For a first approximation, n^,(t) is the dominant additive noise in detectors with 

small internal gain (photoconductors), while n g (t) is the dominant internal noise in 

devices with large internal gain (photomultipliers and avalanche diodes). At present, 

photodiodes and photomultipliers are limited to wavelengths less than approximately 

0.9 pm, avalanche detectors are limited to wavelengths shorter than approximately 

5 pm, and photoconductors are limited to wavelengths shorter than approximately 

200 pm. It is reasonable to say that photomultipliers are close approximations to 

ideal photon detectors with a quantum efficiency of less than approximately 30% . The 

6 12 

detailed noise characteristics of available devices have been discussed extensively. 

3.1.2 Preamplifiers 

Although not used extensively, optical preamplifiers of modest gain are available at 
all laser wavelengths (every laser can be used as an amplifier). For an ideal, high-gain, 
four-level preamplifier, the equivalent input noise can be considered to be zero-mean 
and Gaussian with a power density of hv W/Hz per spatial mode . [if the field of view 
is too large, a high gain cannot be realized.] The total noise power per Hz is then 
hv A«/2 where A is the area, X. is the wavelength,- and n is the solid-angle field of 
view of the receiver. The noise levels per mode for lower gain, or three-level, devices 
are larger than hv. 

3.2 EXISTING RECEIVER STRUCTURES 

Existing receivers employ optical amplifiers, local oscillators, and energy detectors 
in one of the following system configurations: heterodyning, homodyning, or direct 

All power densities are unilateral or "measured." 
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detecting. Since the general characteristics of these systems have been discussed exten- 

10 14-24 

sively, we review them only briefly. ’ This review will provide a basis for the 

subsequent discussion of some important unresolved questions pertaining to all of the 
channels that we consider. 

3.2.1 Heterodyning 

At present, heterodyning is accomplished by adding a strong local- oscillator field 
to the received field, detecting the sum with an energy detector, and passing the detected 
waveform through a bandpass filter. If the received and local-oscillator fields are 
linearly polarized plane waves propagating in the same direction, and if the local 
oscillator is strong enough, the output signal will, in effect, be a frequency-translated 
replica of the incident "waveform" plus a zero-mean, stationary, Gaussian noise of 
power density hv, where v is the carrier frequency. The assumption here is that the 
bandwidth of the received field is small relative to the carrier frequency. 

If the local-oscillator field is a plane wave, but the received field is not, the sys- 
tem, in essence, extracts those plane-wave components of the received field whose 
angle of propagation differs from that of the local oscillator by no more than an angle 
\/D, where K is the wavelength, and D is the diameter of the heterodyning area (which 
is assumed to be circular) . More generally, one can say that heterodyning extracts the 

temporal variation of one spatial mode of the received field, translates it to an inter- 

25 26 

mediate frequency, and adds to it a zero-mean Gaussian noise of power density hv. ’ 

The particular mode that is extracted depends upon the spatial variation of the local- 
oscillator field. 

One of the important characteristics of a heterodyning system is that internal detec- 
tor noise (n (t) and n„(t) in Eq. (1)) is unimportant. It is only the quantum, or zero-point 
fluctuation, noise that contaminates the received field. Thus, for example, if the input 
"signal" power in the relevant mode is P, the output signal- to- noise-power-density 
ratio will be P/hv, provided that the local oscillator is strong enough. 

3.2.2 Homodyning 

A homodyning system can be viewed as a heterodyning system in which the intermedi- 
ate frequency is zero. That is, the original optical field variation is translated to base- 
band. If the system is coherent, that is, the local- oscillator field is in phase with the 
signal component of the received field, the equivalent input noise power density is 
hv/2 W per Hz. Roughly stated, the noise power reduction can be attributed to the phase 
information that is available to the system. 

3.2.3 Direct Detection 

By a direct detection system we mean one composed of collecting optics, a narrow- 
band filter, and a photon detector. Conditioned upon the received field the system out- 
put is described as in section 3.1.1. This description is of limited utility because the 
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received field itself usually contains a noise process, and one desires the output sta- 
tistics unconditioned upon the noise process. In general, it is difficult to obtain these 
statistics, but it is possible to define a signal-to-noise ratio that is sometimes a useful 
measure of performance. 

Specifically, suppose that the received field consists of a signal of instantaneous 
envelope power P g plus a zero-mean background noise of power density P^ W/Hz. Also, 
suppose that the total internal noise of the detector is of mean zero and power density 
P n W/Hz. Subject to some reasonable assumptions, the signal-to-noise ratio per Hz 
referred to the system output will be 


,g 2 p 2 

1 S 


2hvB[G 2 (P s + P b )+P n ] 


(5)- 


where B is the signal bandwidth; q is the quantum efficiency, and G and G are, respec- 

27 28 

tively, the mean and second moment of the detector gain. ’ 

It is clear from (5) that the detector internal noises will influence the signal-to-noise 
ratio unless 

P «G 2 (P+PJ. 
n s b 


At present, it is possible to satisfy this condition for wavelengths shorter than approxi- 
mately 1 pm. For longer wavelengths, the internal noises are significant, even at detec- 
tor temperatures of 4°K. That is, near ideal photon counters are available at wave- 
lengths shorter than 1 pm, but not at longer wavelengths. [Note that we are ignoring 
the important questions of detector bandwidth and quantum efficiency.] It is also worth 
noting, for future reference, that the detector areas available at wavelengths shorter 
than 1 pm are much larger than those available at longer wavelengths. 


3.2.4 System Performance 

We shall now briefly summarize the state of knowledge about the performance of opti- 
cal communication systems operating through a free-space channel. The emphasis will 
be upon important questions that remain to be answered: the corresponding questions for 
other channels will be discussed eventually. 

The performance analysis of heterodyning and homodyning systems is, in essence, 
the analysis of systems employing known signals in additive white Gaussian noise. Such 
systems have been studied extensively during the past two decades and are now well 

understood. Thus the optimum signal processing, the optimum signal sets, and the per- 

29-34 

formance of optimum and many suboptimum systems have all been determined. 

The performance of direct-detection systems is a subject of recent investiga- 
35-37 

tions. Most of these investigations have been addressed to problems in which the 

detector output is a Poisson process when conditioned upon the signal component of 
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the received field but averaged over the background noise. This situation will prevail 
for an ideal photon detector, or a detector whose internal noise is Poisson, if there is 
no background noise . It will also prevail if the predetection noise bandwidth is much 
greater than both the post detection bandwidth and the signal bandwidth, as it often is 
for daytime operation with existing optical filters. These investigations have established 
the optimum post detector processing, optimum signal sets, and system error proba- 
bility for the optimum and some suboptimum systems . 

Some investigations have also been addressed to direct detection of specific signal 
sets when the noise bandwidth is comparable with the signal bandwidth, that is, when 
proper predetection filtering is employed. 38 ” 40 As would be expected, the system per- 
formance improves as the filter bandwidth approaches the signal bandwidth. Although 
appropriate filters do not now exist, one can compensate for the filter inadequacies by 
proper signal design, at least in part. 

We shall not attempt a detailed comparison of heterodyne, homodyne, and direct- 
detection systems. A few general comments may be of value. 

Roughly stated, homodyne detection is superior to direct detection for small- alphabet, 
40 

phase- coherent systems. This is so because simplex signals such as phase reversal 
keying can be employed. For larger alphabets, direct detection with proper prefiltering 
appears to be superior (see Appendix C), but such systems have yet to be implemented. 
The performance of large- alphabet direct- detection systems using wideband filtering can 
be either superior to or inferior to heterodyning and homodyning systems, depending 
upon the background noise level and predetection bandwidth. It is important to note that 
these statements are based upon isolated examples. Much more work must be done before 
the relative merits of different systems can be fully determined — even for the free- 
space channel. 

3.3 LARGE- APERTURE RECEIVERS 

For a given collecting area, the optimum receiver for the free- space channel uses 

diffraction- limited optics. On the other hand, the available and contemplated size of 

41 

diffraction- limited optics is relatively small, so it is possible that improved system 

performance could be achieved by employing larger nondiffraction- limited optics. 

There are several ways in which such a system could be implemented. If, for 
example, direct detection were employed, the system could differ from a diffraction- 
limited system only in the use of a larger detection area, because of the larger point 
spread function of the optics. Ctf course-, this larger area will result in an increased 
noise level. A more refined direct-detection system might employ predetection spa- 
tial filtering to compensate for some of the distortion introduced by the optics. The 
objective of this compensation would be to reject spatial modes that contain only noise. 

Heterodyning with nondiffraction-limited optics' is more difficult than direct detec- 
tion because, as noted in section 3. 2. 1, heterodyning systems respond to a single, or 
at most to a few, spatial modes of the received field. Thus a nondiffraction- limited 
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heterodyning system will not extract all of the available energy in the aperture, unless 
the local- oscillator field is properly matched to the signal field that exists at the mixing 
plane. This matching requires that either the phase distortion or "aberrations" intro- 
duced by the optics be known to, or measured by, the receiver, or that the local oscil- 
lator transverse that portion of the system that is not diffraction- limited. 

Although it may be possible to achieve the required matching by compensation of 
either the signal field or local- oscillator field, it would appear to be more difficult to 
implement such a system than it would be to implement the direct- detection system. 
Under the assumption that both systems can be implemented, their relative performance 
depends upon the background noise level and detector internal noise. 

Transmitted reference systems provide another approach to operation with large 
apertures that are not diffraction- limited. Such systems can be viewed either as 
heterodyning systems in which the local oscillator is located at the transmitter, or as 
direct- detection systems using a particular signaling scheme. For this discussion the 
former view is more useful, since it clearly demonstrates that the "local- oscillator " 
field passes through that portion of the system that is not diffraction- limited. One dis- 
advantage of a transmitted reference system is that a strong "local-oscillator" signal 
is obtained at the expense of signal power. As indicated in Fig. 2, the loss of signal 
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Fig. 2. Heterodyne reception with a transmitted reference. 


power can be reduced, at least in principle, by amplifying the reference field while pre- 

1 3 

serving the phase of the spatial modes. The reduction is gained at the expense of 
increased system complexity and a noisier local-oscillator signal. 

The feasibility, relative merits, and ultimate limitations of all three of the schemes 
above should be explored. At this junction the only firm conclusion is that large- aperture 
receivers with nondiffraction- limited optics should be given careful attention.lt is worth 
noting that this conclusion will be reached for all of the channels that we consider. 

3.4 QUANTUM COMMUNICATION THEORY 

It is apparent from the preceding discussion that quantum effects comprise an impor- 
tant component of the system "noise, " and that the nature of this noise depends in detail 
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upon the system structure. Although this "noise" can be treated by semiclassical means 

once the system structure has been specified, fundamental difficulties are encountered 

in determining the optimum system structure. These difficulties are the subject of quan- 

38-40 42-47 

turn communication theory. ’ We shall now describe the nature of these dif- 

ficulties, briefly review the approaches that have been explored, and summarize the 
outstanding unresolved questions. In general, our remarks are applicable to all optical 
channels, not just the free- space channel. 

It is clear from the discussion of heterodyning and direct- detecting systems that 
optical receivers are subject to a noise whose properties are determined by the 
receiver's structure. Stated alternatively, the receiver interacts with the channel in 
a way determined by the detailed nature of the receiver. This interaction poses no 
fundamental conceptual problems in the analysis of any specific system, but it does raise 
deep issues in system optimization. Specifically, over what class of receivers is the 
optimization to be performed ? Clearly, we prefer to optimize over the class of all 
receivers that are consistent with the laws of physics. The problem then becomes that 
of characterizing the class of such receivers in a manner that is amenable to anal- 
ysis. 

Two approaches to this problem are now available; both are based upon the axioms 
of quantum theory. The first, and perhaps more traditional, begins with the premise 
that every receiver can be characterized by a set of "observables" and that the per- 
formance of the receiver can be determined according to the techniques of quantum the- 
ory (see Appendices D and E). This approach has been rewarding, but it does have 
some shortcomings. In particular, it appears that receivers exist that cannot be 
characterized by a set of observables, for example, a system employing optical pream- 
plification. 

A second approach to receiver characterization has recently been introduced 
(Appendix F). It consists, in essence, of relaxing some of the constraints that an 
"observable" must satisfy. Precisely stated, the first approach characterizes a meat 
surement by a complete set of orthogonal vectors; in the second approach the vectors 
need only be complete. The second approach does encompass optical preamplification 
and, perhaps, may encompass all receivers that are consistent with the laws 
of physics. 

Both receiver characterizations are useful, but fundamental problems are encountered 
in their application. For simplicity, we discuss these problems in the context of the 
first approach; that is, we assume that every receiver of interest can be described by 
a set of "observables." 

First, it is doubtful that every "observable" corresponds to a receiver — even ideally. 
Thus, if one attempts to find the best receiver in the class of those characterized by 
observables, he may be led analytically to an observable that does not characterize a 
receiver. This situation is further complicated by lack of any synthesis procedure. That 
is, there is, at present, no way of determining how a given observable might be 
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"measured." Instead, it is necessary to recognize the observable as that associated 
with a known measurement. These problems are fundamental and must be resolved 
before the performance of quantum- limited communication systems can be fully 
assessed. 

In addition to the conceptual problems just discussed, there are many difficult and 
unresolved problems in the analysis of quantum systems. For example, even for 
binary signaling, the optimum system performance has been determined only for a 
very limited number of waveform sets and channels. For larger alphabets little prog- 
ress has been made, except for a special case that corresponds roughly to orthogonal 
signaling alphabets (see Appendices C, D, and F). 

Although the available results are fragmentary, they do suggest some interesting con- 
clusions and questions. First, in some specific situations the performance associated 
with the optimum observable is superior to that attained by the heterodyning, homo- 
dyning, and direct- detecting systems usually contemplated. Moreover, it appears that 
the optimum observable can, at least in some cases, be implemented in principle 
(see Appendix F). On the other hand, in the problems that have been solved, this supe- 
riority in performance has been rather small — so small as to be of little practical 
importance. 

It is apparent that our understanding of quantum systems is quite meager. This is 
due in part to the limited attention that such systems have received, and in part to the 
conceptual and analytical difficulties involved. In any event, the study of such systems 
represents an important aspect of optical communication and should be continued. In 
particular, an effort should be made to determine whether or not there are systems 
whose performance is substantially superior to that which can be realized by a combina- 
tion of optical amplifiers, local oscillators, and ideal photon counters. 



IV. Atmospheric Channel 


We shall consider optical communication systems involving propagation through the 
clear turbulent atmosphere. The term "clear atmosphere" emphasizes that only the 
effects that are due to random refractive index variations are being considered. Haze, 
smog, clouds, and other atmospheric conditions characterized by concentrations of 
aerosols causing scattering effects are considered in Section V. 

The effects of interest to us are caused by the random variations of the optical 

refractive index in time and space. These refractive index changes are caused by small 

temperature variations in the atmosphere, which in turn are caused by turbulent mixing 

of various thermal layers. Although the typical temperature variations are only 1°K, 

the resulting refractive index variations of a few parts in 10^ cause extensive effects on 

optical radiation fields propagating through the atmosphere. The astronomer viewing a 

star has long been faced with image dancing or wander, image distortion (tremor disk), 
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and scintillation. The image dancing, observed in small telescopes, corresponds to 
variations in the angle of arrival of the energy. The tremor disk, observed in larger 
telescopes, corresponds to the deviations from a plane phase front caused by refrac- 
tive index variations. Scintillation is the term used to describe the variations in inten- 
sity of the image. Thus the random refractive index variations make the received field 
fade in both time and space and also cause extensive phase front distortion. 

At first glance, these effects might appear to severely limit communication through 
the atmosphere at optical frequencies. Indeed, some investigations have shown that sys- 
tems that are near optimum in the absence of turbulence suffer large performance degra- 
dations in its presence. On the other hand, experience at other frequencies and 
preliminary results at optical frequencies suggest that atmospheric turbulence may not 
severely limit either rate or reliability if the system is properly designed. One of the 
major objectives at the Workshop was to determine the issues that must be examined 
before the potential for atmospheric optical communication can be fully assessed. All 
of these issues revolve around the ideas of modern communication theory. 

A prerequisite for the application of the ideas of modern communication theory is a 
statistical description of the effects caused by the channel. For a clear atmosphere at 
optical frequencies, a complete statistical description of the received field, given the 
transmitted field, is required. A goal of the group studying the turbulent atmospheric 
channel was to assess the present state of theoretical and experimental knowledge related 
to this description to establish the areas in which further work would be required before 
an adequate model could be obtained. A second goal was to use the available channel 
descriptions to assess avenues of system development that might ultimately lead to 
efficient systems. 

We shall review the effects of atmospheric turbulence upon optical radiation and the 
areas in which further knowledge is required. Then we shall discuss a channel model 
that is consistent with available knowledge . Finally, we shall use this model to explore 


21 


22 


OPTICAL SPACE COMMUNICATION 


the avenues along which efficient systems might be developed. The emphasis will be on 
a ground-to- space or space-to-ground path, but the results for terrestrial paths will be 
used whenever they are helpful. 

4.1 TURBULENT EFFECTS 

We turn now to a detailed examination of the effects caused by atmospheric turbu- 
lence, that is, by the random refractive index variations. 

4. 1. 1 Theoretical Approaches 

Theoretically, the turbulent irregularities in the atmosphere should cause some depo- 
larization of a linearly polarized light wave. The effect has been calculated and esti- 
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mated to be approximately 10 per km. Saleh was unable to observe the depolarizing 
effect of the atmosphere on a 2.6-km path, even with an equipment sensitivity of -42 dB. 
We conclude, at the present time, that depolarization effects of the atmosphere are of 
no practical interest in the design of optical communication systems. Consequently, the- 
oretical studies of turbulent effects reduce to solving the wave equation, with the refrac- 
tive index described as a random field in time and space. 

50 

By the use of a frozen atmosphere assumption and Fourier transform techniques, 
the problem is further reduced to solving the Helmholtz equation 

V 2 U(?) + k 2 [n(?)f U(?) = 0, (6) 

where U(r) is the complex amplitude of one frequency component of the field, k = 2tr/k, 
and n(r) is the random refractive index at an instant of time. The time behavior is intro- 
duced by assuming that the refractive index field, which is frozen in space, is blown 
across the optical path by a normal wind component. The temporal behavior is thus 
specified in terms of the spatial behavior that is already known. The validity of the 

frozen atmosphere assumption depends on the fact that the refractive index changes are 
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very slow relative to the motion induced by the wind. 

Several approximate solution techniques have been tried on Eq. (6) and all have limi- 
tations . One of the most popular and controversial of these is the Rytov solution tech- 
nique used extensively by Tatarski,^ Chernov, and others. This approximation has 
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recently been under severe attack, and must certainly be modified for very long 

and/or very turbulent paths to account for certain experimental facts. A modification 
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has already been proposed by Tatarski, and others" are under active development. 
While the Rytov approximation technique has many critics and some proven failings, it 
does have some important virtues. The most important of these is that it has been in 
agreement with most of the available experimental evidence. Furthermore, while it 
is known to break down for very turbulent or long paths, it should be applicable for 
almost any ground-space path. Moreover, it is by far the most satisfactory theory 
now available. 

In the Rytov technique it is assumed that the solution to Eq. (6) is of the form 
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U(r) = U Q (r) exp y(r), (7) 

where LMr) is the solution when n(r) = n( r) and the bar indicates an expectation opera- 
tion, 

y(r) = x(r) + j<i>(r), (8) 

with x( r ) the normalized log amplitude, and 4>(r) the phase difference between the unper- 
turbed solution and the actual solution. The technique involves substituting Eq. (7) in 
Eq. (6) and linearizing the resulting equation for y(r), which can then be solved. The 
importance of y(r) in the theory has made it occupy a central role in experimental work. 

To study y(r), it is necessary to have a description of the random refractive index 
n(r). This description is now considered to establish some needed notation. In terms 
of N = (n-1) 10 6 , the deviation of the refractive index from unity in parts per million, 

the refractive index fluctuations are proportional to the temperature fluctuations, with 
57 

the relationship 


6N = 


79 

V - 1 



(9) 


where p is the pressure in mb, T is the temperature in °K, and V = c /c , the ratio 

IT V 

of specific heats, is approximately equal to 1.4 for air. Variations in the mean of N(r) 
with r correspond to large-scale or systematic refractive effects that are not being 
considered here. It follows from the Kolmogorov theory of turbulent flow at large 
Reynolds numbers, particularly as applied to temperature fluctuations by Obukhov 
and Corrsin, that the second-moment description of the spatial variations of N is given 
by the relation 58-80 


D N (r) = C N 


iri 2/3 


( 10 ) 


for r < L 


o 


where L is the outer scale of the turbulence and is of the order of meters, 
o 


Here D-^(r) is the refractive index structure function defined by 


D N (r) = [N(R+r)-N(R);r, 


( 11 ) 


where r = |r | and the bar indicates an expectation operation over the ensemble of atmo- 
2 

spheres. C, T is the refractive -index structure constant that depends strongly on mete- 
™ . 2 
orological conditions, time of day, and altitude, among other things. C M has been 

6163 ^ ' 2 

measured for various conditions on horizontal paths. The variation of with 

altitude has been discussed, 8 ^ but actual measurements are badly needed, as we shall 
discuss in sections 4. 1.2 and 4. 1.4. 

We shall see that this second-moment characterization of the refractive index vari- 
ations appears to be sufficient to answer most questions of interest. A less adequate 
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second-moment model that describes the spatial variation of N(r)in terms of a Gaussian- 

shaped correlation function has also been used, for mathematical simplicity, in some 
51 

theoretical studies. Certain detailed questions involved with the time evolution of the 
effects caused by the refractive index variations may require a more complete statistical 
model than the second-moment characterization now being considered. 

After this short survey of the over- all theory, we now turn to a more complete 
examination of the theoretical and experimental results associated with various aspects 
of the received field. 

4.1.2 Intensity Scintillation 

Intensity scintillation caused by turbulence has long been noted by astronomers, and 
it is the most noticeable feature of laser radiation received over an atmosphere path. 
The turbulence causes the received field to fade in both time and space; that is, the 
beam exhibits a pattern of light and dark spots that change with time. To first order, 
the pattern seems to move across the beam. 

To examine the effect of this phenomenon on optical communication systems, it is 
necessary to characterize the random amplitude variations conditioned on the values of 
the transmitted field. This characterization is considerably simplified if the frequency 
coherence bandwidth of the turbulent effects exceeds the bandwidth of the signal that is 
being considered, that is, if the amplitude and phase of each frequency in the trans- 
mitted signal undergoes the same turbulent change. This appears to be a reasonably 
good assumption for most signal bandwidths. While there are not direct measurements, 
both theoretical and experimental evidence indicate that the envelope and phase vari- 
ations caused by the turbulence are highly correlated over large bandwidths for ground- 
space paths. 5 Stellar scintillation data indicate selective color fading over bands of 

o o o A4 

1000 A and essentially identical intensity traces for 200 A bands separated by 833 A. 
Pulse transmission studies with a GaAs laser over a 1 . 8-km path have shown that the 

q 6 8 

coherence bandwidth is greater than 10 7 Hz. 

With the frequency coherence assumption, it is only necessary to characterize the 
scintillations that result for a transmitted field containing a single frequency. These 
scintillations obviously depend on the spatial character of the transmitted field; that is, 
whether the field is a plane wave, a spherical wave, or a collimated or focused beam. 

Statistics : Both theoretical arguments and measurements support the contention that 
conditionally the received field is a log-normal process. For Gaussian beam and plane 
or spherical wave propagation with reasonable path lengths, it is easy under the 
Rytov approximation to use the Central Limit theorem to show that at any instant of 
time the set of x(r)'s associated with any set of r's (at least those in the central portion 
of the beam for that case) are jointly Gaussian random variables. Here x(r> is the nor- 
malized log amplitude discussed previously. If this result is extrapolated to claim 
that x(r,t), the quantity that describes the turbulent effects on the field envelope, is a 
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Gaussian random process of time and space, then the envelope of the received field 
is a log-normal process. 

Experimental measurements of the first-order density of the intensity of the received 
field have been made over horizontal paths in both the spherical wave and beam cases 
and they tend to confirm the log-normal result. Note that a log-normal envelope 
distribution implies a log-normal intensity distribution. There is some evidence that 

69 

the intensity distribution is not log-normal for very turbulent situations (or long paths). 

For collimated beams, measurements have tended to be concentrated at the center of 

the beam. There is some evidence that near the edges of the beam the field amplitude 

6 8 

variations may have a more Rayleigh or Rician character. This may be due to the fact 
that scattering processes tend to cause the dominant effects in these weak intensity 
regions. While the intensity statistics near the beam edge may be of interest during the 
acquisition phase, they will generally not be important after the terminals are locked in 
and tracking one another. 

The only intensity statistics measurements associated with a ground-space path are 
those from GEOS- II. These data are currently being reduced, but preliminary results 
seem to be in agreement with the log-normal hypothesis. There have been no measure- 
ments of the higher order densities in an attempt to verify or disprove the log-normal 
process assumption. Also, most measurements have been made with small- aperture 
receivers. As the receiver diameter is enlarged to include many regions with inde- 
pendent fading, one would expect the resulting output signal to show Gaussian fading. 

62 VO 

There is some experimental evidence and also a theoretical argument to indicate 
that this effect will either not be observed or will occur only with exceedingly large 
apertures. It is worth verifying this experimentally as large collectors become avail- 
able. Also, while one would not expect any effect, measurements of the intensity sta- 
tistics with large transmitting apertures would be desirable. 

If we accept the premise that the received intensity is a log-normal random process, 
its statistical description reduces to a specification of its mean and covariance. This 
specification requires only a second-moment description of the refractive index varia- 
tions. Even if the intensity is not log-normal, the mean and covariance functions are 
important in the design of efficient communication systems. Thus we now consider them 
in some detail. 


Mean Intensity and Average Beam Divergence: The mean intensity for plane- and 

spherical- wave propagation is normally calculated by using a conservation- of - energy 

71 — 

argument. This amounts to an assumption that I = I q ; that is, the average intensity 
is equal to the intensity in the absence of turbulence. We shall show that this assump- 


tion, when coupled with the assumption of log-normal statistics, leads to the result that 



i.e., the mean of the normalized log amplitude is the negative of its variance. 


In optical-communication applications, we are often concerned with the propagation 


of well- collimated (perhaps even focused) beams. The phase fluctuations induced by tur- 


bulent irregularities in the atmosphere cause portions of the beam to propagate in 
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slightly different directions and hence produce a spreading of the beam with a resultant 
loss in average antenna gain. The turbulence also causes small changes in the average 
direction of propagation of the beam, with the result that in the receiving plane the cen- 
ter of the antenna pattern is not fixed in space but rather moves about in time. Thus, 
for beam propagation, both of these effects will be reflected in the means and vari- 
ances associated with the intensity profile in the receiving plane. The mean intensity 
profile can be considered to be the average antenna beam pattern. The intensity vari- 
ance reflects variations from this average beam pattern, because of beam motion and 
spreading, as well as scintillation. 

A rigorous analysis of the intensity mean and covariance functions associated with 

beam propagation is quite difficult. All theoretical work in the past has been confined 

to Gaussian beams. Beckmann has presented an approximate geometrical approach to 

72 

the beam propagation problem. The Rytov technique has been used by Kon and 
73 74 

Tatar ski and by Schmeltzer to study Gaussian beam propagation in the turbulent 
75 

atmosphere. Gebhart has recently used Schmeltzer's results to study the mean 
intensity. This study is limited, however, to points near the nominal position of the 
beam center. A technique suggested during the study, relying on atmospheric reciproc- 
ity, transforms the calculation of the mean intensity for the case of a uniformly 

illuminated transmitter aperture into the calculation of the mean signal power associated 

T 6 

with an equivalent heterodyne receiver (see Appendices G and H). 

The ratio G(r) of the average antenna gain in the atmosphere to that with no atmo- 
sphere (the normalized mean intensity profile) can be written in terms of the mean and 
variance of x(r), the normalized log amplitude of Eq. (8), as 


G(r) 


G , (r) 

atmos 


G. , ,(r) 

ideal' ' 



(r)+x(r) 


)] 


( 12 ) 


The result depends on the assumption that x is a Gaussian random variable. Even 

granting this assumption, the available numerical results for x(r) are only for horizon- 

— 74 7 5 

tal (homogeneous) paths and r approximately at the center of the beam. ’ Better 

2 

data on the variation of the atmospheric structure constant, C^, with altitude are needed 

before it is worth while to try to extend these numerical techniques to ground-to- space 

paths. Even more fundamental, the x calculation that is being discussed here applies 

only to Gaussian beam shapes and is not easily extended. 

A study of the average beamwidth for collimated Gaussian beam propagation, which 

7 3 

is valid in the near field of the transmitter aperture, has been carried out. The pre- 
dicted beam radius is 

R = constant C.. T £*/ 3 (kaln 0)*^ L (13) 

eit JN o 


where 
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2 2 
, 2L N 
n 2 m 

9 - "72 2 
k a 


N = 5.91 l 
m o 


L = path length 

a/2 = standard deviation of the Gaussian intensity distribution 

i = inner scale of turbulence . 
o 

Note that R grows slightly more rapidly than linearly with path length, because of the 

1 /2 en 
(In 6) ' term. 

There are only a few experimental measurements associated with average beam 

spread and mean intensity. The separation of beam motion into beam wander and beam 

spread and the separation of total intensity fluctuations into variations of the mean and 

deviations from the mean are obviously arbitrary. The separation has normally been 

accomplished in experimental situations by an arbitrary division according to time 
77 

scale. Ochs, at ESSA Research Laboratory, has measured beam wander, beam 
spread, and mean intensity over a 15-km path. An array of 4 photomultipliers was 
used to provide beam- wander information to the transmitter. An integration time 
of 0.5 sec is used in determining beam wander. Every 15 min the vertical-beam 
profile (beam spread) is measured by sweeping the beam past a fifth photomulti- 
plier at a rate of 1.9 prad/sec and recording the resulting intensity with a 0. 5-sec time 
constant. The beam spread is measured as the time interval between the 3 dB points 
divided by the angular sweep velocity. Finally, the peak average intensity is recorded. 
The results of these measurements are shown in Fig. 3. The values are seen to depend 
heavily on the time of day. Also, the mean intensity during daylight hours does not vary 
as the reciprocal of the square of the beamwidth. This may be due to increased extinc- 
tion by haze, asymmetry of the beam, or hot spots (concentrations of energy within 
the beam). It would be of value to know the actual cause. 

Other beam-pattern experiments have been conducted at ESSA Research Laboratory 

78 

over horizontal terrestrial paths. These experiments at the HeNe wavelength include 
both fast- and long- exposure photographs of the received beam pattern over paths 
varying from 550 m to 145 km. The long-exposure photographs show that the atmo- 
sphere limits the extent to which a beam can be focused and causes the half-power beam- 
width to grow faster than linearly with path length. In particular, transmitter apertures 
larger than 11 cm caused no further focusing effect, and the beam diameter followed 
the relation 


D = alf. 


(14) 
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15 KILOMETER PATH 
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Fig. 3. Measurements of beam wander, beam spread, and mean 

77 

intensity for a 15-km path. (Measured by Ochs. ) 


where a and (3 are 4.5 X 10 6 and 1.2, respectively. The results here are in rough 

73 

agreement with the predictions of Eq. (13). 

The only average antenna gain data available for laser transmission over a ground- 
space path is that from the Surveyor YII laser pointing experiment. These results seem 
to imply that the beam is broadened by a factor of two or three by atmospheric turbu- 
lence, but this is at best a crude estimate. 

To the extent that the atmosphere is reciprocal, it is reasonable to expect that the 
average size of star images observed with a terrestrial telescope can be used to esti- 
mate the average spreading of a beam transmitted from the Earth along the same slant 
path. The image size in the focal plane of the telescope is a measure of the angular 
spectrum required to represent the distorted wavefront of the stellar energy at the 
telescope objective. This fact, which will be considered again when phase fluctuations 
are considered, seems to imply that a reasonable estimate of average beam spread can 
be obtained by convolving the transmitter beamwidth with the average stellar angular 
spectrum. Image diameters of stars observed with large telescopes are usually within 
the range of 5-20 prad at night. There are no corresponding data for the daytime hours. 

It is clear that more experimental data on average beam spreading and mean 
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intensity for ground -to -space paths are needed. A possibility for obtaining more data 
is to use the GEOS- II Argon detector that is now in orbit. Because of the pointing 
accuracy that is possible, the frequency limitation of the detector at 50 Hz, and the low 
orbit with the resulting limited observation time, the following kind of experiment 
looks attractive. An aperture mask would be used to form a pattern of narrow fan beams. 
These beams would then be moved through space at such a rate that the satellite passes 
through the pattern slowly, thereby causing the detector to scan the beam pattern. The 
measurements of the beam pattern and mean intensity at the center of the fan beams 
could then be compared with the theoretical values in the absence of turbulence. If pos- 
sible, simultaneous measurements should be made of the angular spread of stars in the 
vicinity of the satellite to test the hypothesis that beam- spread information can be 
obtained from stellar observations. 

Intensity Variance : Note that with a solution of the form of Eq. (7) the turbulence 
effects on the field envelope (or intensity) are multiplicative, in the sense that the unper- 

v/j* •£ ) 

turbed field envelope is multiplied by e M ’ , where x(h, t) represents the effect of tur- 

bulence. Under the assumption that the Rytov approximation is valid, the variance of x 

has been calculated in the plane- wave, spherical- wave, and Gaussian-beam cases. In 

79 

the plane-wave case the result is 


2 . R , ,7/6 

o- = 0 . 56 k ' 

X 


f 

J 0 


C* (s) s 5/6 ds, 


(15) 


80 

while the result of the spherical- wave or collimated Gaussian-beam cases is 


//tv r<-5/6)k 2 Re 


J ds C^(s) 1 Re g(z, s) 


5/6 


- — g(z, s) 


5/6 


, (16a) 


where C^s) is the refractive- index structure constant defined in Eq. (10), z is the path 
length, is the standard deviation of the amplitude of the Gaussian beam, and 


2(z-s) s - ike 


g(z,s)= — 


z - ika 


2 ‘ 


(16b) 


The integrations are along the optical path. The result is applicable at positions near 
the beam center and reduces to the spherical- wave case when a — 0. Indeed, for 
Gaussian- shaped laser-beam propagation over path lengths such that kc^/z « 1, the 
resulting variance is in good agreement with the spherical- wave result, which can 
be written for a homogeneous path as 8t ^ 

0.124 


(17) 
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These equations can be used to calculate answers, once C, T (s)is specified along the opti- 

1N 2 

cal path. These calculations have been carried out for specific choices of C M (s); in par- 

79 80 

ticular, results are available for a homogeneous path, ’ and for a particular choice 

81 82 

of atmospheric vertical profile. ’ The numerical values of the atmospheric profile 

results are open to question, because of the lack of a verified atmospheric model and 

recent changes in the postulated model. ^ In particular, the dependence for plane- 

and spherical- wave propagation, ^ and the dependence for plane- and spherical- 

8 3 

wave propagation over short enough homogeneous paths seem to be confirmed. The 
wavelength dependence, of course, is important in frequency selection. 

Experimental results over very long and/or very turbulent atmospheric paths point 
out the inadequacy of the Rytov approximation results for these cases. Several experi- 
menters have noticed that when the turbulence is severe enough and the path sufficiently 

long, the log-amplitude variance saturates and remains constant as the path length is 
83 85 86 

increased. ’ ’ Tatar ski has proposed a modification of the theory to account for 

56 

this effect. While it is important to better understand this saturation effect and its 
implications, it is unlikely to occur on any reasonable optical path associated with a 
ground-space link. To obtain an expression for the intensity variance from the log- 
amplitude variance, it is necessary to use the log-normal statistic assumption and to 
know the mean of x. The result is 


2 7 2 
*1 =1 


4<r 

e X -1 


I = I e 
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2 

X+o - 

X 


(18) 


(19) 


with I the intensity in the absence of turbulent effects. The value of x has already 

been discussed. As indicated, for plane- or spherical- wave propagation a conservation- 

— _ g 

of-energy argument can be used to argue that I = I , with the result that X = -o' . 

u X. 

Stellar measurements have provided some experimental data on the intensity vari- 
ance at visible wavelengths during the night. In particular, measured values of a^/T 

(sometimes called the coefficient of intensity variance) for starlight are approximately 

48 

0.4 for small apertures and zenith paths. This corresponds to of approximately 

0.2. Intensity scintillation measurements have been made, using a detector on GEOS- II 

with an Argon laser on the ground. The transmitter has a beam divergence of 0.7 mrad, 

and the detector aperture has a 1-cm diameter. The value of <r /T measured over a 

1 84 

160-Hz bandwidth was 0. 3, which is roughly in agreement with stellar observations. 

The intensity (or log- amplitude) variance depends on both the transmitting and 
receiving apertures. The variance expressions discussed above were for receiving 
apertures that are small relative to the intensity correlation distance. As the 
receiving aperture becomes larger, it includes several independent intensity values and 
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a law of large numbers argument can be used to show that the variance decreases. For 

the case of plane-wave propagation along either a horizontal homogeneous path or a path 

from space to earth using a particular atmospheric model the amount of the reduction 
71 

has been calculated. The results show that for reasonable log- amplitude variance val- 
ues (<0.6), the intensity variance decreases by a factor equal to the coherence area 
divided by the receiver area, where the coherence area has a diameter approximately 
equal to the coherence diameter defined in the next section. 

There are a lot of experimental data on the aperture averaging effect for stellar 
48 87 

sources. ’ These data show reasonable agreement with theory. Experiments also 

6 2 

have been carried out with laser sources and terrestrial paths. These measurements 
seem to indicate that a beam source has different correlation properties from those 
attributed to a plane wave and that the received intensity may have larger correlation 
distances and hence less aperture averaging. These correlation properties will be 
discussed. 

The variation of intensity variance with transmitter aperture is more complicated. 
Here, we are concerned with beam transmission and a very small (relative to the cor- 
relation distance) receiver aperture. It appears that the intensity variance will depend 
on the aperture illumination, as well as on its size. In the case of a uniformly illumi- 
nated aperture the assumption of atmospheric reciprocity permits the calculation of the 
received-field variance in terms of the output of a heterodyne detector whose aperture 

area is the same as that of the transmitting aperture area. The variance of the output 

88 

of a heterodyne detector has been studied. In the case of a Gaussian- shaped aper- 

74 

ture illumination Schmeltzer's beam- propagation results have been used to examine 

the effect of beam size on the intensity scintillation. This has been done for focused and 

collimated beams on homogeneous paths and for collimated beams on an Earth- to- Space 
80 81 

path. ’ The results tend to show that the resulting log-normal variance can be 

written as a product of the spherical-wave variance and an aperture- dependent factor 

that is a function of ka Q /z. For a homogeneous path, the aperture -dependent factor does 

not vary too much for collimated beams, but decreases rapidly for focused beams when 
2 

k a /z becomes larger than one. For the vertical transmission path and a collimated 

o 2 

beam, the aperture-dependent factor decreases with increasing values of ka Q /z, and for 
large values goes approximately as the inverse of the aperture area. 

There appear to be no detailed experimental measurements of the transmitter aper- 
ture effect. The theoretical variation of scintillation for a collimated Gaussian beam 
and a vertical path suggest that an interesting experiment would be to vary the size of 
the beam that is used with a ground-based transmitter and monitor the scintillations 
at a spacecraft-borne receiver. 

Intensity Correlation: As we have indicated, the received field, after propagation 
through the turbulent atmosphere, exhibits a spatially varying intensity. A particularly 
important parameter of this variation is its correlation distance, the distance over which 
the field envelope is correlated. It is possible to define this parameter in many 
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essentially equivalent ways; here, we shall define it as the value of separation for which 
the spatial covariance function of the envelope goes to zero for the first time . 

Again, most of the theoretical results are in terms of the log- amplitude of the field. 
The assumption that this quantity is a Gaussian process means that lack of correlation 
implies statistical independence. This is a particularly important feature that will 
be exploited in the design of some of the receivers discussed in section 4. 3. 

The spatial covariance function of the normalized amplitude, C (p), has been studied 

A. 

by using the Rytov approximation. Note that in general the spatial covariance function 
would have to be written 

C x (rrr 2 ) = [x(r 1 )-x(r 1 )][x(r 2 )-x(r 2 )]. (20) 

For many situations of interest, however, the covariance across the receiver plane can 
be written in terms of the single scalar 



The intensity covariance function can be written in terms of C (p) by using the log- 

A. 

nromal property. It has the form 

Cj(p) =I 2 {exp[4C x (p)]-l} (21) 


when the field is homogeneous over the set of p values of interest. Note that, while the 
intensity and log- amplitude correlation distances are equal, the curves have greatly dif- 
ferent shapes, with the shape of the intensity curve dependent on the amount of turbu- 
lence present, that is, the log-amplitude variance. 

In the case of plane- and spherical- wave propagation over homogeneous paths, it is 
possible to write C (p) in the form 

A. 

C X ( P ) = Af(p/(4z/k) l/2 ), (22) 

where f is a function that depends on the propagation mode but is insensitive to the 

amount of turbulence. The first zero of f is defined as the correlation distance, L^, 

and is approximately 0 . 76V XL for the plane-wave case, and 1 . 8 n / XL for the spherical- 
81 82 89 

wave case. ’ ’ Here L is the propagation path length. 

In the beam case, an expression like Eq. (22) is no longer sufficient. A rather com- 
plicated integral expression for C (p) is available for the case of Gaussian-beam propa- 

X. 

gation (and hence spherical- wave propagation in the limit of zero source width) over 
74 

any path. While this expression is difficult to evaluate, it has been suggested that an 
expression of the functional form 

c x ( p ) = °'^g 1 (p/(4zA) 1//2 , p/(2z/ka Q )) = o- x g 2 ^p/(4z/k) 1 ' /2 , fi). 


(23) 
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2 . 82 

where G = ka Q /z, may be appropriate in this case. Note that this expression involves 
the additional scale length provided by aperture diffraction effects. 

Experimental measurements over terrestrial paths tend to confirm these 
results.^ 0, When £2 » 1, even the beam propagation results are in reasonable 

agreement with the plane- and spherical-wave theory. When Q becomes smaller and 
less than one, however, there seems to be the possibility of correlation over larger dis- 
tances of the order of 2z/k a . There is some experimental support for this in terms 

o 

of studies of the intensity variation with aperture diameter. The dependence of the 
correlation distance on the source size is an important consideration for terrestrial opti- 
cal communication systems. For a ground- space link it is probably less important 
because the satellite receiver will normally intercept energy over an area that contains 
at most a few, and normally only one, coherence areas, that is, areas over which the 
intensity is essentially constant. Also, the value of Q associated with the down link will 
normally be large. Neither theoretical nor experimental results are available for the 
log- amplitude covariance functions on paths such that the intensity has saturated. 

If the turbulence is not homogeneous along the propagation path, the results are 

changed. If the variation is a smoothly decreasing (or increasing) function along the 

2 

path, as in some models for the variation of with altitude, then the major modifica- 
tion to the results will be to replace the path length with a scale height associated with 

2 2 
the variation of C^. Thus, for the assumption of an exponential variation of with 

altitude, the correlation distance depends on the parameter (4h sec 0/k) 1 ^ 2 , where 

o 80 82 

h is the scale height of the atmosphere, and 0 is the zenith angle of the path. ’ On 

the other hand, anomalies in the path will tend to further distort the covariance function 

from its homogeneous path form. For spherical- wave or Gaussian beam propagation, 

the effects of various turbulence regions in changing the shape of the covariance function 

depend on their position in the path. It has been suggested that this dependence of C (p) 

2 

on the turbulence profile along the path be used for indirect measurements of C N along 

91 iN 

the path. This would involve measuring C (p) and then trying to invert an appropri- 

X 74 

ately modified form of an equation of Schmeltzer. Further investigation would be nec- 
essary to establish the feasibility of carrying out this inversion procedure on this equa- 
tion. Such a procedure using stellar or laser sources would have value in monitoring 
the amount of turbulence at the tropopause, and might even have some value in the detec- 
tion of clear air turbulence if it were sensitive enough. 

Measurements along inhomogeneous terrestrial paths have shown the effects 
described above. Measurements between two mountain peaks, which emphasized the 

turbulence at the two ends of the path, showed correlation distances corresponding 

8 5 

to the positions of the two effective turbulence regions. About the only data avail- 
able on spatial intensity correlation for ground- space paths are stellar data. Both direct 
and indirect measurements seem to indicate that the correlation distance is between 7 cm 
and 10 cm for a zenith path. ’ Indirect results tend to indicate that the correlation 
distance is proportional to the square root of the secant of the zenith angle. 
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Frequency Spectrum of Sc intillation: If the frozen- atmosphere assumption that was 
discussed at the beginning of section 4.1.1 is assumed valid, the spectrum of the log- 
amplitude (or intensity) fluctuations can be obtained from the spatial covariance function. 
Indeed, a joint time- space covariance function can be defined as 

— „ ( v-p \ 

C x (p, t) A[x(r 1 ,t 1 )-x(r 1 ,t 1 )][x(r 2 +p,t 1 +T)-x(r 2 +p,t 1 +T)] = C I p + — — T j , (24) 

where v is the wind velocity vector normal to the optical path. The spectrum of x is 
obtained by setting p equal to zero in Eq. (24), that is, making the argument vt, and then 
Fourier transforming. This has been carried out by Tatarski for plane-wave propaga- 
tion over homogeneous paths and he has obtained a normalized form that is only a func- 
tion of the variable f*J VL/v, where v is the normal wind velocity. 50 Thi^ result has 

50 93 

been tested experimentally by Tatarski and Ryznar. ’ 

Intensity scintillation spectra for laser sources and terrestrial paths have been mea- 
94 95 68 93 

sured by many people. ’ ’ ’ These measurements show that almost all of the 

power is confined to the first few hundred Hertz. These spectra seem to be of the 

form Sj(f ) = A e a ^, where a depends on the atmospheric conditions, for example, wind 

speed. The spectrum seems to be independent of receiving aperture for collimated 
9 5 

beam propagation. 

48 87 

Many scintillation spectrum measurements have been made for stellar sources. ’ 

In general, for stars near zenith and with small apertures, the spectrum exhibits a break 
point at approximately 100 Hz and decreases to essentially zero by 500-1000 Hz. The 
spectra are receiving aperture sensitive, with the high frequencies decreased as the 
aperture increases. ^ This is due to more averaging of the rapid spatial changes 
that are responsible for the higher frequencies. 

Spectral measurements made with a detector on GEOS- II have shown a peak in the 
spectrum near 100 Hz. This may be caused by the relative spacecraft- ground motion 
acting as a pseudowind. The measurements in general confirm the lowpass nature of 
the scintillation spectrum. 

4.1.3 Phase Fluctuations 

The phase-front distortions caused by atmospheric turbulence are responsible for 
the stellar image dancing and distortion that are so familiar to astronomers. The first 
of these effects is due to changes in the angle of arrival of the received energy, a tilt 
of the phase front, and the second is due to phase-front distortions that are more com- 
plicated than a simple linear tilt. These random spatial variations of the phase are 
often the limiting factor in system performance. This is true for imaging systems 
and also for some field measurement communication systems^’ (see sec- 

tion 4.3.3). 

The purpose now is to examine the characterization of this phase-front distortion. 

In particular, a characterization that is an adequate description for the design and 
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analysis of optical communication systems is required. We have shown in sections 4.1 .1 
and 4.1.2 that for any instant of time this corresponds to characterizing the spatial ran- 
dom process <j>(r) defined by Eqs. (7) and (8), and in general we are concerned with the 
time-space process <j>(r,t). For most communication system applications it is the vari- 
ation of <(>(r, t) with r that is the important effect. Thus the major emphasis will be on 
the spatial behavior of the phase process. 

Statistics: The same Central Limit theorem argument that is used to imply that the 
set of x(r)'s associated with any set of r's are jointly Gaussian random variables (see 
section 4.1.2) also implies that the corresponding set of <j>(r)'s are jointly Gaussian ran- 
dom variables. Again, this result is often extrapolated to claim that <j>(r, t), the phase 
fluctuation process, is a Gaussian random process of time and space. Thus, the pro- 
cess is characterized by its mean and covariance functions. 

The theoretical form of the phase statistics associated with the field collected 
oyer a large aperture is unknown. It is perhaps a reasonable conjecture that the 
resulting temporal phase process is a Gaussian process. In many communication 
applications the quantity that appears is <ji(rj,t) -<j>(r 2 ,t), the phase difference between 
two points. If 4>(r,t) is a Gaussian process, then this difference is also a Gaussian 
process. 

The measurement of the absolute phase variation, a direct measurement of 
<j>(r, t) relative to a known reference phase, requires the use of a stable transmitting 
laser and a field detector, for example, a heterodyne detector with a very stable local 
oscillator. Such measurements can be used to study the statistics, mean, vari- 
ance, and spectrum of <j)(r,t) at a single point (actually over an aperture that is small 
relative to a phase coherence area) or these same quantities for larger apertures. 
Only a few measurements of this kind are available and they have been somewhat 
limited by the available laser stabilities. Because of the difficulties of holding two 
oscillators stable relative to one another and even the stabilities that can be easily 

achieved with a single laser, the same laser has been used for both the source and local 
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oscillator and relatively short two-way paths. The results thus far are mostly 

qualitative. No measurements of the probability density of 4>(r , t) for either small or 
large apertures are now available; again, this is due in part to source stability limita- 
tions. Thus there has been no direct experimental verification of the hypothesis that 
the phase process is Gaussian, and there are no measurements to show what happens 
to the phase statistics as the receiving aperture is. increased. 

Measurements of the statistics of the phase difference between two points in the 

102 

receiving aperture have been made. These measurements do not present the laser 
stability problems discussed above for the case of absolute phase measurements. These 
measures were made over a 4. 5-km path with an He-Ne laser. They show that the phase 
differences have a first-order Gaussian density. 

Mean Function and Variance: The mean function is a measure of the changes in the 
average optical path length. In most communication situations the mean of <j>(r, t) is not 
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a function of r for points in the receiving plane over the range of r of interest. 
Also the time variation of the mean function is normally slow relative to sig- 
naling times and to some extent systematic. Thus <j>(r,t) is often un im portant 
in the design or analysis of communication systems. 

An integral expression is available for the mean and variance of the phase function 

74 

for the case of Gaussian beam propagation. This, of course, also provides a solu- 
tion for the spherical wave case in the limit of the transmitting aperture shrinking 
to zero. The evaluation of this integral presents the same difficulties as the evaluation 
of x". Also, the turbulence description of Eq. (10) must be modified before it can be 
used in these calculations; that is, a suitable description for large and small values of 
p must be specified. One expects the phase variance to vary inversely as the square 
of the wavelength and linearly with the quantity /q C n (r) dX, where the integral is alon-g 
the path. 

Experimental measurements of the phase mean and variance present the same dif- 
ficulties as the measurement of the absolute phase statistics. At the present time, 
few measurements of the phase mean are available. The absolute phase data that are 
available often show a phase drift corresponding to either path- length changes or oscil- 
lator instabilities. Even when this drift is removed from the data, short-term (few sec- 
onds) measurements indicate a variance of many radians. 

Phase Structure Function and Coherence Area: The spatial variation of the phase 
can be described in terms of the phase- structure function 

D^r.r+p) = | 4>(r )-4>(r+p | 2 . (25) 


This is a particularly convenient description when the refractive index variations are 

described by Eq. (20). Often the structure function is only a function of p = |p| and can 
* 

be written as D^p). Note that D^(p) grows as the correlation between the phase values 
decreases, and when D^(p) is equal to twice the variance, the phases are uncorre- 
lated and with the Gaussian assumption independent. 

The phase- structure function has been calculated for plane-wave propagation. In 
79 

this case the result is 


D + (P) = Ak 


2 5/3 
P 



C^(r) dx, 


(26) 


When the phase process is homogeneous in the receiving plane, that is, the mean 
and variance are not functions of r, at least in the region of interest, the structure 
function and covariance function are related by 

D <t>( p) = 2 [ C cp< o) - C 4,<P>]' 

where 


C^p) = [4>(r )— 4>(r )] [4>(r+p)— 4(r+p)]. 
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where the integration is along the propagation path, and A is 1 . 46 for £ « p «'/Xl and 

2.91 for \f \ L « p « L . The parameters £ and L are associated with the turbulence 
o o o 

and are of the order of millimeters to centimeters and meters, respectively. The path 

length is L, and k = 2-rr/\. A series expression for D,(p), which is good over the whole 
82 * 

p range, also exists. The results for homogeneous paths and for a particular atmo- 
2 

spheric profile of variation are available in graphical form. They are of the same 

general form as Eq. (26), except that the exponent is larger for small values of p/\T 4L/k. 

In the spherical- wave case, a series expression for D,(p) has been calculated for 
89 * 

homogeneous paths. Again the results are available in graphical form as a function 
of p/V4L/k. D^(p) has a functional form similar to the plane-wave result. 

For Gaussian beam propagation the phase- structure function results are more com- 
plex. While an integral expression that is valid for most Gaussian beam situations 

74 

exists, it has not been evaluated. Tatarski has published a result that is applicable 
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in the near field of the transmitter. The expression has been evaluated for a homo- 
geneous path. In this case if r and r + p are assumed to lie on a line passing through 
the center of the beam, the expression for D^(r, r+p) can be written 


D,(p,R) = 0.46 k 3 aC? T n 2 / 3 p 2 In 20 2 
<b N m 
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a/2 = standard deviation of Gaussian intensity distribution 
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P = IPI 


This relation is valid for n/Xl » £ , g « 1 , and 0 » 1 . Note that it predicts a p variation 
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rather than the p ' variation associated with plane- and spherical- wave propagation. 
Also the path- length and wavelength dependence differs from that for the plane- and 
spherical- wave cases. 

The correlation function between phase variations at the beam center and at a point 

a distance p from the center has also been studied under the assumption that the 

refractive-index variations can be specified in terms of a Gaussian- shaped correlation 
103 

function. This work yields explicit results, but their applicability to the atmosphere 

is open to question because of the choice of refractive- index correlation function. 

For very turbulent situations (or very long paths) it has been suggested that the 

phase- structure function approaches a constant limit independent of the path length; that 

56 

is, it no longer increases with increasing path length. The calculation was carried out 
for plane- wave propagation, and the p and k dependence was not affected. 

Perhaps the most important single parameter associated with the spatial variation 
of the phase is the coherence area (or coherence diameter). This is the region of the 
receiving plane over which the variations of <j> are within some specified limits. If the 
validity of Eq. (26) is assumed, the coherence diameter would be defined as 


D 


coh 


Ak‘ 


Acj> 

/ 0 L c n ( ^ dx 


13/5 


(28) 


where Acj> is the allowed rms phase variation over the area. Often the phase coherence 

2 

is specified in terms of the parameter (Acj> « 6. 88) 


D 


eff 


0.0588 X 2 
J 0 L C^(r) dx 



(29) 


which plays a fundamental role in describing the performance of heterodyne communi- 
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cation systems ’ (see section 4.3.2). Values of this parameter (actually an essen- 

7 6 

tially equivalent parameter) have been calculated for a particular atmospheric model. 

For visible wavelengths, the values of D™ for a space-to- earth path are a few cen- 
timeters, while they are approximately thirty times larger at 10.6 pm. 

Few direct measurements of D.(p) or the coherence diameter have been carried out. 

9 

One set of measurements over a 4. 5-km horizonal path with a collimated 2-in. diam- 
eter He-Ne source have shown that the phase is quite coherent over apertures of 0.5 in. 

102 

with only a small amount of linear phase-front tilt. One-inch apertures, however, 
showed higher order phase disturbances. Measurements of the phase difference 
between two small apertures for separations of 2-18 in. showed little dependence on the 
separation. These results would lead one to conclude that the phase variations are 
already essentially uncorrelated for two 2-in. separations. These measurements also 
indicate a standard deviation of approximately 8. 5 rad for the phase difference. Indirect 
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measurements of D „„ have been carried out by studying heterodyne receivers oper- 
611 ^ ^ 

ating in the atmosphere. These results seem to be in general agreement with those 
that would be predicted by Eq. (26). 

Phase Front Representation: The form of the spatial variation of the phase has been 
studied both theoretically and experimentally. One representation for the entire received 
process, including both the amplitude and phase variations, is in terms of a plane-wave 
decomposition, that is, an angular spectrum. (This is discussed in section 4.2.) 

Decomposition of the phase-front variations into various geometrical shapes has also 
104 

been considered. The combination of these shapes (linear tilt, spherical deformation, 
hyperbolic deformation) in essence form the first few terms of a two-dimensional 
Taylor series expansion of 4>(r). This study showed that if the phase- structure function 
had the form of Eq. (26), the linear phase-front tilt would be the most significant form 
of distortion for small apertures. Indeed, the value of can be increased by a fac- 
tor of 3. 4 by using phase-front tracking to remove the effects of this linear phase vari- 
ation (see section 4.3.2). Corrections for the higher order distortion terms (spherical 
and hyperbolic) were not nearly as effective in terms of compensating for the phase dis- 
tortion. 


The power pattern of the energy collected over a circular aperture has been con- 
sidered for both tracking and nontracking operation, under the assumption that Eq. (26) 
105 

is valid. This study determined the average power pattern and compared the beam- 
widths and gains for tracking and nontracking apertures. The tracking beamwidth is 
approximately 30% narrower for apertures such that the rms phase difference between 
the center and the edge of the aperture is a few radians. 

Several experimental studies of the phase-front distortion over various apertures 
have been carried out by using interferometry.'* Most of these experi- 

ments have been over relatively short paths, although one experiment, which involved 

the use of a Fresnel biprism to form an interference pattern between the field collected 

1 02 

over two parts of the aperture, was carried out over a 4. 5-km path. The results of 
these experiments are largely qualitative. They do tend to confirm the importance of 
the linear phase-front tilt for small apertures and the importance of other phase dis- 
tortions as the aperture gets larger. 

Phase-Time Behavior : The "frozen atmosphere" assumption (see section 4.1.1) ties 

together the temporal and spatial phase changes. Thus the spatial variations of <j>(r) 

already considered can be used to predict the spectrum of the phase. The frequency 

spreading, caused by the atmosphere, of the output of a large- aperture field- 

88 

measurement (heterodyne) receiver has been studied. This result includes effects 
attributable to the amplitude fading, as well as the phase effects. If the phase effects 
tend to predominate, then this calculation is associated with the bandwidth of the spec- 
trum of the phase process. The results show that the atmospheric frequency spreading 
decreases with receiving aperture. 

Measurements have been made of the phase spectrum over paths up to several 
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kilometers.'*' These measurements show that the phase spectrum is at most a few 

hundred cycles wide, being only slightly wider than the amplitude spectrum. The 

temporal structure function of the phase difference between two small apertures has been 
102 

measured. The results are in agreement with the frozen- atmosphere assumption 
and indicate bandwidths of hundreds of Hz. 

Angle of Arrival: The variation of the average angle of arrival of the energy incident 
on the receiving aperture is the quantity of interest here. This is obviously just the 
linear phase-front tilt of the field over the aperture, which has been discussed above. 
This quantity is sufficiently important in its own right in such applications as acquisition, 
phase-front tracking receivers, and adaptive transmitters that it will now be considered 
in more detail. 

Under the assumption of plane-wave propagation, the variations of the average angle 

of arrival of the energy over a receiving aperture of diameter D has the mean- square 

t 87 
value 


(Ac ) 2 


= 1.027 AD' 



C 2 (r) dx. 


(30) 


where the phase- structure function is assumed to have the form of Eq. (26). Here A a 
is actually the variation of the propagation direction, from the normal to the aperture, 
along one of two orthogonal axes. Based on the phase statistics one would expect A a to 
be a Gaussian random process with zero mean and a spectral width, similar to that for 
the total phase spectrum. The total variation in the angle of arrival can with good accu- 
racy be described as composed of two orthogonal components like A a. It is thus a 
Rayleigh random process with twice the variance of Aa. The phase- structure function 
for Gaussian beam propagation, which is given by Eq. (27), can be usedto study the 
angle of arrival in the beam case. Tatar ski computes the effective angle of arrival 
over the entire received beam, but the limitation to path lengths in the near field of 

the transmitter effectively restricts the receiver size to be similar to the transmitter 
73 

size. The result is 


(Aa) 


eff 


0.1 c^i 2 / 3 


In 26 


(31) 


where the parameters and conditions for validity are specified below Eq. (27). This 
result is different from Eq. (30) both in its path length and its wavelength dependence. 
In this case the variance goes as the log of the path length. 

Measurements of the average angle of arrival have been made at the Marshall Space 

108 109 ° 

Flight Center. ’ These measurements were made with a 6328 A He-Ne laser over 

a range of 3200 m. At the transmitter, the beam diameter was 3 mm and beam diver- 
gence was 2 mrad. The angle of arrival was measured at the focal plane of a variable- 
aperture telescope with an image dissector angle error sensor. These measurements 
show that each component of the angle of arrival is essentially a Gaussian process. The 
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rms deviation in angle of arrival was measured as a function of receiver aperture diam- 
eter. The results were: 6.3 prad rms using a 5-in. receiver; 4.6 prad rms using a 
12-in. receiver. The frequency response of the angle sensor was 150 Hz. When a spec- 
tral analysis was performed on the received fluctuations it was found that the spectral 
power density fell to 10% of the DC spectral power density at 10 Hz and to 1% at 100 Hz. 

Measurements were also made in a double pass of the laser beam to a corner reflec- 
tor and back, over 6,400 m. There was no significant difference in the power spectrum, 
and the rms angle deviations appeared to be approximately doubled. Note that these 
results are in approximate agreement with Eq. (30). 

Stellar measurements can be used to study the variations in angle of arrival for a 
48 87 

Space-to- Earth path. ’ These results indicate that the rms variation varies lin- 
early with the secant of the zenith angle and has a magnitude of 1 arc second for average 
seeing conditions. For bad seeing it may be as large as 10 arc seconds. 

4.1.4 Beam Time Behavior 

The mean and variance of the intensity, for collimated beam propagation, have been 
considered at various positions in the receiving plane. The mean calculation corre- 
sponds to a calculation of the ensemble average antenna gain, which is the time aver- 
age antenna gain if the process is assumed to be ergodic. The variance calculated at 
any point includes the effects of both beam wander, changes in the position of the center 
of the received beam, and beam profile changes, as well as the scintillation effects. In 
some situations, particularly those associated with the acquisition phase of the com- 
munication problem and the use of adaptive transmitters, it is important to consider 
separately the time behavior of the variations in beam center and profile. Thus impor- 
tant questions are the following. 

1 . What are the statistics associated with the beamwidth when the beamwidth is 
measured relative to the instantaneous beam center ? 

2. What are the statistics associated with the position of the center of the beam? 

3. What is the time scale of the changes in beam center position and beam profile? 

4. To what extent can the average antenna gain be increased and the variations in 
antenna gain decreased by compensating the transmitted wavefront for the distortion it 
will suffer (for example, by using a "self- steering" phased- array transmitter tracking 
a beacon) ? 

Few theoretical or experimental results that are applicable to these questions are 
available. This is due in part to the difficulty of separating the various intensity effects 
in both analysis and measurement. The variation of the beam center is directly 
related to the average direction of propagation over the beam profile. The angle-of- 
arrival results from section 4.1.3 can be used to study this motion at a satellite by con- 
sidering the angle of arrival of the beam in a plane at the "top" of the atmosphere. Some 
feeling for the importance of the beam motion in causing antenna gain variation can 
be obtained under the assumption of instantaneous atmospheric reciprocity (see 
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Appendices G and H). Then the gain of a uniformly illuminated transmitter could be 
evaluated in terms of the intensity variations of an equivalent heterodyne receiver. 
One would then draw the following conclusions. 

1 . The gain is a log-normal process (if the heterodyne output is shown to be log- 
normal (see section 4.3.2)). 

2 . The transmitter aperture size should be limited to the diameter over which the 
received field is phase-coherent; that is, essentially uniform phase, if the gain fluctua- 
tions are to be small. 

3. At 10.6 pm, and possibly at 0.5 pm, the transmitting aperture size can be 
increased by roughly a factor of three, without increasing the gain fluctuations, by 
tracking the instantaneous phase-front tilt of a beacon and locking the transmitter 
to this signal. (Note that the applicability of phase-front tracking of a beacon to trans- 
mitter performance depends on the point-ahead angle required and on whether the atmo- 
sphere is reciprocal for two paths with this angular separation.) 

Item 2 can be interpreted to mean that beam wander is unimportant in this case and 
item 3 would be a simple case of transmitter adaptation to increase the antenna gain 
(make the beam width smaller) and still keep the effects of beam wander on gain varia- 
tion small. These results, however, are far from a complete answer, in that they are 
only applicable to the uniform transmitter illumination case, they need experimental 
verification, and the interpretation is somewhat subjective. 

4. 1. 5 Areas Requiring Further Investigation 

We shall list those areas of the channel characterization problem that seem to need 
further theoretical or experimental work and in some cases suggest experiments. The 
list should not be considered as complete, but rather as indicating ideas and areas that 
are important for optical communication work, but have not been sufficiently studied at 
the present time. Some ideas may have been neglected in the press of time. Also, future 
communication system work may point out other channel characterization ideas that need 
to be studied and may indicate that some of those now thought to be well understood 
require further work. This may be particularly true when trying to analyze some of the 
communication systems discussed in section 4. 3. Most of the ideas mentioned here 
have been considered in our channel characterization discussions. The emphasis here 
is, in general, on the Earth-Space channel. 

Theory: Further theoretical work is needed in the following areas. 

1 . Turbulence Process 

Examine further the applicability of the Rytov solution to the turbulent optical propa- 
gation problem. Specifically, effort should be concentrated on determining whether this 
approximation is adequate for providing an adequate statistical model for the communi- 
cation situations of interest. If necessary, alternative techniques applicable to beam 
propagation should be developed. 
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Examine the variation of the turbulence structure constant with atmospheric condi- 
tions and altitude. In particular, try to determine conditions that will lead to very 
turbulent layers at particular altitudes. 

Examine the communication implications, if any, of a more complete description of 
the actual three-dimensional structure of the turbulence (temperature fluctuations). The 
temperature fluctuations are non-Gaussian and skew. 

2. Mean Intensity and Average Beam Divergence 

Examine the average beam intensity pattern associated with uniformly illuminated 
apertures. 

Determine the average beam divergence and mean intensity for an Earth- to- Space 

path for both a Gaussian beam and a uniformly illuminated aperture. Consider the 

implications of the results in terms of optimal energy usage. This work will require 

2 

a realistic theoretical model of the variation of with altitude. 

Examine the applicability of stellar measurements in determining average beam 
divergence on a ground-to- space path. 

3 . Intensity 

Compare the intensity variance throughout the beam profile for both Gaussian 
beams and uniformly illuminated apertures and examine the communication implications. 
Look also for any theoretical evidence of a breakdown of the log-normal model off beam 
center . 

Examine the variation with transmitter aperture size and illumination of the inten- 
sity covariance function. This is particularly needed for uniform aperture illumination. 

Study the feasibility of using intensity covariance function measurements, or some 

2 

other indirect measurement, to determine the C „ altitude profile. 

Try to obtain a better understanding of the implications of the intensity variance 
saturation effect with particular emphasis on determining when it will be important in 
communication applications. 

Examine the spatial continuity of the temporal intensity variations within the beam. 
This may have implications for the control systems associated with tracking receivers. 

4. Phase 

Examine the dependence of the angle-of- arrival characteristics — variance, statis- 
tics, and receiving aperture dependence - on the transmitting aperture size and illumi- 
nation. 

Determine the dependence of the phase- structure function, particularly the phase- 
coherence diameter, on the transmitting aperture size and illumination. 

5. Beam Time Behavior 

Examine the magnitude, statistics, and time scale of any rapid variations in the 
beam antenna pattern, including beam wander. 
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Examine the relationship between temporal variations in the beam antenna pattern 
and temporal variations in the intensity at a fixed point in the receiving plane. In par- 
ticular, consider the possibility of adaptive techniques at the transmitter to reduce tem- 
poral variations in the antenna pattern and potentially permit the use of increased antenna 
gain while keeping the intensity variations small. 

Experiments: Many experiments and measurements are needed to verify and assist 
theoretical work. In many cases an Earth-Space path is desirable, at least ultimately, 
for the experiments. Ground-based experiments will be very useful, however, in 
establishing experimental techniques and better defining experimental goals. The exper- 
imental work should include the following areas. 

1 . Turbulence 

2 

Measure the height distribution of C^., the refractive index structure constant. The 
measurements should include direct temperature measurements, probably with balloons, 
and, if feasible, the development of some remote sensing technique . 

2. Mean Intensity and Average Beam Divergence 

Make measurements of the beam spread and the variations in beam spread. Quan- 
tities of interest are the beam divergence, the statistics, magnitudes, and time scale 
of any variations in beam spread, and the relation between beam spread and the mean 
intensity at the center of the beam. 

Measure the beam intensity profile and its dependence on aperture illumination. 

Make daylight measurements of stellar image sizes if it can be established that the 
measurements can be used to predict the average beam divergence. 

Measure the beam divergence on a ground -to -space path, using the Argon detector 
on GEOS- II. (This experiment is discussed in section 4.1.2.) Simultaneous stellar 
image measurements should be made to determine their applicability to the beam- 
divergence problem. 

3 . Intensity 

Measure the distribution and variance of the intensity variations throughout the beam 
profile. 

Measure the intensity statistics and variance for large apertures as a function of 
receiving aperture size and transmitter aperture size and illumination. 

Obtain the statistics associated with stellar scintillation, if the required results are 
not available in published works. 

Measure the intensity covariance function for stellar sources in the visible and infra- 
red spectral regions. 

Measure the dependence of the intensity covariance function on the transmitting aper- 
ture size and illumination. 

Measure the intensity time- space covariance function to determine the applicability 
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of the frozen- atmosphere hypothesis, with particular emphasis upon determining the con- 
ditions under which it breaks down. 

Measure the amplitude frequency coherence bandwidth. Such a measurement could 
be made directly with two lasers operating at different wavelengths if one can be tuned 
over an appreciable frequency range. Pulse measurements could also be used in a more 
indirect manner . 

4. Phase 

Measure the statistical distribution and variance of the absolute phase fluctuations 
caused by the turbulence. 

Measure the statistics of the phase- difference process associated with two small spa- 
tially separated apertures to check the Gaussian hypothesis. 

Measure the phase- structure function to determine its spatial separation dependence 
as a function of transmitter aperture size and illumination. 

Make measurements to determine whether the phase- structure function saturates 
with respect to its path-length dependence; that is. Does the linear path length depen- 
dence persist even for very turbulent and/or long paths? If there is saturation, mea- 
sure the phase -structure function. 

Measure the phase-coherence area, suitably defined, for ground- space paths. 

Measure the spectrum of the phase difference between two small apertures as a 
function of their spatial separation. Also measure the time- space structure function to 
check the frozen- atmosphere assumption relative to the phase variations. 

Measure the spectrum associated with the absolute phase variations caused by the 
turbulence. 

Measure the frequency coherence bandwidth of the phase process. 

5. Beam Time Behavior 

Measure the temporal behavior of the beam antenna pattern. Quantities of interest 
are the statistics, magnitudes, and time scale of variations in the beam position (beam 
wander, and the beam divergence or profile). The relationship between these quantities 
and the intensity variations at the. nominal beam center is of interest. 

Measure the beam wander to determine if the temporal variations are spatially con- 
tinuous . 

6. Phase and Intensity 

Measure the correlation between phase and amplitude variations. In particular, make 
measurements to determine under what conditions, if any [x(r j )+x(r,,)] and [4>(r^ J-^r^)] 
are independent quantities. 

7. Adaptive Systems 


Make measurements to test the concept of atmospheric reciprocity, including its 
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angular dependence. The measurements should be made for various aperture sizes. The 
possibility of using the Argon detector on GEOS- II in conjunction with a simultaneous 
stellar measurement should be considered. 

Make measurements of the increase in antenna gain that can be achieved by servoing 
the transmitter to follow the linear phase- front variations measured by a receiver 
tracking a beacon. In this case, the receiver and transmitter must share a common 
aperture. 

4.2 CHANNEL MODELS AND REPRESENTATIONS 

We shall now attempt to unify the results of section 4. 1 by considering a formal 
mathematical model of the signal portion of the field at a receiving plane, given a par- 
ticular transmitted field. Various approximate representations will then be discussed. 

A statistical channel model has previously been proposed, and it still seems to 
be useful. * *** Its validity, while predicated on the correctness of the log-normal 
statistics, is not directly dependent on the Rytov approximation. 

For simplicity, assume that the transmitted field is linearly polarized and is not 
spatially modulated. The polarization assumption can be modified by introducing a sec- 
ond orthogonal polarization component. These components can be considered separately 
because the turbulence causes negligible depolarization. It will be seen that any time- 
varying spatial modulation could be introduced in the same way that the antenna pattern 
(a fixed form of spatial modulation) is handled, that is, by considering it as a time- 
varying antenna, pattern. 

The problem is now reduced to a scalar field problem and the transmitted field is 
described simply as a function of time, as if it were an infinite plane wave. It is con- 
venient to represent this field and also the received field in terms of their complex 
29 

envelopes. That is, the transmitted and received (signal portion) fields are assumed 
to be of the form 


e t (t) = Re 


E t (t) e 


j2irvt 


(32) 


€ g (t, r) = Re 


E (t, r) e 
s 


j2TTVtl 


J’ 


(33) 


where E,(t) and E (t, r) are complex envelopes, and v is the optical carrier frequency, 
t s 

The received field is now described in terms of E t (t) by the relation 


E g (t,r) = E t (t) e 


y(t, r) 


(34) 


where 

Y(t, r) = x(t. r) + jcKt, r). (35) 

The random process y(t, r) is assumed to be a complex Gaussian random process. It 
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represents the effects of the antenna pattern, as well as those of the turbulence. As a 
complex Gaussian process, y(t, r) is completely described statistically in terms of its 
mean and covariance functions. The covariance function is described in terms of the 
covariance functions of x(t, r) and 4>(t, r) and the cross covariance function of x and 4>. 
The antenna pattern of the transmitted radiation obviously affects the means and covari- 
ances in a fundamental way and, even for the same turbulence conditions, they will 
change greatly as the actual transmitted field changes in form from a plane wave to a 
spherical wave or a beam. 

The adequacy of the suggested model rests heavily on the assumption that the coher- 
ence bandwidth of the turbulent channel exceeds the bandwidth of E t (t). (In section 4.-1 
we indicated that this is a reasonable assumption for many signal bandwidths.) Fur- 
thermore, the log-normal process assumption is being fully exploited. This has 
not been adequately checked, but at the present time it seems to be the only rea- 
sonable alternative. While this model provides the desired complete statistical 
description of the received field, it obviously places a very heavy burden on the 
mean and covariance calculations. In many cases it may be very difficult to 
obtain analytical expressions for these functions, and hence certain detailed anal- 
yses may not be possible. On the other hand, the model does contain, in a 
very compact form, the relevant features of the received field; consequently, 
it is useful in problem formulation and in understanding the over-all character 
of the problem. 

To complete the description of the total received field, it is necessary to describe 

the additive background radiation arising from natural and artificial sources. When the 

receiver field of view does not include the sun, most of the natural radiation in the 

wavelength region from the visible up to 3 pm is due to reflected or scattered solar 

radiation. The noise magnitude depends on the solar altitude and receiver attitude. In 

the visible range on a clear day, the background noise for a receiver on the Earth with 

—1 8 2 112 

a field of view centered at the zenith is approximately 10 W/cm - ster - Hz. In 
many situations the important background radiation above 3 pm is the thermal radiation 
of the terrestrial environment whose spectral shape is reasonably modeled by «280 °K 
black-body radiation. 112, 11 3 Around 10 pm the background level is 10 W/cm 2 - 
ster - Hz. There is some variation in the background level above 3 pm with atmospheric 
conditions (clear or overcast), but little dependence on the time of day. At the 
visible wavelengths, the amount of background is clearly time- dependent and greatly 
reduced at night. In some applications direct sunlight, reflections from planets 
or the moon, or stellar radiation may be important. Information on these sources 
is available. 111 

Normally, it will be appropriate to model the background radiation field as an 
unpolarized field whose components in the plane of the receiver are independent. 
Moreover, the complex envelope of each polarization component can be considered to 
be a zero-mean complex Gaussian random process with independent components that 
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are stationary in time and space. Furthermore, it is reasonable to assume that the 
time- space covariance function of the real and imaginary parts are identical, and that 
the spatial covariance function is only a few wavelengths wide. Finally, the time spectra 
of the real and imaginary parts can reasonably be considered to be flat over most band- 
widths of interest. 

At this point, the channel characteristics can be summarized. The model corre- 
sponds to a fading channel with little time or frequency dispersion. It differs from the 
more familiar Rayleigh fading channel, in that the amplitude variations have a log- 
normal rather than a Rayleigh distribution, and the phase variations have a Gaussian 
rather than a uniform distribution. The variance associated with the phase process may 
be large. 

The model discussed above presents severe and largely unsolved analysis problems. 
Using it to obtain optimum systems is very difficult. One reason for this difficulty is 
the unsolved problem of obtaining a suitable mathematical representation for the sum 
of log-normal and Gaussian processes. Also, the very compactness of the model may 
tend to obscure certain features of the received process that good suboptimal com- 
munication systems might exploit. Thus, while the model is very useful, there is a 
great need for reasonable approximations and alternative representations that retain 
the relevant statistical details but increase the feasibility of both optimum analysis 
and the design of good suboptimum systems. This need for a sufficient understanding 
to make reasonable approximations is the motivation for the detailed consideration of 
the various turbulence effects that we have carried out in section 4.1. 

The problem of formulating meaningful and useful approximations to this model or 
obtaining useful alternative representations is certainly not a closed issue. Several pos- 
sible approaches either have been used or look promising. One approximation that can 
be used to reduce the problem to manageable form is to partition the receiving aperture 
into coherence areas and to consider the field amplitude and/or phase as constant (coher- 
ent) within each area and statistically independent from area to area. At some cost in 
complexity, the areas may differ for the amplitude and the phase. This reduces the prob- 
lem to the consideration of a set of independent phase normal processes, that is, pro- 
cesses defined as the exponentiation of complex Gaussian processes. If the time variation 
is considered in the same way, defining the process to be constant over a coherence time 
and perhaps independent from coherence time to coherence time, then the problem is 
reduced to the consideration of a set of phase normal random variables. Such a 
description obviously involves a lot of approximation, but it does seem to retain the 
essence of the statistical description and it has proved useful in analyzing communi- 
cation systems.** * 


In some situations, it is appropriate to consider the background field as a constant 
intensity. The applicability of this approach depends on the relative bandwidths of the 
optical predetection filter and the final electronic receiver processing. 
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A process model that looks potentially useful and should be given further theoretical 
and experimental study is a representation in terms of the angular spectrum of the 
received process. This might take the form of trying to represent the field over the 
receiving aperture in terms of a set of plane waves characterized by random direction 
vectors and Gaussian- shaped intensity distributions with independent amplitudes. A 
representation of this kind would be useful in the analysis and synthesis of systems. 

An interesting aspect of the representation problem is the relation between aper- 
ture plane representations and representations in the focal plane of a lens. This 
obviously provides a tie between the ideas discussed above. The number of degrees 
of freedom should be preserved by the linear lens transformation. It appears that, 
at least for a large lens, the degrees of freedom, which are represented as coher- 
ence areas in the aperture plane, are represented by Airy disk-sized areas in the focal 
plane. The number of Airy disk- sized areas in the focal plane with significant energy 
would then be equal to the number of coherence areas in the aperture plane. This rela- 
tionship should be examined theoretically and experimentally. 

The present discussion suggests the need for additional theoretical and experimental 
work with respect to the channel representation problem. Theoretical work of interest 
should include the following investigations. 

1 . Study of the ergodicity question for the received field process and determination 
of the communication implications of non-ergodicity. 

2. Study of the relationship between aperture and focal plane representations of the 
received field, particularly with reference to the form in which the degrees of freedom 
appear. 

Study of the representation of the received process in terms of an angular 
spectrum of plane waves, perhaps with nonuniform intensities. This study should 
include an examination of the implications and validity of such a representation both as 
an exact model and as a reasonable approximation. Particular emphasis should be 
placed on the statistics associated with the representation. 

Experimental work should include measurements to determine the applicability of 
an angular spectrum representation for the process and the statistics associated with 
such a representation. Such measurements might be made by performing simultaneous 
measurements of the amplitude and phase processes associated with an array of detec- 
tors appropriately located in the focal plane of a receiving lens. 

4.3 COMMUNICATION SYSTEMS 

Here we consider communication systems for use with the turbulent atmospheric 
channel. The emphasis is almost completely on digital communication systems. While 
a fundamental quantum detection viewpoint would be pursued, ideally, it seems more 
fruitful at this time to concentrate on systems that use direct (energy) detection or some 
form of heterodyne (field) detection. The use of optical preamplifiers before the final 
detection operation will also be considered. Several novel communication techniques 
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that can be used with any detection scheme, and with analog, as well as digital, systems, 
are also discussed. All of the systems and techniques considered here have promise in 
some applications and merit further investigation. 

4.3.1 Direct-Detection Techniques 

We are concerned here with atmospheric systems that use direct (incoherent or 
energy) detection. This type of detection process was discussed in section 3.2.3 in 
the context of the free-space channel. Since many of the results discussed there are 
directly applicable to the turbulent- channel case, the major emphasis here will be on 
those additional features that arise because of the turbulent channel. 

A fairly general model for the direct-detection problem will be introduced. It will 
provide a framework within which to discuss various important features of the problem. 
Various promising receiver structures, the signal design problem, and useful experi- 
ments will then be considered. 

A Model for Direct-Detection Systems: The model shown in Fig. 4 is useful in dis- 
cussions of direct detection techniques for the turbulent channel. For simplicity, a 
scalar field, linearly polarized, has been assumed. As indicated in Fig. 4, D signals, 
represented by the components of the vector rjt), are available to be processed. Each 



Fig. 4. General model for structured atmospheric channel receivers. 


component is the sum of a conditional Poisson process, representing the photodetector 
output as described in section 3.1, and a signal- independent Gaussian noise process 
representing the noise associated with the detector and its associated circuitry, for 

sjc 

example, the electronic preamplifier. The inclusion of several photodetectors allows 
the possibility of using separate detectors for various spatial regions or modes of the 
received field. The input to the system, the portion in dashed lines, is the transmitted 
signal. The field e g (t, r) represents the signal portion of the field at the receiving plane 


sjc 

In general, for nonideal devices the photodetector output is a conditional filtered 
Poisson process and there is also an additive signal- independent shot process (dark 
current). (See section 3. 1 . 1 .) 
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after it has propagated through the turbulent channel. The relationship between the 
transmitted field and e (t, r) was discussed in section 4,2. The field e, (t, r) represents 
the background noise. 

The model allows for linear processing on the field before the energy detection oper- 
ation; the additive noise source e n (t, r) accounts for the possibility that this processing 
introduces noise. The field processing may, for example, be time filtering with an opti- 
cal bandpass filter, spatial filtering with some optical system, amplification with an 
optical preamplifier, or some combination of these operations. This processing might be 
used to decompose the field into various spatial modes. The processing noise e n (t,r) is 
especially needed to model adequately the optical preamplifier case (see section 3. 1 .2). 

As discussed in section 3.1, the process r(t) depends on the intensity (instantaneous 
power) of the field at the detectors, e^t, r), and on the Gaussian signal- independent noise 
associated with electronic operations on r(t). The direct detection technique is most 
attractive when the additional noise caused by the electronic operations is negligible 
compared with the signal-dependent noise. This is often the case when the detector has 
internal gain, for example, a photomultiplier, and thus with the present technology is 
most attractive in the visible and near infrared wavelength region. On the other hand, 
the situation may also occur if the linear field operation includes sufficient amplifica- 
tion — a technique most applicable in the longer wavelength region. 

Assuming, for the present, that the detectors are noiseless (the additive Gaussian 

noise can be ignored), consider the noise caused by the detection operation on e^ft, r). 

Several situations are of interest. One is when the background field e^ft, r) is sufficiently 

large relative to the signal field € (t,r) that it dominates e (t,r), with e (t, r) = 0. In this 

s o n 

case the "shot noise" in r(t) will depend essentially on the filtered version of e^tt, r) and 
will be signal- independent. If the rate is high enough relative to subsequent processing 
time constants, it may be appropriate to consider this shot noise as a Gaussian process. 
Also, if the input optical filter bandwidth is large relative to the bandwidth associated with 
the subsequent processing, it may be appropriate to consider the background intensity as 
a constant. Another case is when an optical preamplifier is used. Normally, for this 
case, e n (t,r) is the dominant input noise to the detectors. Finally, the signal field, espe- 
cially after the linear processing, which will normally limit the field of view and contain 
an optical bandpass filter, may be larger than any noise field. In either of the last situa- 
tions, the "shot- noise" processes r(t)must normally be considered in detail in any anal- 
ysis. 

As noted in section 3.2, the forms of the optimal processors are available for known 
signals in the presence of a noise whose bandwidth greatly exceeds the post detection 
bandwidth, that is, where the background noise can be treated as a fixed inten- 
sity.* 14-116, 36 g ome p er f orm ance results for the case in which the background inten- 
sity is large relative to the signal intensity are also available. 1 1 ^ f n the problem of 
interest here, however, the signal fields are random, and it seems that only a few spe- 
cial situations are understood. Before discussing specific direct-detection systems, we 
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first enumerate the questions that would ideally be asked. For digital signaling with 
a probability- of- error criterion the relevant issues, in terms of the model are these. 

1 . The optimum receiver structure, that is, processing onjr(t), and the resulting 
performance . 

2. Optimization of the linear field processing within some set of allowed operations 
to minimize the probability of error. For example, only passive devices may be allowed 
for this processing. 

3. Optimization for a given receiving aperture of the number of photodetectors and 
of any spatial field processing performed by the detectors. 

4. Optimization of the signal set with respect to the minimum probability- of- error 
criterion. 

It is necessary, in general, to consider the joint optimization problem. 

Analog communication and parameter estimation problems are also of interest. 
Among the important problems that should be considered are the FM communication 
problem and the synchronization problem (time estimation) for digital communication 
systems. FM systems look attractive because of the relatively small effect that the 
turbulence has on the frequency. Some of the more promising digital systems, such 
as PPM, require accurate timing, and thus good synchronization, even though one ter- 
minal may be moving relative to the other. Of course, optimizations of the type speci- 
fied for digital signaling are the desired result. 

Several receiver structures, suggested by the available analysis results, 'by the 
channel properties discussed in sections 4. 1 and 4. 2, and by previous experience with 
channels that have similar properties, will now be considered. In general, they will 
be considered in order of increasing complexity and presumably performance. As indi- 
cated above, few theoretical results are available, but the structures considered are all 
sufficiently promising in terms of simplicity or performance to warrant further study. 
Specifically, their performance and the trade-off between complexity and performance 
should be established. Of course, other promising implementations may also exist. 
The systems considered here reflect the interests and experience of those at the 
Workshop. 

The property of the received field that most directly affects the design and perfor- 
mance of the receivers considered here is the turbulence- induced temporal and spatial 
variations in the magnitude of the complex envelope rather than its phase variations. 
Thus the receivers that we now consider will try to combat the effects of these envelope 
magnitude variations. 

Single Detector (Photon Bucket): The simplest receiver consists of a single detec- 
tor followed by electrical processing of the detector output. Such a receiver has been 
considered in section 3. 2. 3. For the case of a known signal and background noise 
which is modeled as a fixed intensity, the detector output processing is known to be 
essentially a weighted counting operation and the error performance has been evaluated 
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for certain signal sets. When the signal is a random process, as it is for the 

atmospheric channel, much less is known about the optimum receiver and its perfor- 
mance. The following theoretical results provide some insight into the characteristics 
of a good system of this type. 

The maximum- likelihood receiver for the problem of detecting a spatial Gaussian 
process that does not vary over the transmitted pulse length in the presence of an addi- 
tive Gaussian noise field that is spatially and temporally white has been considered in 

119 

the absence of any real detector effects. The receiver computes a quadratic 
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Fig. 5. Single detector (photon bucket) receiver for atmospheric channel. 


functional of the received data. By extrapolating from this analysis, we conjecture that 
a receiver of the form shown in Fig. 5 will provide good performance for the detection 
of Gaussian fields, even when the real detector statistics are considered. In Fig. 5, 
the optical filter approximates, for the pulse signal, the temporal matched-filter oper- 
ation that should be performed on the field. As it is now possible to build optical fil- 

O O 

ters with a 1 /20 A bandwidth, at least at 6328 A, the optical filter may be a satisfactory 
approximation to a filter matched to the pulse. Such filters have a field of view limited 
to approximately a degree, but this should not present a problem for the turbulent chan- 
nel. The spatial filter has a transmission function that is roughly analogous to that of 
a field stop, which limits the receiver field of view. In the Gaussian case, the spatial 
transmission of this filter is 

_ / TR(p) \ 1/2 

t(p>= — 

\N 0 + TRtp)/ 

when the covariance function of the signal field, R(p), is narrow relative to the aperture 
diameter and T is the pulse length. Here the tilde denotes a Fourier transform. The 
exact time processing on the detector output, of course, depends on the signals used, 
but it seems reasonable to expect that it will be some form of weighted counting. 

Some theoretical results, again with real device characteristics ignored, seem to 
indicate that if the signal is spatially phase-normal rather than Gaussian, that is, it has 
log-normal amplitude and Gaussian phase processes, more complicated processing 
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involving additional quadratic functionals is required. These results, which assume 
that the aperture is large relative to the correlation distance, would call for two quadra- 
tic functionals when only phase variations are important, and three when both phase and 
amplitude variations must be considered. It is possible that the form of these results 
is due in part to the analytical approximations that were required to make some prog- 
ress with this very difficult problem. Next, an ad hoc receiving structure based on the 
Gaussian detection results discussed above will be considered. 

On an intuitive basis, if only a single detector is used, it would seem that the receiver 
shown in Fig. 5 should also be a good choice in the turbulent channel case, that is, when 
the signal field is phase-normal. While no theoretical results are available to specify 
the filter for this case, it is reasonable to expect that the filter should be matched to the 
turbulence so that only the detector area where appreciable signal energy is expected 
is exposed. In terms of the channel models of section 4.2, this filter could be designed 
by considering the plane wave or angular spectrum representation of the field over the 
aperture or, for sufficiently large apertures, by using the focal plane degrees of free- 
dom model previously suggested. Note that these representations provide an indica- 
tion of the physical size of the detector. Also note that "poor" receiving optics, where 
the quality is judged relative to the field perturbations caused by the turbulence, will 
cause an unnecessary increase in the size of this field stop and hence admit more back- 
ground noise. The time processing, of course, is in general different in detail from that 
for the Gaussian field, although it is still probably some form of weighted counting. 

No general performance results are available for direct detection with the turbulent 

channel. Some results are available for binary on-off signaling when the aperture is 

limited to a single amplitude coherence area, and for orthogonal signaling with opti- 

121 122 

mum quantum detection. ’ Also, the variation of error probability with log- 

amplitude variance has been considered for a PCM/PL system operating with a con- 

118 118 
stant background intensity. It has been assumed by Peters and Arguello that the 

intensity fluctuations for large apertures are log-normally distributed, which still 

requires further confirmation. The same assumptions could be used to get other special 

case results. Some experimental direct- detection measurements for a turbulent channel 

have been m'ade^^ and others are in progress. 

Array of Detectors : The potential of using more than one detector can be visualized 
by using the aperture plane representation of the received field in terms of a set of 
coherence areas. For the present case, the field amplitude is the important quantity, 
and it is the amplitude coherence diameter that is important. The use of a detector array 
permits a more optimal use of the spatial diversity available in the signal field than 
may be possible with the system of Fig. 5. One possible form for the array receiver is 
shown in Fig. 6. The two lenses and spatial filter form a field of view filter to limit 
the background noise. The spatial filter field stop, as in the case of the single detector, 
would be adjusted to the angular spread of the signal field. The output of this lens 
system is an image of the aperture with reduced background noise. As in the 
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single-detector case, the optical bandpass filter acts as a matched filter on the field 
and further reduces the noise. 

The array of detectors is so chosen that each individual detector diameter is equal 
to the amplitude correlation distance. [Note that additional optics could be used to 
reduce the required detector size (minify the image).] The output of each detector is 
individually processed before the outputs are combined to make a decision. For digital 
signaling with an alphabet of size M, the time processing associated with each detector 
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Fig. 6. Direct-detection array receiver for atmospheric channel. 

channel will have M outputs, where each is associated with one of the possible trans- 
mitted signals. The outputs corresponding to each particular signal are added together 
to form the decision variables. The exact processing on each diversity channel is not 
known. In general, it will depend on both the signal set and the variance of the log- 
normal amplitude fluctuations. Intuitively, the reason that individual detector processing 
permits more efficient use of the spatial diversity is that it can perform an operation 
on the data to suppress the effect of those detectors, and hence aperture areas, that 
are dominated by noise; that is, in forming the decision variables the outputs of the 
high intensity signal areas are weighted more heavily than the outputs from the areas 
with severe fading. Thus the processing attempts to estimate the fade associated with 
each diversity path using the data from a single baud time and uses these estimates in 
forming the decision variables. 

The performance of this type of receiver has not yet been studied. A particularly 
important result that should be obtained is the gain in performance, over that of a single 
detector, which can be attained by the use of a detector array with individual processing. 
It is not clear how substantial the gain will be . 

Channel Measurement Receiver : A channel measurement receiver is an extension 
of the detector- array receiver. In the array receiver the decision procedure weights 
the contributions of various portions of the collecting aperture according to the gain 
associated with that portion. A more efficient operation would, of course, be possible 
if the gains were actually known. This suggests that the receiver should try to esti- 
mate the log-normal channel amplitude space-time process and then use this esti- 
mate in carrying out the signal processing. This approach has proved to be a useful 
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technique for other channels. It is potentially very attractive for the turbulent channel, 
because of the very slow (1 msec) fading rates that allow long time averages to be used 
in the channel measurement procedure with a resultant decrease in measurement noise. 
A pilot tone can be used to facilitate the channel estimation, for example, an unmodu- 
lated signal at a nearby wavelength, but it is also possible ordinarily to use the infor- 
mation signal directly for this estimation. The last procedure is often more convenient, 
and for some channels more efficient in terms of total signal power. 


INPUTS FROM OTHER 
DETECTOR CHANNELS 
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Fig. 7. Direct-detection channel measurement receiver for atmospheric 
channel. (Only a single diversity element is shown and binary 
signaling is assumed.) Jj(t) and J ^ (t) are the intensities asso- 
ciated with the two messages. J Q is the background intensity. 

As an example of the ideas involved, consider the receiver shown in Fig. 7. It is 
assumed that binary signaling is being used, but the modification for the M-ary case is 
_evident. Only a single one of the spatial channels is shown in detail; there would in gen- 
eral be as many channels as in the array receiver just considered. The detector would 
be preceded by the field of view and wavelength filters shown in Fig. 6. It should be 
emphasized that even the limited structure shown here is an ad hoc configuration based 
on conjecture rather than analysis. The heart of the system is the gain estimator box 
that uses its knowledge of the possible signals and the background level to form a 

2 X j ( t ) 

noisy estimate of e , the channel intensity gain. The key to its operation is an 

averaging or filtering of the signal over a time corresponding to the channel amplitude 

coherence time. This gain estimate is then used in the matched-filter (weighted counting) 

* 

receiver that would be optimum if the gain were exactly known. 

Channel measurement receivers look promising and certainly warrant further inves- 
tigation. The most desirable form of the receiver and its performance should be deter- 
mined. Emphasis must be given to the effect of signal sets on performance; 

❖ 

These ideas are based upon an exposition which P. E. Green submitted at the Work- 
shop (see Appendix I). 
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Signal Design: The signals for a direct- detection system must be designed to permit 
discrimination between the signal states on the basis of energy alone. Any form of 
amplitude modulation is certainly applicable, as is polarization modulation and frequency 
shift keying. In the last cases, predetection processing is required to convert the infor- 
mation into an energy form. For binary signaling with polarization modulation, this 
prior processing might take the form of a Wollaston prism to spatially separate the 
two polarizations; for frequency- shift keying, a set of bandpass optical filters would 
probably be used. 

For direct detection with real devices the performance in general depends on the sig- 
nal shape, as well as on its energy. Direct- detection performance for digital signaling 
with the free- space channel is in general improved by the use of pulsed operation with 
high peak intensities. None of the characteristics of the turbulent channel would seem 
to change this situation. The signal set that minimizes the probability of error is not, 
in general, known. On the other hand, characteristics of the turbulent channel are such 
that most of the signal sets that have been considered for the free- space channel can 
still be used effectively. Pulse position modulation which is attractive for the free- 
space channel should still be very attractive for the turbulent channel. This reflects 
the fact that the time spreading of the channel is very small, and that the frequency 
coherence bandwidth is very large. The fact that the turbulence causes almost no depo- 
larization makes polarization modulation attractive. The very small frequency 
spreading caused by the atmosphere makes frequency- shift keying look attractive. This 
technique is limited by present modulation technology. 

Experiments: Several of the more important experiments for direct- detection 
receiving systems await further theoretical work. This is particularly true of experi- 
ments associated with arrays of detectors and with channel measurement. A better 
understanding of the optimum system implementation and of the design trade-offs is 
required before the ideas can be meaningfully tested on the turbulent channel. It 
should be emphasized, however, that such experiments will be necessary to assess 
the many approximations involved in the system models. Furthermore, the informa- 
tion about the intensity spatial covariance function provided by some of the experiments 
discussed in section 4.1 will be required for an analysis of these systems. 

On the other hand, some meaningful experiments could be carried out immediately 
with the single- detector receiver. These include performance measurements 
(probability- of- error measurements) for various turbulence and background conditions 
to study the optimization of the field stop (spatial filter) and the effects of aperture size. 
Also, measurements of the intensity distribution for large apertures are needed to ver- 
ify the extent to which it remains log-normal. This information is needed for the the- 
oretical performance calculations. The experimental situation should, of course, 
simulate the space- ground path and include a study of the effects of transmitter aperture 
size and type of aperture illumination. In making error performance measurements, 
care should be taken to optimize the processing of the detector output signal with 
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respect to the digital signal set that is being used. Signal sets should not be limited to 
binary. For example, M-ary PPM should be considered. 

4.3.2 Heterodyne Detection Techniques 

We now consider turbulent- channel communication systems that use heterodyne 
(field) detection. This type of detection was discussed in section 3.2.1 for the free- 
space channel and some of those results are directly applicable to the turbulent chan- 
nel. The emphasis here will be on the features that arise because of the turbulence. 

Thus we focus on techniques to overcome the effects of the phase-front distortion and 
amplitude fading. 

We shall consider heterodyne detection within a fairly general framework. This 
will permit an examination of some of the relevant features. Subsequently, we shall 
consider important and/or promising receiver structures, signal design, and experi- 
ments that would be of interest. 

Model for Heterodyne Detection Systems: The model of Fig. 4 can be applied to 
heterodyne detection by replacing the field e n (t, r) by a heterodyning field e^/t, r )- The 
actual heterodyning field is sometimes introduced ahead of the linear field processing, 
but this causes no fundamental change because this processing is known and linear. 
Although the complex envelope of the heterodyning field often has no spatial or time var- 
iations, the additional generality provided by these variations is sometimes required. 

As discussed in section 3.2.1, it is sometimes appropriate to approximate the vector 
process r(t) by a term that depends on the field e Q (t, r) (and €j/t, r)) and an additive 

— 12 5 

Gaussian noise term that depends only on e^/t, r). In this case, all of the detection 
and communication theory results for random signals in additive, signal- independent, 
Gaussian noise become applicable. The situation is still complicated, however, by the 
fact that the random signal process is log-normal (or log-normal derived). 

Another source of noise may arise when heterodyne detection is used. This is the 
noise associated with the heterodyne field e^/t, r). If this field is generated by a laser 
local oscillator, it is usually reasonable to assume that it contains no undesired varia- 
tions. If, however, this field is obtained by amplifying part of the transmitted signal, 
this noise will be very significant. 

The phase-front distortions caused by the turbulence are extremely important when 
heterodyne detection is used because the detector output corresponds to an integral of 
the field over the detection surface. If the field is not spatially coherent over the detec- 
tor surface, the varying field phase will cause large variations in the detector output 
signal strength and severely affect performance. Thus a primary concern in designing 
heterodyne receiver systems is attempting to minimize the effects of the phase-front 
distortions and still use large receiving apertures to increase the collection of signal 
power, and hence error performance or communication rate. While the phase-front 
distortions will be most influential in determining the form of heterodyne receivers, the 



ATMOSPHERIC CHANNEL 


59 


amplitude scintillations will be very important in determining their performance. Thus 
both effects will be important in the examination of the receivers that follow. 

We shall now consider several heterodyne receiving structures, beginning with the 
optimization and performance of the simplest heterodyne receiver and preceding to 
more complex and powerful systems. The objective is to survey those ideas that look 
most promising either from the point of view of ease of implementation, or perfor- 
mance, or some trade-off between these two factors. These receivers, of course, 
try to exploit the channel characteristics discussed in section 4.1. In some cases, 
the receiver structure is based on analysis, while in others it is suggested by the chan- 
nel characteristics and previous experience with similar channels. 

Single Detector: We consider, first, the receiver that uses a single photodetector. 
This detector would be preceded by optics to limit the field of view to the angular spread 
expected in the signal portion of the field. If the receiving aperture is sufficiently small, 
so that the phase is essentially constant over it, the detector output will correspond to 
a log-normal signal in independent additive white Gaussian noise. The error perfor- 
mance of such a system, for binary orthogonal signaling and baud times less than 

the channel amplitude coherence time, has been calculated numerically and also 

1 2 1 \ 2 6 > 2 

bounded. ’ These results are a function of the parameter tr^, the variance of the 

normalized log amplitude of the received field. As would be expected, if the value of 
2 

<r is nonzero, more average signal power is required to achieve the same performance 
X 2 

as in the nonfading, c r = 0, case. The situation just considered is essentially true for 

A. 

any reasonable receiver aperture on a spacecraft, because the field there normally has 

large amplitude and phase- coherence areas. 

For the heterodyne receiver on the ground, as the receiver aperture increases, signal 

variations caused by the phase variations of the received field start to become important. 

This effect has been studied both theoretically and experimentally.^’ * 23 ’ * 2? ’ * 2 ® 

These studies show that the average signal power saturates as the receiver aperture 

increases, and that for any aperture size, it will not exceed the power that would be 

63 76 

associated with an aperture of diameter (or r Q ) ’ (see section 4.1.3). The curve 
of average signal power versus diameter has a knee at r Q and is approximately 3 dB 
below the saturated value for an actual receiving diameter of r . For a homogeneous 
path, r Q is 63 ’ 76 

r = 1.2 X 10- 8 X 6 / 5 L- 3 / 5 c; 6 '' 5 , 
o p N 


where X. is in microns, and L is the path length. It has also been evaluated for an 
inhomogeneous space-to-ground path using an assumed model for the refractive index 
variations. Again, it depends on wavelength to the 6/5 power. In summary, the aver- 
age signal power for a single heterodyne detector is limited, and a larger receiving 
aperture can be effectively used at 10.6 pm than in the visible wavelength range. Studies 
have shown that performance is actually degraded when the receiving aperture diam- 
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eter approaches or exceeds r ^,121 s ig na q variance increases rapidly with 

aperture size beyond r Q ; the error performance obtained if the aperture is assumed to 
contain two independent phases is significantly degraded. To the extent that the resulting 
signal is accurately modeled as a log-normal process, the signal variance results can be 
used with the previously mentioned results for the detection of log-normal variables in 

white Gaussian noise to evaluate the performance of these larger aperture receivers. 

70 123 

While some theoretical and experimental results are available, ’ the use of the log- 
normal model for the received signal as the aperture size increases should receive fur- 
ther investigation. These considerations tend to indicate that the aperture diameter 
should be r Q or less, and that good error performance may require that it be limited to 
perhaps r / 2. 

The results just presented are for the case of a uniform local- oscillator field. The 

effective aperture area can be increased by using a local oscillator whose intensity is 

matched to the average phase properties of the atmosphere, that is, e^ft, r) now actually 

varies with r. This is equivalent to antenna beam shaping and the gain that would be 

expected is in general small (perhaps 10%). This technique for improving system per- 

129 

formance has been studied. (This technique is to be distinguished from instantaneous 
phase front matching by the heterodyning signal where the gain may be extremely 
large. We consider this in section 4. 3. 3.) 

Single Detector with Angle- of- Arrival Tracking : The fact that the predominant 
phase-front distortion for small apertures is a linear tilt has been discussed in sec- 
tion 4.1.3. Thus, one way to increase the effective aperture for a heterodyne 
receiver, over that which can be used with a single static detector, is to remove 
the effect of this instantaneous phase-front tilt. This can be accomplished by tracking 
(measuring) the instantaneous angle of arrival and matching the heterodyning field to 
it. The range of values and frequencies associated with this angle of arrival have been 
considered in section 4.1.3. Possible techniques for accomplishing the tracking and 
introducing the local oscillator have been considered. 

The performance improvement that can be attained with tracking heterodyne 

receivers has been considered theoretically.*^’ 130 ipj^egg studies consider the 
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resulting power pattern and the increase in receiver aperture or average signal 

104 105 

power that can be obtained by tracking. ■ The results show that. the effective aper- 
ture diameter is increased by a factor greater than three, and that the resulting average 
signal power is increased by more than a factor of 10. Thus, while the static receiver 
is limited to a diameter such as r Q /2 by the signal fluctuations, it is likely that a diam- 
eter of the order of 1 . 5 r can be used in the tracking receiver. This last statement 
is based on the conjecture that the signal fluctuations associated with a tracking receiver 

are of the same order of magnitude as those associated with a static receiver whose 

* 

aperture is approximately one third as large. No theoretical or experimental results 
concerning the signal fluctuations of a tracking receiver are available at the present 

* 

This conjecture was suggested by D. L. Fried 
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time. To the extent that the signal fluctuations are reasonably modeled as a 
log-normal process and that their variance can be predicted from the static receiver 
case, the performance of the tracking receiver can be obtained from results discussed 
in conjunction with the static heterodyne receiver. These conjectures would lead to the 
conclusion that the same error performance can be obtained with a tracking receiver 
as with a static receiver, with a reduction of more than 10 dB in signal power. 

Array of Detectors: It may often be desirable to use receiving apertures that are 
larger than those that can be effectively used by a single detector receiver. One technique 
to increase the aperture size that can be used with heterodyne detection is the use of an 
array of detectors with separate processing on each detector output. Such a receiver 
would be modeled by Fig. 4 when e n (t,r) is replaced by the heterodyning field, e^(t,r), and 
each of the detector outputs r.(t) results from a phase- coherent portion of the signal field. 
Following the discussion of section 4. 2, the detectors may be located in the image plane 
of preceding optics, which limits the field of view to that occupied by the signal energy, or 
they may be located in the focal plane of the optics. In the former case each detector area 
must be limited to the size of the imaged phase- coherence area (which is assumed to be 
less than or equal to the amplitude- coherence area), while in the latter case the detector 
size is determined (at least for large apertures) by the diffraction pattern of the receiving 
optics. In either case, the total number of detectors is essentially the same. 

The optimum receiver must process the vector output r(t) in such a way as to make 
a minimum probability of error decision in the digital communication case. This 
requires that the output of each detector, each component of r(t), be separately pro- 
cessed. The optimum use of spatial diversity and the requirement for separate pro- 
cessing of each detector output has already been discussed in section 4.3.1 in conjunction 
with direct- detection receivers, but the individual detector processing is more impor- 
tant, indeed crucial, in the heterodyning case. The form of the processing that is 
required has been studied under the assumptions of large phase variance, equal- 
amplitude and phase-coherence areas, statistical independence of the amplitude, and 
.phase values associated with the various coherence areas, and decision times short 
relative to the channel coherence time.**'*' The resulting receiver, which is min-max 
relative to phase, requires at each detector output a set of matched filters followed by 
envelope detectors and a nonlinear operation. The receiver is shown in Fig. 8. There 
must in general be a filter matched to each signal in the allowed signaling set (cor- 
relation operations could also be used). The processed outputs associated with the 
same signal are added together after the nonlinear operation to form the decision var- 
iables. The receiver chooses the message corresponding to the largest sum. The 

required nonlinear operation has the form of a biased squaring operation which depends 

2 

on cr , the variance of the normalized log amplitude of the received field, and the 

X 2 22 

signal- energy- to-noise spectral density ratio per diversity path. 

The performance of this diversity receiver, with the same assumptions, has 
also been studied.*** It has been shown that the probability of error, Pr [e], for 
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SAMPLE 
AT t-T 


h,(ll ■Re[s 1 “(T-t)e' 12 ’'"] 


Fig. 8. Min-Max phase receiver for heterodyne detection with an array 
of detectors. (The S.(t) are the complex envelopes of the signal 
waveforms . ) 

a signaling alphabet of M signals, which are orthogonal at the receiver, can be written 
in the form 


B 2 -KE < Pr [e] < 2 KE , 

where K is the number of information bits per transmitted waveform, and E is called 
the "system reliability. 11 The exponent E depends on p, the average received energy- 
to-noise spectral density ratio per information bit, a , the averaged received energy- 

P 2 

to-noise spectral density ratio per diversity path, and cr . the variance of the normalized 

AW X 

log amplitude of the received field. The coefficient B also depends on K, p, and a . 

A reasonable understanding of system performance can be obtained by using the expo- 
nent KE as a performance measure. It can be shown that, for fixed values of K and 
0 ^, KE increases at least linearly with increases in the number of coherence areas con- 
tained in the receiving aperture. ’Thus the error exponent increases at least linearly 

with increasing receiver area. The fact that KE depends on a makes diversity con- 

P 

trol desirable. While the spatial diversity is not completely controllable by the sys- 
tem designer, some control can be exercised by varying the transmitter power or antenna 
gain; also time diversity can be used to augment the available spatial diversity if this 
is necessary. If appropriate values of can be achieved, the minimum value of p at 
which arbitrarily reliable communication is possible is In 2; that is, the turbulent chan- 
nel, when used with orthogonal waveforms, has the same channel capacity as an 
equivalent nonfading white Gaussian noise channel. 
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Numerical techniques have been used to study, for large values of (3, the KE vs a 

2 P 

curves as a function of <r . The results show that a broad range of a values will yield 

X P 

relatively good performance, but that too much diversity can greatly degrade the per- 
formance. Furthermore, the results show that, even for the optimum value of a which 

2 P 

is a function of ( r , the error probability for a fixed alphabet size (fixed K) and a 
X. 

receiver of fixed area is increased by the turbulence. Thus, for reasonable turbulence 
levels and near optimum a values, 3-4 dB of additional signal power at the visible fre- 
quencies, and somewhat less in the 10.6-pm portion of the spectrum, are required to 
achieve the error probability that would result on a nonfading channel. Similar numer- 
ical techniques could be used to study the behavior of KE for other ranges of |3. 

The receiver structure and performance results obtained above should also be 
applicable to the nonheterodyning case in which the background noise is very large; 
that is, it is so large that it dominates the signal and the filtered shot process at the 
detector output looks Gaussian. Similar results are probably also applicable when 
an imaging optical preamplifier is used in the receiver front end (See sec. 4.3.3, "Ave- 
"Avenues for System Development," p. 65). 

Detector Array with Angle-of -Arrival Tracking: One way to reduce the required 
number of separate detectors (and resultant processing), for a given receiving aperture, 
would be to track the instantaneous angle of ai'rival of the received field. This type of 
operation has previously been discussed for a single detector and is extended here to 
the array case. If the total receiving aperture is not too large, so that the major 
phase distortion over the entire aperture is a tilt, the number of detectors required in 
the array can be reduced by tracking the phase-front tilt with the entire array. The indi- 
vidual detectors can now be approximately three times larger in diameter without 
significantly degrading the performance. The relatively large phase-coherence diam- 
eters at 10.6 pm might make this system attractive at that wavelength. If the receiving 
aperture is too large (>3 or 5 r Q or some other criterion), the phase distortion over the 
entire aperture is poorly approximated by a linear tilt and the idea of tracking with the 
entire array is no longer attractive. In this case the tracking would have to be done by 
individual detectors (or subarrays). Note that phase-front tracking may reduce the 
required number of detectors by a factor of approximately ten. 

It should be emphasized that the over- all system that is being discussed here is 
essentially a cascade of a tracking system(s) and an array-processing system. Because 
of this, the comments on phase-front tracking and array processing made in the two pre- 
ceding sections are directly applicable. Also, as with the case of a single tracking 
detector, no direct performance calculations are available, but the performance can 
perhaps be inferred from the static- array case; the conjectures required for this to 
be the case are similar to those for the single tracking detector. 

Array with Channel Measurement: If the time-space variations of the amplitude 
and phase of the received field were known, the optimum receiver would use this knowl- 
edge to help it make its decisions. For the receivers discussed above no attempt has 
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been made to obtain this type of channel knowledge with the exception of the angle- of- 
arrival tracking receivers, which do measure and exploit a knowledge of the field angle 
of arrival, although not directly in the decision process. Thus a logical extension of the 
array receivers discussed above is a channel-measurement receiver which attempts 
to measure both the amplitude and phase associated with each coherence area in the 
receiving aperture, and uses this measurement in its decision process. Such a pro- 
cedure is attractive for the turbulent channel because the amplitude and phase changes 
associated with each spatial diversity path are very narrow-band processes (<1000 Hz). 
The narrow-band property means that accurate measurement of these processes should 
be feasible. The gain of this receiver over the simple array receiver is that the deci- 
sions can be based on channel measurements made over times long relative to the basic 
decision time for the signal set that is being used. If this decision time approaches 
the channel coherence time (=1 msec), then the measurement receiver will not provide 
any gain; on the other hand, if high rates and reasonable signal set sizes are to be used, 
the decision time will be short relative to the channel coherence time. This type of 
receiver also exploits more fully the instantaneous properties of the phase process than 
was the case with the array receiver which assumed phase independence from coherence 
area to coherence area; this latter assumption is clearly not true and leads to a subop- 
timal receiver. It should also be noted that the angle- of- arrival tracking and channel 
measurement ideas can be used simultaneously; again this would reduce the number of 
detectors required. 

One implementation of the type of channel-measurement receiver described above 
is considered in Appendix I. The required processing is developed as an extension of 
the array receiver previously described. There are many open questions about the exact 
receiver structure and the signaling schemes that are appropriate for this_application. 

In particular, the signals should be such that the channel measurement can be effectively 
carried out without dependence on the signal sequence that is actually transmitted. As 
pointed out in Appendix I, the relative stability between the transmitting and local- 
oscillator lasers will be a crucial parameter, at least in the system implementation dis- 
cussed there. This requirement may, however, be somewhat relaxed by the use of one 
or more phase-lock loops, and other realizations of the basic measurement idea such 
as those discussed in section 4.3.3 essentially avoid this particular difficulty. 

The performance of such channel measurement receivers is unknown at the present 
time. The evaluation of this performance will involve both questions of how well the 
channel parameters can be measured and how well the receiver can make its decisions 
when using noisy channel estimates. 

Signal Design: The signal design question for heterodyne reception essentially 
reduces to the signal design question for radiofrequency channels in which the chan- 
nel is slowly fading and the noise is additive, white, and Gaussian. The only dif- 
ferences are that the fading is log-normal and the phase Gaussian. If, however, a large 
phase variance is assumed, the channel behaves as though the phase had a uniform phase 
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density. As well as the normal modulation techniques, polarization modulation is also 
readily available at optical frequencies. If it is used, care must be taken in the receiving 
optics and heterodyning system to ensure that no incidental polarization changes are 
introduced. 

Experiments: As with direct detection, the form and advisability of some hetero- 
dyne system experiments must await further theoretical work, for example, the mea- 
surement receiver. There are several experiments, however, which can and should be 
performed soon. These are concerned in part with the statistics of the resulting 
fading signal after heterodyning. In particular, the fading statistics associated with a 
single detector, as a function of receiving aperture diameter, should be measured. This 
information will assist the evaluation of the performance of such receivers when the 
receiving aperture exceeds a coherence diameter. Similar measurements of the fading 
statistics should be made for a phase-front tracking receiver. These measurements 
are needed before the static and tracking receivers can be compared with any con- 
fidence. It should be emphasized that the statistical measurements, probability density, 
and perhaps joint density which are being discussed here would provide much more 
information than the few signal-to- noise ratio measurements available, at present. Also 
measurements should be made to evaluate the performance of heterodyning arrays. 

4.3.3 Avenues for System Development 

Several means of improving optical system performance, over that previously con- 
templated, have already been discussed, for example, proper use of the spatial diversity 
within an aperture, and estimation of channel parameters for use in tracking or adap- 
tive receivers. Another means of improving system performance may be to employ 
an adaptive transmitter, that is, a transmitter that measures the state of the channel 
and bases its operating strategy upon these measurements. 

Several kinds of transmitter adaptation are of interest. For example, one might 

attempt to estimate the path loss and adjust the transmitted data rate to ensure a given 

level of performance. Or, one might attempt to estimate the spatial phase degradations 

that the atmosphere will introduce, and predistort the transmitted phase front to com- 

1 31 

pensate for these distortions. .The latter adaptation, at present, is somewhat specu- 
lative and falls in the category of advanced research; the former adaptation has been 
used on some nonoptical channels and is better understood. We now briefly discuss it 
further. [These ideas were introduced at the Workshop by P. E. Green.] 

Variable-Rate Communication: Variable -rate communication is of interest because 

the amplitude fluctuations on the atmospheric channel are log-normal and slowly 

varying. The log-normal distribution is very skewed and values that far exceed 

the mean are not improbable. Astronomers find that with small apertures sudden 

intensity increases of up to 400% above the mean level, with durations of a hundreth 

48 

of a second, are not unusual. The objective of variable-rate communication is to 
exploit these periods of high channel gain. 
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To be specific, we consider a ground-to- space communication system. The space 

receiver will normally observe only a single coherence area, so it is reasonable to try 

to match the transmission rate at the earth terminal to the instantaneous channel gain. 

❖ 

The improvement is large for Rayleigh fading and should be even greater for the more 
asymmetrical log-normal form of fading. Of course, the use of a variable-rate system 
requires that the channel gain be monitored. For the ground-to- space channel, under 
the assumption of instantaneous atmosphere reciprocity, this gain can be monitored by 
measuring the gain of a signal coming from the spacecraft. (The concept of instan- 
taneous atmosphere reciprocity is considered in Appendix G.) The effect of the point- 
ahead angle requirement may limit the utility of the technique at visible frequencies. It 
does not appear to be a problem at the far infrared (10.6 p.), however, even at syn- 
chronous altitudes (see D. L. Fried, Appendix H). 

The variable -rate communication system, while an attractive idea, must be 
examined further to establish its potential and applicability. In addition to the reciproc- 
ity concept, the system-performance improvement, in terms of rate gain or power 
savings, and the associated queing times, with the resultant memory requirements, 
require investigation. 

Transmitted Reference Systems: Transmitted reference techniques provide another 
promising approach to the problem of channel measurement and receiver adaptation. 
Briefly stated, one first observes that, if the spatial phase of the received field were 
known, a near optimum receiver would heterodyne the received field with a local- 
oscillator field possessing the same spatial phase. One then observes that the required 
"local- oscillator " field can be obtained by transmitting a reference, or pilot tone, along 
with the signal. Precisely stated, if the reference and signal fall within a coherence 
bandwidth, they will both experience the same spatial phase distortion in the channel 
and hence the received "pilot tone" will have the desired spatial phase and may be used 
as a "local oscillator." 

The difficulty with the transmitted reference scheme is that a strong local oscillator 
can be attained only at the expense of signal power. An appealing variation on the 
scheme employs an optical amplifier to bring the level of the received pilot tone to the 
desired value. This amplifier must, however, preserve the spatial phase of the pilot 
tone; that is, it must be an "image amplifier." Such a system is discussed further in 
Appendix I. Its feasibility should be further investigated. More generally, the ques- 
tion of optimum quantum receivers for transmitted reference systems should be 
explored. 
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For example, to acheive a 10 error probability one would have to use 47 dB or more 
average transmitter power with a fixed rate system than with a variable-rate system 
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whose bit rate was varied to maintain constant received signal-to-noise ratio. 



V. Scattering Channels 


One objective of the Workshop was to determine the questions that must be answered 
and the data that must be obtained to assess the feasibility of optical communication 
through clouds, fog, and haze. Since this possibility has not been widely recognized, the 
questions that emerged can only be viewed as an initial attack on the problem. If the 
answers to these questions are encouraging, more comprehensive investigations will 
be in order. 

The common characteristic of channels containing clouds, fog, and haze is that they 
involve an appreciable amount of scattering of the incident radiation. Although many of 
our comments apply to all of these channels, we shall frame our discussion in the con- 
text of water clouds. 

That the utility of an optical communication link is reduced by the presence of clouds 
is well known. Indeed a system that is designed to operate efficiently in their absence 
will often be inoperable in their presence. Preliminary results suggest, however, that 
useful communication from space to ground through clouds may be possible with a prop- 
erly designed system. In particular, the effect of clouds upon the optical field at a 
receiver is grossly similar to the effect of the turbulent atmosphere; in both situations 
the field is distorted in time, frequency, and space. The magnitudes of the effects differ 
appreciably, however. For example, the turbulent atmosphere may spread the beam- 
width of a highly collimated pencil beam to tens of microradians — a thick water cloud 
may spread it to tenths of a radian. Also, the coherence bandwidth of the atmosphere 
appears to be at least tens of GHz — that of a thick cloud may be of the order of MHz. 
Finally, the fading statistics of the field received through the clear turbulent atmosphere 
are probably log-normal; for a thick cloud they are probably Gaussian. 

5. 1 COMMUNICATION CONSIDERATIONS 

The differences just discussed are important; but the gross similarity between cloud 

and atmospheric channels suggests that the potential for communication through clouds 

is primarily limited by the ratio, P r /N Q , of the total average collected signal power to 

the noise power per Hz that would be collected if the receiving aperture were diffraction- 
111133 

limited. ’ This conclusion is perhaps even better established for the cloud chan- 

nel than for the turbulent atmosphere because channels with Gaussian statistics have been 

134 34 135 

more thoroughly investigated than have those with log-normal statistics. ’ ’ 

In particular, if large bandwidths and high peak powers can be employed, the system 
capacity is just 1.44P r /N Q bits/sec; that is, the system error probability can be made 
arbitrarily small (at the expense of system complexity) for any information rate less 
than 1.44P r /N Q bits/sec. This capacity is, of course, just that which would be obtained 
in the absence of a cloud for the same ratio of total average collected received power to 
noise power density per Hz per diffraction-limited antenna. 
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Fig. 9. Model for propagation through 
a thick cloud. 


Given the preceding conclusions, the next 
question is, To what extent does the pres- 
ence of a cloud reduce the power that can 
be collected by a given receiving aperture ? 
The answer, at present, is not precisely 
known, but in the absence of absorption the 
amount of power collected by a receiver 
located "near" the cloud is often at least a 
few per cent of that which could be collected 
in the absence of the cloud. To see how this 
can occur, we consider the particular case 
of a thick nonabsorbing cloud. 

As illustrated in Fig. 9, we suppose 

2 

that a beam of power density p watts/m 
and area is normally incident upon the 
upper surface of a thick planar cloud. A 
fraction, 1 -t, of the power is reflected and, 
in the absence of absorption, the remaining 
fraction, t, traverses the cloud emerging 
from the bottom or sides. For simplicity, 


we ignore the radiation emerging from the sides. Furthermore, we suppose that the 
cloud is so thick that its lower surface behaves as a Lambertian radiator that is "uni- 


formly" bright over some circular area A The power density from this area is then 

A l t 2 

“A - p W/m . 

2 

Now suppose that the receiving optics has a circular collector of area A f parallel 
to the cloud surface and a field of view with cone half-angle <j> which is completely filled 


by the area A. The power collected from an incremental surface area ds of the cloud 
will then be 


, , .,2 A, ptA cos 0 

1 / cos 6 \ P_ r 

2ir\ 2 / A 2 


ds, 


(36) 


where 0 is the zenith angle to the element. If the axis of the optical system passes 
through the center of A,,, the total collected power will be 

_ A A pt 

P r = a [l — cos <j)], (37a) 

2 

provided that 
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For purposes of comparison, the power collected in the absence of the cloud would 
be 



if the transmitted beam did not diverge appreciably in propagating the distance h. Note 
that the ratio of the powers is 

P At 

~ =-r^-[l - cos 4>]. 

P r A 2 

In the particular case of an infinite plane wave incident upon the cloud (so that = A = 
oo) and a receiver with a 2ir sr field of view (so that 4> = 90°), P r /P r is approximately 
equal to t. Since t often exceeds 10%, we conclude that, subject to our many assump- 
tions, the presence of the cloud would reduce the collected power by at most a factor of 
10. Although such a loss is certainly significant, it may be acceptable in some circum- 
stances. 

Of course, the loss increases, for a given field of view, as the distance from cloud 
to receiver increases and also as absorption in the cloud increases. Moreover, it varies 
with wavelength and cloud type. Thus the crudity of the preceding arguments cannot be 
overemphasized; however, they do suggest that the possibility of optical communication 
through clouds should be further investigated. 

Our ignorance of the channel characteristics is so great that it would be premature 
to discuss system structure in any detail, but a few general observations are in order. 
Since very small coherence areas are expected for clouds located near the receiver, an 
efficient receiver must properly exploit the " diversity" contained in its aperture just as 
a system must for the turbulent atmosphere. In fact, the entire discussion of receivers 
for the atmosphere is relevant to clouds, although the reduced coherence bandwidth 
(increased time dispersion) of clouds may limit the utility of transmitted reference sys- 
tems. Also, because of this dispersion, special signaling techniques may be necessary 
when operation at high information rates is required. 

A practical difference between the use of atmospheric and cloud channels relates to 
the field of view requirements 'and the size of the coherence areas. As noted previously, 
the receiver field of view may have to be much larger for a cloud channel than for an 
atmospheric channel — perhaps tens of degrees. Also, the coherence areas of a cloud 
channel will often be much smaller than those of an atmospheric channel. Consequently, 
the already difficult problems of implementing a near-optimum atmospheric system 
become even more pronounced for a cloud system — at least when the cloud is near the 
receiver. 

In summary, efficient systems for operating through clouds may be feasible, but 
they will probably require the development of new devices. These devices are sim- 
ilar in kind to those needed for atmospheric systems, but the requirements are more 
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severe. Based upon crude extrapolations of available data, it appears that such sys- 
tems may yield adequate performance in some situations. A thorough investigation 
of the channel characteristics is, however, a prerequisite to any serious contempla- 
tion of its use. Some of the channel characteristics to be investigated will now be 
considered. 

5. 2 CHANNEL CHARACTERIZATION 

To fully assess the communication potential of scattering channels, it is necessary 
to know the statistical properties of the field in the receiving aperture for the modulation 
scheme of interest. Although it has not been verified, the multiple scattering that occurs 
in a thick cloud will probably result in a received field that behaves as a zero-mean 
Gaussian random process. Assuming this to be true, the channel characterization prob- 
lem reduces to that of determining the correlation function of the received field. Although 
complicated by the spatial fading that occurs, this problem is similar to that encountered 
with other fading dispersive channels; hence, we shall not pursue the general approach 
to the problem here. 134, 34, 135-138 

A complete determination of the correlation function referred to above is desirable, 
but experience with other channels has shown that its detailed behavior does not greatly 
influence the communication performance levels that can be realized. In fact, a set of 
approximately four parameters usually provides sufficient information to grossly esti- 
mate the performance levels that can be achieved. Were it not for the spatial distortion 
of the wavefront, three parameters would suffice: the total average received power col- 
lected by the aperture, the average time dispersion suffered by a short pulse, and the 
average frequency dispersion suffered by an unmodulated carrier. In the presence of 
spatial distortion, an estimate of the number of coherent areas contained in the receiving 
aperture is also required. Stated alternatively, one must know the number of essentially 
independent plane waves that are required to characterize the field over the receiving 
aperture, that is, the angular spectrum. These parameters, of course, are just those 
discussed in connection with the turbulent atmosphere. 

Clearly, the complete correlation function of the field, on the one hand, and the four- 
parameter description, on the other, represent extremes of channel characterization. 
Between them lie a variety of compromises in completeness and complexity. One char- 
acterization that may be particularly useful in the development of efficient systems con- 
sists in specifying the absolute propagation delay and Doppler shift, as well as the time 
and frequency dispersion, associated with each of the plane waves that are contained in 
the angular spectrum of the received field. 

The development of adequate characterizations for cloud channels will involve an 
application of scattering theory. We now briefly review the state of that theory, the 
information about the cloud that is required to apply it, and the physical properties of 
clouds, fog, and haze. 
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5. 2. 1 Scattering Theory 


The scattering and absorption of a plane, monochromatic wave by a homogeneous 

139 

sphere can be described exactly by using a theory developed by G. Mie. The theory 
yields absorption and scattering cross sections for arbitrary polarization of the incident 
beam, as well as the distribution of the light scattered out of the beam, that is, the scat- 
tering function. Efficient computer programs for evaluating the relevant expressions 
have been written. The cross sections depend on the refractive index of the particle with 
respect to the surrounding medium, the wavelength of the incident light, and the particle 

diameter. The more difficult case of many-layered spherical droplets with different 
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refractive indices has been solved by several investigators. ’ The resulting equa- 
tions can also be programmed, with additional computer time being the main price paid 
for the additional complexity. For more realistic models of atmospheric particles (that 
is, irregular, possibly many-layered, particles), the scattering properties are not cal- 
culable, but they can be measured. 

The single -particle properties can be integrated for a particle size distribution (if 
known) to arrive at the monochromatic volume scattering properties of the medium. Thus 
for an optically thin medium with plane -parallel boundaries in which' only single scat- 
tering is important, the reflected and transmitted intensities are simple functions of the 
medium's optical thickness and geometry. The optical thickness of a medium along an 
incremental path length ds is d-r^ = P ex ^. ds, where P ex ^. is the volume extinction coef- 
ficient, the sum of the volume scattering and absorption coefficients (|3 = P sc + P a bs^‘ 

For optically thick media, in which multiple scattering must be considered, the 
transfer problem is extremely complicated. Except for a few simple cases, solutions 
to the exact equations of radiative transfer are difficult to obtain. The problem becomes 
even more complicated when polarization is taken into account because the light beam 
must be expressed in terms of the four-component Stokes vector notation and the scat- 
tering phase function becomes a four-by-four matrix. In general, polarization should 
be considered, since an initially unpolarized beam will be partially polarized after 
a single scattering event. Its importance in communication applications is still not 
clear. 

At least two alternative means of approaching the multiple-scattering problem exist. 

One means consists in modeling the scattering events by a random-walk, or Monte Carlo, 

calculation. This method can be applied to nonhomogeneous media and is not limited to 

plane -parallel boundaries. Polarization can be included in the model quite easily; also, 

142 143 

nonstationary problems can be solved. Several authors ’ have used the Monte 

Carlo method to calculate multiple-scattering effects in haze and clouds. Their results 

for multiple scattering in a plane -parallel, Rayleigh (or molecular) atmosphere compare 

144 

very well with the exact results of Chandrasekhar and of Coulson, Dave, and Sekera. 

Another means of approaching the multiple-scattering problem is to model a thick 

138 

cloud as a concatenation of thin clouds within which only single scattering occurs. 



72 


OPTICAL SPACE COMMUNICATION 


This approach involves approximations that preclude high-precision results, but it has 
the advantage of yielding closed-form results that may be suitable for communication 
system analysis and design. Further correlation of the results obtained from this 
approach with those obtained by Monte Carlo techniques and by experiment will be 
required to establish their suitability for communication applications. 

In summary, several techniques for predicting the transmitted (and reflected) inten- 
sities from haze, fog, or clouds are at hand. Also, the scattering properties of spher- 
ical particles can be described exactly if we know the particle size distribution, and the 
real (n) and imaginary (n 1 ) parts of the index of refraction. More information on cloud 
particle composition, size distribution, ana particle shape may be required to apply 
these results to real clouds. Alternatively, it may be worth while to try to assess the 
significance of these unknown parameters by examining extreme cases covered by 
existing models. If the results indicate a lack of sensitivity to them, our present knowl- 
edge of cloud characteristics may suffice. This state of knowledge will now be briefly 
summarized. 

5. 2. 2 Physical Characteristics of Clouds, Fog, and Haze 

Clouds are assemblies of ice crystals and/or water droplets with sizes ranging from 
1 to 50 pm. The liquid water content of clouds varies from 0. 1-1.0 gm/m for clouds 
above 2 km, to as much as 10 gm/m for low-level clouds. Detailed knowledge of clouds 
is incomplete, but certain general characteristics are known. At altitudes where the 
air temperature is less than -40 °C, clouds in thermal equilibrium with the surrounding 
air are composed of ice crystals rather than water droplets. Samples show that the ice 
crystals in cirrus and stratus clouds are predominantly prismatic columns with hexag- 
onal bases averaging 40 pm in width and 100 pm in length. The particle density in 
stratus and cirrus clouds is small (0. 1-30 particles/cm ). There is no information 
about the properties of high-altitude cumulus clouds. Low-altitude cumulus cloud 
samples indicate that they contain many more particles (100 particles/cm ) and that the 
droplets are smaller, ranging in diameter from 5. 0 pm to 30 pm. On occasion, these 
lower clouds can be composed both of ice crystals and water droplets. The cloud struc- 
ture is far from homogeneous, and the cloud properties can vary considerably from point 
to point in the cloud. 

Haze and fog, being more accessible to ground-based instruments, have been 

studied extensively for many years, with a large variety of instruments and techniques. 
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According to Rozenberg ’ haze and fog are characterized by smaller particle 
sizes. Haze is composed of finely dispersed formations of particles with diameters that 
are less than 1. 0 pm and with an average particle size that is close to 0. 25 pm. Dry 
haze is observed at relatively low humidity (< 80%) , foggy haze at moderate to high 
humidity (>85%), and true fog is observed only at saturation. Roughly stated, for dry 
haze the visibility (or meteorological range) is greater than 3 km; for foggy haze, 
visibility is between 1 km and 3 km, and for true fog the visibility is less than 1 km. 
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Haze and foggy haze are relatively stable atmospheric states. The transition from 
haze to foggy haze is practically continuous, and is usually accompanied by the disap- 
pearance of relatively fine particles and the appearance of larger particles (presumably 
by particle growth and agglomeration). 

5. 2. 3 Complex Refractive Index for Water and Ice Particles 

The complex refractive indices of water and ice play an important role in scattering 

theory. The real and imaginary parts of this index are denoted n and n' , respectively. 

147 

The values of n and n 1 are known "reasonably" well for liquid water. Recent 

experimental results at 10. 6 pm suggest, however, that the tabulated values may not be 
sufficiently accurate for path loss estimates at that wavelength, because of the strong 
influence of n' on the results of theoretical calculations. Direct measurement of n 1 is 
difficult, but an approximate value can be obtained by measuring the intensity of the 
radiation scattered at several well-chosen angles by an artificial water cloud (see sec- 
tion 5. 3). The uncertainties in the published values of n and n 1 for ice are probably 
greater than those for water, so theoretical calculations of scattering by ice clouds are 
questionable. 

There is some experimental evidence that relatively thin ice clouds absorb appre- 
148 

ciably. Results of experiments on the reflection and transmission of ice clouds and 
frost layers indicate that near 2. 8 pm the reflectance is low and the transmittance is 
greater than at neighboring shorter or longer wavelengths. This phenomenon, known as 
the Christiansen filter effect, is caused by a decrease in n on the short wavelength side 
of a strong absorption band. Transmission measurements on artificial ice clouds would 
determine if the transmittance near 2. 8 pm is great enough to be of value in a communi- 
cation system. The possibility of increased transmittance on the short wavelength 
side of the strong absorption band near 13 pm should also be investigated. Since the 
Christiansen filter effect requires a rather sharp absorption band and the liquid water 
absorption bands do not appear to be as sharp as the bands of ice, a useful Christiansen 
effect in a water-droplet cloud is less probable than in an ice cloud. 

5. 3 PROPOSED EXPERIMENTS 

We now describe experiments that would provide some of the required information 
just discussed. The proposed experiments are relatively simple and are designed with 
two objectives in mind. The first objective is to obtain enough data on the scattering 
and absorption properties of real clouds to determine whether optical communication 
through clouds is feasible, and if so, to determine the order of magnitude of the attenu- 
ation to be expected from various types of clouds. The second objective is to obtain some 
of the information about n, n 1 , and cloud composition that is now lacking. No advances 
in the state of the art are required to perform the experiments. 


74 


OPTICAL SPACE COMMUNICATION 


5.3.1 Laboratory Experiments on Artificial Clouds (Ice and Water) 

The scattering phase function for various sized single spherical particles of water 
and ice should be measured at several wavelengths and all scattering angles. This will 
provide information on n and n 1 and will give a check on the theory for single particles. 
Additional experiments should be performed on single ice crystals of different sizes and 
shapes in order to obtain information about their scattering properties relative to 
spheres. 

The extinction and scattering at several angles by artificial clouds of water droplets 
and ice should be measured at several wavelengths in the atmospheric windows, partic- 
ularly in regions where one might expect the Christiansen filter effect. The results will 
provide much needed information on n and n' and a check on the theory for a wide range 
of optical thicknesses. 

5.3.2 Natural Cloud Experiments 

One experiment of this kind is designed to measure the insertion loss through typical 
clouds. Specifically, a heliostat, or other sun-tracking device, should be used to form 
an image of the sun on a monochromator, or other bandpass filter. The major atmo- 
spheric windows from the visible to 25 pm should be included, and the field of view of 
the instrument should be smaller than the sun. The radiant intensity in the direction of 
the sun would be measured through a clear sky and through a variety of cloud cover. 
Measurements should also be made while viewing clouds within 10° or 15° of the sun. 

The experiment will provide information on the attenuation coefficients at different 
wavelengths. Although the cloud composition will not be known accurately, it may be 
possible to establish approximate cloud conditions under which optical communications 
are feasible. By making some realistic assumptions about cloud composition, rough 
checks on the theory and on the values of n and n' can be obtained. 

Another useful experiment would illuminate a cloud with a narrow beam. The energy 
source could be a visible wavelength laser, mounted in an airplane or on a mountain top. 
Receivers on the ground would record the signals received through the cloud. Given this 
configuration, or its converse, the following measurements suggest themselves. 

1. With a cw source, the average energy distribution on the underside of the cloud 
should be mapped, using a scanning energy- detecting receiver and several receiver 
fields of view. The fluctuations in energy should also be examined. 

2. Again with a cw source, using two energy- detecting receivers, the correlation 
of the energy from two small areas on the underside of the cloud should be studied as 

a function of the spacing between these areas. For this measurement, the area should 
be as small as possible; that is, the receiver field of view should be as small as is con- 
sistent with usable signals. For thick clouds the expectation is that the two energies will 
be uncorrelated if the areas do not overlap. 
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3. With the energy source pulsed, and two energy-detecting receivers, the time- 
dispersion effects of multiple scattering through the cloud should be measured. As a 
very rough estimate, one would expect multipath dispersion through a 1-km thick cloud 
to be between 0. 3 psec and 3 psec. On this scale, a source pulse of 100 nsec duration 
with receiver time resolution of 200 nsec would be adequate. 
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Abstracts of Briefing Sessions 
NASA Optical Communication Program 
H. L. Anderton 

Two of the responsibilities of NASA Headquarters are(i)to develop a coordinated pro- 
gram using the best talent in the field centers, and (ii) to obtain financial support from 
the United States Congress for the program. To get the needed financial support requires 
a careful presentation of the technical work of the field centers in language that makes 
sense to the nontechnical members of Congress and still gives an accurate picture of the 
work that is to be done. The condensation of material that results from these attempts 
sometimes leads to extremely simplified statements that might lead-a dedicated scien- 
tist to wonder whether the work was being adequately covered and, especially, whether 
the proper caveats were being included. (Some examples of Congressional presenta- 
tions were shown.) 

In spite of cutbacks in the over-all NASA budget for 1969 and subsequent years, there 
are indications that the optical communication technology area will be given increasing 
emphasis and funding. The rationale for this is the great promise afforded by the laser 
for increased data rates in space communication and the need for a large Space telescope 
for astronomical observations. (Typical funding figures were shown.) 

High-Capacity Optical Communication System 
T. Kinsel 

The large potential bandwidths available at optical frequencies suggest the possibil- 
ity of an optical communication system capable of operating at very high data rates. 
Active optical devices that are available at present have relatively limited bandwidths, 
however. The purpose of this session is to describe a design of the transmitting and 
receiving terminals that circumvent the usual limitation of available optical devices. 

The basic single -channel system operates in the binary pulse code modulation (PCM) 
mode with the binary states represented by the presence and absence of optical energy. 

A phase-locked laser driven by an external clock frequency is used as the source and 
an optical modulator, operating as a gate, transmits or rejects individual optical pulses. 
The performance of an experimental system utilizing a helium-neon laser operating at 
a bit rate of 224 M bit/sec, with the modulator constructed of lithium tantalate, is 
described. 

A first stage for increasing the information capacity of the system uses time domain 
by time division multiplexing. Since the optical pulses from a phase-locked laser are 
generally quite narrow relative to their period, a number of independent channels can 
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be interleaved in time -4 if a Helium-Neon laser is used, and 32 if an Nd:YalG laser is 
used. At the present time, higher ultimate data rates are obtainable if this process is 
done in the optical domain, by using a device called a sinusoidal polarization switch. This 
device is constructed with an electro-optic element driven by a sinusoidal timing signal, 
and systematically varies the states of polarization of a series of pulses which is coupled 
through it. This device is used in both multiplexing and demultiplexing, and the design 
of a 24-channel time-multiplexed system will be sketched. Although the data rate of such 
a system is approximately 6. 7 G bit/sec, each active optical device is operating within 
its bandwidth. The performance of a two-channel time-multiplexed system using a 
helium-neon laser will be described. 

A further stage of multiplexing may be done in the frequency domain, by using tech 7 
niques of frequency division multiplexing. Here, independent transmitting terminals, 
which may be time-multiplexed, would be constructed at different carrier frequencies, 
by using different laser sources or by parametric oscillation. The outputs would 
be combined into a single polychromatic beam. 

A final stage of multiplexing may be done in the spatial domain, by using a multi- 
mode beam waveguide. One way of doing this would be to launch a number of Gaussian 
beams into the waveguide with different launching parameters. At the receiver plane 
these would be resolved as separate beams. 

Combining all three stages of multiplexing simultaneously indicates the possibility 
of constructing, with available or nearly available devices, an optical communication 
system operating at 10 bits/sec. 


Direct-Detection Communication Systems 
J. R, Kerr 

Direct-detection optical communication systems will be signal-shot-noise limited 
when large (bandwidth X SNR) products are required, provided that detectors -with-gain 
are available and background levels are not too high. Thus, with the development of 
microwave -bandwidth photomultipliers and avalanche photodiodes, direct detection has 
a clear advantage over heterodyne systems at visible and near- infrared wavelengths. 
This is especially true when system complexity and atmospheric effects are considered. 
For instance, lasers operating in the FM mode comprise ideal transmitters for such 
systems, and single-frequency operation is unnecessary. Alternatively, mode-locked 
lasers represent simple pulsed sources for optical PCM applications. Studies show that 
the most attractive system alternatives, at present, are (i) direct detection with doubled- 
YAG or Argon laser transmitters, and (ii) heterodyne detection at 10. 6-p wavelengths 
with stabilized CO,, lasers. We must recognize, however, that (SNR bandwidth) char- 
acterizations of these alternatives may be inadequate from a communication theory 
standpoint, and more sophisticated configurations may evolve. 
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Optical Heterodyne Communications Experiments at 10. 6 pm 

F. E. Goodwin 

A 10.6-pm optical heterodyne communication system that is in operation over a 
30-km (18 mile) path was described. Television signals have been transmitted over the 
beam by the use of optical FM, achieved with internal electro-optic phase modulation 
of the transmitter oscillator. A frequency deviation of 3 MHz is thus obtained for a 
driving voltage of 1 kV, as compared with the fundamental frequency noise of the system 
of ~1 kHz. For clear-day operation, average carrier-to-noise ratios at the receiver 
approach 60 dB, and demodulated video signal-to-noise ratios average 50 dB. 

Optical propagation experiments reveal amplitude scintillation noise in excess of 
20 dB, represented by fluctuations of +10 dB and -10 dB from the average carrier level. 
Other observations show that the 10. 6-p.m beam has a remarkable ability to penetrate 
fog. Data related to scintillation, visibility, temperature, humidity, and optical absorp- 
tion were presented. 


Gas Laser Sources 
W. B. Bridges 

In this presentation the salient characteristics and state of the art of Helium-Neon, 

Argon-ion, and carbon dioxide lasers are reviewed. A comparison is made on the basis 

of a figure of merit for ground use of the form (Power )/(Wavelength) and ((Power) X 
2 2 

(efficiency) )/((Wavelength) X (length)) for spacecraft use. The best performance fig- 
ures for the three laser types in multimode, TEM Q0 , and single-frequency operation are 
given for laboratory, commercial, mil-spec and space-qualified units. Specific examples 
of lasers under development at the Hughes Aircraft Company are described: a space- 

qualified Helium-Neon laser, an airborne Argon-ion laser, and a high-power pulsed 
ground-based CO., laser oscillator-amplifier system. 

Solid-State Lasers 
J. E. Geusic 

The solid-state laser is now quite competitive with its gaseous counterparts. Con- 
tinuous operation at output powers over 100 W has been achieved with Nd:YalG at effi- 
ciencies as high as 1. 5%. Laser efficiency and life (1000-5000 h) is limited by the 
tungsten-iodine pump lamps in present use, and it is expected that employment of lumi- 
nescent semiconductor pumps now under development will raise the efficiency by a factor 
of at least five and permit even longer life operation. The solid-state Nd:YAlG laser is 
quite versatile in operation because, in addition to cw operation, it can be operated in 

a continuously pumped repetitively Q-switched mode (peak power >1 kW, average power 

3 4 

>1 W, pulsewidth 0. 1 psec, pulse repetition rate 10 -10 p. p. s) or a mode-locked mode 
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(pulsewidth 30 X 10 sec). The Nd:YAlG laser output, which is at 1. 06 p, can be con- 
verted with 100% efficiency to 0. 53 p by utilizing the nonlinear barium sodium niobate 
(Ba_NaNb O ) inside of the laser cavity. The use of barium sodium niobate with the 
YAlG laser provides a practical solid-state green coherent source. Wavelength tun- 
able coherent sources are now under development, with barium sodium niobate used 
as the nonlinear element of a parametric oscillator. Such a continuous, tunable 
oscillator has been built in the near infrared, with the doubled YAlG green source 
as pump, and oscillation from 0. 95 p to 1. 15 p has been demonstrated. Continuous 
parametric oscillators covering the entire visible and near infrared are now prac- 
tically achievable. 


Optical Modulators 
J. E. Geusic 

Various physical phenomena can be utilized for modulating an optical beam. The 
electro-optic effect, the magneto-optic effect and the acousto-optic effect are three phe- 
nomena that have been seriously considered for high-speed optical modulation. Of these 
the electro-optic effect has been studied and developed most extensively, and there are 
good theoretical reasons to believe that in principle it has the best potential. For 
example, in comparing the electro-optic effect and the acousto-optic effect it is evident 
that the electro-optic effect has greater potential, since the velocity of light is approx- 

5 

imately 10 times greater than the velocity of sound. 

The analysis of electro-optic modulator devices is quite straightforward. One result 
of this analysis is that (in an optimally designed modulator) electrical power that is 
required to modulate light is proportional to the modulation bandwidth B, proportional 
to the wavelength cubed, and inversely proportional to the square of the electro-optic 
coefficient (under the assumption of a linear effect). For most electro-optic materials, 
the electro-optic coefficient is only a weak function of wavelength over the visible to near 
infrared, and one general conclusion is that it is easier with most materials to build 
modulators in the visible range than it is at longer wavelengths. At present, the best 

electro-optic materials are the niobates (LiNbO^, Ba 2 NaNbg0^g ). With these 

materials, modulation bandwidths of approximately 500 MHz have been achieved in 
the visible, with the expenditure of ~1 mW/MHz of bandwidth. The electrical Q's 
of the materials are generally 100 or less, and this imposes a practical maxi- 
mum in the power that can be delivered to the modulator, for the finite Q's 
cause heating effects that produce optical nonuniformities. At present, modulation 
bandwidths of several GHz in the visible are not practically achievable with pres- 
ent materials, and in the infrared, particularly at 10 p, even 0.5 -GHz bandwidths 
are unpractical. New and better electro-optic materials are definitely needed and 
the improvements can be in the form of higher electro-optic coefficients and/or 
higher electrical Q's. 
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Tracking in Optical Communication Systems 
E. C. Fraser 

Three different problems are involved in the operation of an optical communication 
system: pointing — aiming on the basis of a priori information before acquisition is 

accomplished; following — maintaining the sensitive axis aligned with the direction of the 
incoming signal; and data collection — for orbit or trajectory determination. To accom- 
plish these tasks, several fundamentally different system configurations have been 
considered. An open - loop system points each of the terminals on the basis of externally 
supplied pointing information. Such a system is therefore always operating in an acqui- 
sition mode and, therefore, cannot take advantage of techniques such as beam narrowing 
following lock-on. In a closed-loop system, once acquisition has been achieved by open- 
loop pointing, each terminal determines the pointing error of the other terminal's trans- 
mitter, and relays this information back to institute pointing corrections. Pointing errors 
can be detected by such techniques as receiver arrays or scanning of a spatially encoded 
beam. The principal difficulty with closed-loop systems is their inability to respond to 
rapid disturbances or medium fluctuations that results from the long transport delay 
associated with transmission paths of interest (approximately 1 /3 sec round trip to a 
synchronous satellite and approximately 16 min round trip to a spacecraft in the vicinity 
of Mars). In a cooperative system each terminal's transmitter is used as a beacon at 
which the other terminal can aim. Thus all tracking by each terminal is done only at 
that terminal, and maximum tracking performance can be achieved. It appears, there- 
fore, that for most applications a cooperative -system configuration will be most appro- 
priate. 

Modern control theory dictates that for optimum tracking performance of each 
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Fig. A-l. Optimum tracking control system. 
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terminal the tracking control system should take the form shown in Fig. A-l. The pur- 
pose of the estimator is to develop, using algorithms evolving from statistical estima- 
tion theory, minimum variance (least noisy) estimates of the system state (both target 
and telescope) from all measurement data obtained from the system. These state esti- 
mates are then used by the controller to generate optimum (in a minimum mean square 
following error sense) control signals for driving the telescope. The high-speed beam 
steering loop shown following the telescope dynamics must be provided in those cases 
in which the field of view (or beamwidth) is smaller than components of the path distur- 
bances that are too fast for the main telescope to follow. This situation arises, for 
example, if one were to attempt tracking through the atmosphere with a fractional second 
of arc field of view. 

The primary feature of the estimator-controller system configuration is that it 
resolves the dilemma of tracking loop bandwidth selection by separating the tasks of 
tracking -signal processing and telescope -pointing control. To remove line-of-sight dis- 
turbances, the estimator can act as a very lowpass filter, accepting only those signal 
components pertinent to the target trajectory. On the other hand, the controller can inde- 
pendently realize a high-gain wide -bandwidth tracking loop so as to achieve crisp tran- 
sient response and immunity to wind torques and other load disturbances. 

The implementation of the complete estimator-controller system requires consid- 
erable computation capability and complex instrumentation. In many applications the 
system performance requirements are not sufficiently stringent to require this degree 
of sophistication. For these cases there are some simplifying approximations that result 
in systems of considerably reduced complexity. In fact, the conventional autotrack sys- 
tem configuration can be derived from the estimator-controller formulation for the 
appropriate set of systems specifications and approximations. 

Precis of Presentation on Effects of Atmospheric Turbulence 
on Optical Communication 

D. L. Fried 

The following presentation covers a large part of our present quantitative knowledge 
concerning optical effects of atmospheric turbulence, particularly as related to com- 
munication equipment in the space-to-ground and ground-to-space mode of operation. 

The quantitative knowledge that we have is based almost entirely on theoretical calcula- 
tions with a minimum of experimental backup. The theory generally takes the form of 
generalized solutions for instantaneous optical instrument performance in terms of 
instantaneous propagation conditions, from which the statistics of instrument perfor- 
mance are calculated in terms of the statistics of propagation. This is coupled with a 
solution for instantaneous optical instrument performance in terms of instantaneous 
propagation conditions, from which the statistics of instrument performance are cal- 
culated in terms of statistics of propagation. This is coupled with a solution for the 
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instantaneous propagation results in terms of instantaneous turbulence conditions, from 

which the statistics of propagation are computed in terms of the statistics of turbulence. 

The statistics of turbulence, of significance in optical propagation, are derived from the 

Kolmogorov theory of the inertial subrange, which asserts that the refractive-index 

structure function, that is, the mean-square difference of refractive index at two points 

2/3 

at distance r apart, is proportional to r ' . From this, various propagation statistics 
are computed, such as the log-amplitude variance, the phase and wave -structure func- 
tions, the mutual coherence, etc. Fairly good evidence exists for the validity of the 
Kolmogorov theory in the inertial subrange in the atmosphere. Data for substantiating 
this for the lower atmosphere were presented in a series of slides. 

A variety of calculations for log-amplitude covariance and wave -structure function 
for infinite plane-wave and spherical-wave sources was presented in a series of slides, 
all based on the Kolmogorov hypothesis for atmospheric turbulence. The results are 
all derived by using the Rytov approximation. The most pronounced feature of the log- 
amplitude covariance functions is that they go to zero in a distance of the order of the 
square root of wavelength times path length. For vertical propagation down through the 
atmosphere, based on our limited knowledge of the vertical distribution of turbulence in 
the atmosphere, the equivalent path length is approximately 3200 m. 

Next, our understanding concerning the performance of a laser transmitter operating 
in a ground-to-space mode was presented. Calculations, thus far, make available esti- 
mates of scintillation of transmitter antenna gain as a function of diameter. Mean trans- 
mitter antenna gain is available from some work by F. Gebhart to be published shortly 
in the Journal of the Optical Society of America. (Not included in this presentation, but 
a significant result derived during the period of the conference, is the fact that trans- 
mitter antenna gain on an instantaneous basis can be shown to be identical to heterodyne 
receiver gain operating over the same path. Since heterodyne receiver gain is a well- 
studied quantity, this provides us with a substantial body of knowledge applicable to 
transmitter antenna gain.) 

Very simple results for aperture averaging were presented in a single slide. The 
results show that aperture averaging may be considered as the consequence of a set of 
independently scintillating regions on the collection aperture, each region having an area 
of approximately wavelength times path length. 

The problem of general interpretation of wavefront distortion was discussed from the 
point of view of whether or not any of the distortion features have recognizable and poten- 
tially useful aspects. It was pointed out that a large part of the wavefront distortion can 
be interpreted as a linear variation of complex phase. To the extent that the complex 
phase is mostly ordinary phase with very little log amplitude, this corresponds to a 
tilted wavefront, which can be "tracked out," In the absence of log-amplitude variance, 
tracking the tilt can increase coherent aperture diameter by a factor of approximately 
8 dB. In the presence of strong log-amplitude effects, however, the gain is almost neg- 
ligible. 
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Results were presented for long- and short-exposure imaging through the atmo- 
sphere. Pointed out were the possible gain of short-exposure imaging and the general 
tendency of resolution to increase uniformly with diameter until a critical diameter is 
reached, at which point the resolution saturates, that is, stops increasing with increasing 
diameter. The significance of these results is principally in terms of their relation to 
the optical heterodyne receiver, imaging resolution being directly proportional to antenna 
gain and thus to achievable signal-to-noise ratio in a heterodyne receiver. Equivalent 
results showing the same saturation effect were presented for the heterodyne receiver 
average signal-to-noise ratio. The knee of the curve, corresponding to the onset of sat- 
uration, is located at a diameter denoted by r . Graphs presenting values of r Q for var- 
ious propagation paths were shown. Results concerning atmospheric modulation effects 
in heterodyne receivers were presented, indicating the very rapid, very large increase 
in modulation noise as the receiver diameter is increased beyond r . Doppler-shift 
effects were evaluated. It was shown that under all conceivable conditions, Doppler 
shifts were entirely negligible. 

Finally, a series of curves relating binary channel performance in the presence of 
scintillation to performance in the absence of scintillation. The calculations for these 
curves are based on the fairly well-established fact that atmospheric scintillation almost 
always appears to show a log-normal distribution. Significant loss of effective signal- 
to-noise ratio was found for fairly moderate levels of scintillation. 

Effect of Atmospheric Turbulence: 

Experimental Results from the ESSA Research Laboratories 

R. S. Lawrence 

The optical effects of turbulence are measured with mobile equipment so that a vari- 
ety of paths can be instrumented. Paths in the Boulder (Colorado) area ranging in length 
from 25 m to 150 km have been used. Whenever practical, small, high-speed resistance 
thermometers are used to provide simultaneous direct measurements of the magnitude of 
the refractive-index variations. The thermometers have a thermal capacity much lower 

3 

than that of 1 mm of air and a frequency response exceeding 1 kHz. 

Beam spread observed on paths 5-15 km in length is normally ~10 prad at night, and 
often exceeds 100 prad at midday. For short paths, where the turbulence is uniformly 
distributed, intensity scintillations display the log-amplitude variance and covariance 
expected from Tatarski's theory. On longer paths, where the height above ground is great- 
est near the center of the path, log-amplitude variance is reduced and the covariance 
function is distorted in qualitative agreement with the concept of reduced turbulence at 
midpath. Quantitative comparisons with the results from airborne temperature sensors 
are being made to test the usefulness of the covariance function as a measure of the dis- 
tribution of turbulence along the path. 

The saturation effect, where the intensity of scintillations fails to increase with 



APPENDIX A 


93 


increasing strength of turbulence, has been demonstrated to begin at a log-amplitude 
variance of 0. 6 on a 1 -km path, 2 m above flat ground. 

Temporal variation of phase over paths shorter than 100 m is measured with a 

reversible fringe -counting Michelson interferometer. Resulting power spectra contain 

only approximately 1% as much bandwidth as the corresponding intensity fluctuations. 

Measurements of modulation phase over paths of several kilometers show that the usual 

rms fluctuation of 1-sec averages of phase path relative to the 30-sec moving average 
-9 

is found to be 3 X 1 0 of the total path. 

Statistical Communication Theory at Optical Frequencies 
E. V. Hoversten 

This talk had several purposes: (i) to discuss the role of statistical communication 

theory and indicate the modeling questions that must be answered before the theory can 
be applied, (ii) to discuss the differences between the communication models and prob- 
lems at optical and microwave frequencies, and (iii) to illustrate the ideas in the first 
two topics in terms of the turbulence atmosphere channel and discuss some specific 
results for this channel. 

The abilities to establish ultimate performance limits and to provide insight into how 
to build 'good' systems, whose performance approaches these limits, were cited as the 
contributions of statistical communication theory. The types of questions that the theory 
can answer, for example, optimum receiver structure, error performance and capacity, 
and the need, if the full power of the theory is to be realized, for complete statistical 
models, particularly of the channel, background noise, and receiver processing, were 
discussed. The differences between optical and microwave communication that were 
considered included the importance of the spatial aspect at optical frequencies, the back- 
ground noise at optical frequencies, the different noises associated with the detector 
processing, that is, the quantum aspect at optical frequencies, the different channel sta- 
tistics encountered, and the differences in the operations that can be performed with 
available devices. 

An approximate statistical model was discussed for the turbulent atmospheric chan- 
nel and the ideas above were illustrated in terms of this discussion. The form of the 
optimum processing for heterodyne detection with a large receiving aperture was con- 
sidered. This led to a discussion of diversity receivers and the optimum use of spatial 
diversity. The discussion was in terms of digital signals whose duration was short 
relative to the channel coherence time. The error performance that could be achieved 
with the above-mentioned processing for M-ary orthogonal signaling sets was discussed. 
It was indicated that the error probability had the form 

B2~ KE Pr [e] « 2~ KE , 

where K is the number of information digits per transmitted waveform, and E is the 
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system reliability. The dependence of E on the aperture area, turbulence, and energy - 
to-noise ratio per diversity path was discussed. The results which were considered are 
discussed in section 4. 3. 2 of this report. 

Optical Propagation through Clouds 
A. Holland 

The light reflected or transmitted by fog, haze or clouds depends on the number and 
size of the water droplets or particles present, on the absorption and scattering cross 
sections of the particles and on the angle of incidence of the incoming radiation. The 
scattering or absorption cross sections of individual droplets of water can be calculated, 
exactly by using classical Mie theory if the particle diameter and refractive index are 
known. These "monochromatic" single -particle properties can be integrated over a par- 
ticle size distribution (if known) to give the scattering properties of a volume element 
of the medium. For most cases of interest these monochromatic properties can be used 
to represent the medium under moderately broadband illumination (under the assumption 
that the complex index of refraction does not vary rapidly with wavelength). Then for 
optically thin media, in which only single scattering is important, the reflected and 
transmitted intensities are simple functions of the optical thickness and transfer geom- 
etry. 

On the other hand, for optically thick media for which multiple scattering must be 
considered, the transfer problem is extremely complicated. With the exception of a 
few special cases, solutions to the exact equation of transfer are difficult to obtain. It 
is possible, however, to model multiple-scattering problems by using a random-walk 
or Monte Carlo technique. The method has been applied recently to the problem of mul- 
tiple scattering by clouds and fog. 

We present a review of the known properties of haze, fog, and clouds, an outline of 
Mie theory and the equations of radiative transfer, and finally a brief discussion of some 
computer results for haze and clouds. The entire discussion is restricted to unpolar- 
ized light — a scalar theory. Consideration of polarization complicates the discussion 
slightly, and for that reason it has been neglected here. There is no doubt, however, 
that polarization should be accounted for in any multiple-scattering calculation, since, 
in almost any real medium, an initially unpolarized beam of light will be partially polar- 
ized after a single scattering event. 

Quantum Detection Theory 
C. W. Helstrom 

Detection theory evaluates communication receivers on the basis of the average prob- 
ability P g of error and seeks that receiver structure for which P g is minimum. In a 
binary communication system the receiver must decide periodically, on the basis of its 
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input, between two hypotheses H q and H^. When the input can be described by classical 
physics, the optimum receiver bases its decision on the likelihood ratio. When, as in 
optical communications, it is necessary to describe the input data quantum-mechanically, 
the receiver is required to choose which of two density operators p Q and best repre- 
sents the input optical field. It does so by measuring a projection operator, II. If the 
outcome of the measurement is 0, hypothesis H q is chosen; if it is 1, is chosen. The 
best projection operator is the one for which the average probability of error is mini- 
mum; it can be expressed in terms of the eigenvectors of the operator ( 1 — ^>.p j - £p Q , 
where £ is the prior probability of hypothesis H q . The detection of a coherent signal 
in the absence of any noise was discussed; it can be represented as a choice between 
two coherent states. The average error probability so attained is lower than that for 
either counting the number of photons or for using the "classical" detection operator, 
which corresponds to matched filtering. The detection of a coherent signal with ran- 
dom phase occupying many modes of the receiver can also be treated; the error prob- 
ability is expressed in terms of the Laguerre distribution. In instances when the 
optimum detection operator is difficult to calculate, a threshold detector, which is its 
limiting form at low signal-to-noise ratio, can often be specified. 
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Optical-Communication Space Experiments 
J. L. Randall 

The following subjects were discussed: (i) an Optical Guidance Sensor system for 
automatic (unmanned) rendez-vous and docking of two spacecraft, (ii) a precision pointing 
and tracking system for laser communication, (iii) a 30 megabit/sec optical communica- 
tion system, (iv) a 10. 6-p coherent communication system, and (v) OTES (Optical 
Technology Experiment System) proposed optical- communication experiments. 

The rendez-vous and docking guidance system utilizes a GaAs laser and incoherent 
diode as light sources and an image dissector for tracking to provide range, range rate, 
angle, and angle rate which are necessary to perform the rendez-vous and docking 
maneuver. The system is light-weight and requires low power. The system is being 
built by ITT for NASA. 
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The precision pointing and tracking system for optical communication is a prototype 
of an optical communication system that is being tested to study the problems of space 
communication with emphasis on the pointing output. Thus far, all problems appear 
solvable insofar as pointing is concerned. This system was built for NASA by Perkin- 
Elmer. 

o 

The 30-megabit PCM optical communication system uses a 6328 A He-Ne laser with 
incoherent detection. It uses a binary bit code and operates either on-off or in a dual 
channel capacity where the ones and zeros are marked by two circular states of polariza- 
tion. Error rates per bit have been measured under varying conditions. Results indi- 
cate that the system operates near that which can be theoretically predicted with the 
Poisson statistics of photoelectrons. The system was built for NASA by ITT. 

The 10. 6-p communication system utilizes coherent detection by mixing the signal 
with a 10. 6-p local oscillator. The signal is impressed on the transmitter laser by 
frequency-modulating the laser cavity with a piezoelectric transducer. Two voice chan- 
nels are transmitted in the way. The system operates within about a factor of two of 
that predicted. That is, the measured value and calculated value of S/N B agree within 
a factor of two. The system was built by Honeywell for NASA. 

A brief discussion of some of the space experiments that are being studied on the 
OTES was held. These experiments under consideration have as a goal the solving of 
all of the practical problems associated with operating an optical communication system 
in space. This would include communication through the atmosphere and from ground 
to space. The results of these experiments would be the proof of the feasibility of opti- 
cal communication in space and to provide the data necessary to design communication 
systems for mission applications. 
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Designing "Real" Optical Communication Systems 
T. Kinsel 

A communication theorist can generally derive a block diagram for an optimum com- 
munication system designed to operate with a specific transmission channel. This block 
diagram will specify various operations to be performed on the signal waveform. The 
problem considered here is the embodiment of such a block diagram either with existing 
devices or physically possible devices. For an optical communication system a useful 
way for the experimentalist to attack this problem is to consider the location of the inter- 
faces where information conveyed by electrons is converted to photons or vice versa. 
There 'are two such interfaces, one at the transmitter and one at the receiver. It may 
be asserted that the location of these interfaces is variable, within limits, and their opti- 
mum location in a real system will depend upon the characteristics of the various 
devices. It would seem that the performance of "real" systems would generally be 
improved by separating the interfaces as much as possible. In other words, it is 
desirable to do as much of the signal processing by optical techniques as possible. 

An example of a system that represents an attempt to maximize the separation of 
these interfaces is the time-multiplexed optical PCM system as described by the author. 
At the transmitter terminal the interface is located at the array of optical modulators 
which is near the front end of the transmitter. Following this the processing is done 
optically using active devices. At the receiving terminal the interface is located at the 
array of photodetectors. Before this, the process of recovering the separate channels 
is done optically. 

Let us also consider the diversity- channel receiver suggested by Kennedy and 
Hoversten to operate with a turbulent channel. Here the input to the receiver has a dis- 
tribution of phases and propagation directions. A block diagram of the optimum receiver 
is shown in Fig. B-l . 

One design of a real receiver that performs the operations is shown in Fig. B-2. The 
input beam is mixed with a local oscillator whose output is a combination of plane waves. 
Various pencils of light in the input beam will then be selected and imaged onto an array 
of photodetectors. The photon- electron interface is located at this photodetector array 
which is near the front end of the receiver. Subsequent processing would be done elec- 
trically using either integrated circuit approaches or possibly bulk effects in solids. 

An alternative approach that places the interface at a later stage of the receiver is 
shown in Fig. B-3. This approach is highly speculative and not fully analyzed. It is 
used to suggest possible approaches to the problem and possible areas of device 
research. The mosaic array of photodetectors is replaced by an optically nonlinear 
crystal. The number of diversity channels would be determined by the number of plane 
waves in the local oscillator. The output of the nonlinear crystal at the sum frequency 
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Fig. B-l. Idealized. diversity receiver. 
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Fig. B-2. Discrete realization of the diversity receiver. 
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Fig. B-3. Distributed realization of the diversity receiver 
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would be coupled into a multimode traveling-wave optical amplifier such as has been 
discussed by Geusic and Scovil. The bandpass filter would be supplied by the oscilla- 
tion bandwidth of the gain medium. The output would be imaged onto a single photodetec- 
tor which would provide the squaring process and sum over the different channels. 

There are obviously many problems in actually constructing such a receiver and this 
discussion should serve only as an example of the types of signal processing in the opti- 
cal domain that should be explored. 




Appendix C 


On the Reliability of Quantum-Mechanical Communication Systems 

Jane W. S. Liu 

A model of typical communication systems is shown in Fig. C-l. The input to the 

system in any given signaling interval of duration t is one of M input symbols denoted 

m,, m , . . . m.,. For simplicity, it is assumed that the M input symbols occur with 
* ‘•Vi ^ ^ 

equal probabilities. Thus the information rate of the system is R = — - — nats/sec. Cor- 
respondingly, the output of the system in every signaling interval is also one of M sym- 
bols which is the estimation m of the transmitted input symbol m. The receiver is 
designed so that the probability of error P(e) defined as 

P(e)'= Pr [nW m] (1) 


is minimum. 


When the input symbol is m., the electromagnetic field generated by the signal source 

I t 

is in a state specified by the density operator p^. Therefore, in any given signaling 


TRANSMITTED RECEIVED ESTIMATE OF 



Fig. C-l. Quantum-mechanical model of a communication system. 


interval, the signal field is in one of the M states specified by the density operators , 

t t 1 

p 2’ 1 • ‘ P M' ® ecause the presence of the additive thermal noise field, the total 

received field at the receiver is in one of the M states specified by the density opera- 
tors p j , p^, . . . Pjyj. 

Within the time interval during which the signal is expected to arrive, the receiver 
makes measurements on the received field and decides which one of the M input symbols 
is transmitted on the basis of the outcomes of the measurements. We shall assume that 
measurements made by the receiver are ideal. * An ideal measurement is one in which 
the state of the field immediately after the measurement depends only on the outcome of 
the measurement and not on the initial state of the field. It follows that the probability 
distribution of the outcome of any subsequent measurement on the field does not depend 
on the transmitted symbol. Therefore, any measurement made after an ideal measure- 
ment is useless in estimating the transmitted symbol. 

An ideal measurement can be characterized by a complete set of measurement state 
vectors {j x^)}, where x^ denote K-tuples of parameters that are the possible outcomes 
of the measurement associated with the set (Ix^)}. [To simplify our discussions, we 
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assume that the parameters x^, as well as the eigenspectra of the density operators p^, 
are discrete. It is evident that this assumption imposes no real restriction, since all 
discussions can be generalized by replacing the sums with integrals and probability dis- 
tributions with their corresponding probability density functions.] When the field is in 
the state specified by the density operator p^, the probability of the outcome being x k 
is given by 

PfXjJm^ = P(x k | Pi ) = <x k | Pi |x k ). (2) 

Clearly, measurements of dynamical variables that are represented by commutative 
Hermitian operators are ideal measurements. In this particular case, measurement 
state vectors associated with the measurement are just the simultaneous eigenstates of 
the commutative Hermitian operators; hence, {| x k )} is not only a complete set but also 
an orthogonal one. 

Hereafter, we shall specify the structure and discuss the performance of an optimum 
receiver when the density operators p^ are commutative. As an example, the perfor- 
mance of a system in which signal fields have known orthogonal classical waveforms but 
random phase is derived. The performance of the heterodyne receiver and photon coun- 
ter are also included for comparison. 

Another example of systems in which the density operators are commutative is one 
in which both the signal fields and the noise field are in completely incoherent states. 

A fading dispersive channel can thus be modeled quantum-mechanically. The optimum 
performance of the optimum equal- strength diversity system is derived in Eq. (3). We 
shall not discuss this system here. 

I. STRUCTURE AND PERFORMANCE OF AN OPTIMUM RECEIVER 
FOR EQUIPROBABLE SIGNALS REPRESENTED BY 
COMMUTATIVE DENSITY OPERATORS 

Let us assume that the receiver makes an ideal measurement characterized by the 

set of measurement state vectors { | )} on the received field. Then, the receiver uses 

a randomized strategy to estimate the transmitted symbol as m. with probability p., 

1 JiL 

when the outcome of the measurement is x k _ Clearly, 

Y^M 

j = 2, p jk = i; for a11 -• 

j=i 

It follows that the probability of error defined in Eq. (1) is 
. yAVI y 

p(e) ;1 -i 2 Z pjk< x klpjl x k>- (3 > 

3=1 k 
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In general, it is difficult to find the set {|x k )} and probabilities p^ k that minimize 
the right-hand side of Eq. (3). Here, we shall consider only the case in which the 
are commutative. 

Since the density operators are commutative, they can be written 

"r Z v ii bi><Vil. <- 

i 


where (ly^)} is the set of simultaneous eigenstates of the density operators p^; j = 1, 2, 

. . . M. is the eigenvalue of p^ corresponding to the eigenstate | y^). Substituting this 
expression for p^ in Eq. (3), we obtain 


j=l k i 


(5) 


Let ( | Yq) denote the largest of the M eigenvalues y^, Y 2 p • • • Yjyjq for a given i. 

It is clear that the right-hand side of Eq. (5) is bounded from below by 


p ( e > 541 -M Z ^ji) I <^il x k>< X kK> 

i k 

_ , 1 \ ( max \ 


( 6 ) 


This lower bound of P(e) given by Eq. (6) can be achieved by using a receiver that 

measures dynamical variables of the received field which are represented by the set of 

operators {x} with simultaneous eigenstates | y^ ) . Since the | y^ ) form an orthogonal 

complete set, these operators are Hermitian and commutative. Moreover, it follows 

from Eqs. (2), (4), and (6) that the optimum decision procdure is the maximum- 

likelihood decision rule. That is, p., is set to be one if 

3 K 

P({x k }/Vj) & P({ x k }/m i ) ; i, j = 1, 2, . . . M, 

where (x k ) denotes the set of eigenvalues of the operators in {x} associated with the 
state |y k ). p^ is equal to zero for all i; i 4- j. 


II. PERFORMANCE OF THE OPTIMUM RECEIVER FOR FREQUENCY 
ORTHOGONAL SIGNALS WITH KNOWN AMPLITUDE 
AND RANDOM PHASE 


We shall now discuss the optimum performance of a specific system in which the 
M density operators p^ are commutative. For the sake of simplicity, we shall assume 
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that the signal field is linearly polarized and is not spatially modulated. Therefore, in 
the absence of additive noise, the field over the receiving aperture is a plane wave. Let 
us assume that the signal field propagates along the z-axis. Hence, the relevant com- 
ponent of the electric field operator can be expanded as 


E(r,t) = i 



When the transmitted input symbol is rn., the field present at the receiver in the 
absence of the additive thermal noise field is in a state specified by the density operator 


P j 2 IT 


1 I ? 6(2) ( a k _or ik e ^) l c k>< a ki d2 °: 


( 7 ) 


in the P representation. In Eq. (7), the |a^) are the right eigenstates of the annihila- 
tion operators a^. The are complex quantities known to the receiver, while 4>> which 
corresponds to the absolute phase of the received field, is assumed to be a random vari- 
able evenly distributed over the interval (0, 2 tt). That is, the signal field at the receiver 

is in one of the coherent states II Icr., e^). 

k 

In the presence of the additive thermal-noise field, the total received field is in a 
state specified by the density operators 


p 2ir dcf> P 
J n 2ir J 


II exp 
k 


l a k" ff jk e 


i<j>i 


<n> 


KX' 


d 2 a, 


iT<n) 


(8) 


when the message rm is transmitted. That is, the additive thermal-noise field alone is 
in a state specified by the density operator 

I I 2 


P (n) = n 

k 


I 


exp 


<n> 




2 

d a. 


t ( n ) 


(9) 


where (n) is the average number of thermal-noise photons in each mode of the noise 
field. For simplicity, we assume that over the frequency range of interest, it is a con- 
stant. 

It can be shown that the density operators given by Eq. (8) are commutative when the 

complex numbers <j., satisfy the condition 
I K 


l 


a a . , , = 6 ... a • 
jk ] k J 1 


(10a) 


Hereafter, we shall restrict our attention to the set of M orthogonal signals 
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corrresponding to signals in frequency position modulation systems. In this case, the 

o\, also satisfy the condition 
1 K 


cr cr... = cr., 
]k 3 k l jk 


6 ..; 

J3 


for all k. 


(10b) 


(In the classical limit, signals specified by Eqs. (7) and (10b) have known wave- 
forms that constitute a frequency orthogonal set. The absolute phase of the signal 
field is, however, unknown to the receiver. We assume that it is a random vari- 
able evenly distributed over the interval (0, 2ir). The classical electric field ampli- 
tude of the additive noise field can be considered as a white Gaussian random 
process.) 

As a simple example of the type of signals described above, let us consider the case 
in which the signal field is unmodulated. In this case, the density operator in Eq. (8) 
simplifies to 



( 11 ) 


j.y_ 

That is, when the transmitted input symbol is m., only the j mode of the signal field 
is excited while the other modes remain in their vacuum states. Here, we denote the 
relevant modes of the field as the 1st, 2nd, . . . Mth modes. It is clear that these den- 
sity operators can be written 



where In^, n^. , . n^) denotes the simultaneous eigenstate of the number operators aja^, 
a*a 2> . . . a^a M of the individual modes corresponding to eigenvalues n^, . . . n^, 

respectively. Therefore an optimum receiver for this set of M signals measures simul- 
taneously the variables represented by operators aja^ a*a , . . . a^a M . 

Given that rm is transmitted, the joint conditional probability distribution of the 
observed numbers of photons in modes 1, 2, ... M, respectively, being n^, n , ... n M 
is 
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p <w 


M 


/m.) = 


M 


<n) 


\ n l +n 2 + ' ‘ ' n M 


l+<n)/ \ 1 + < n)y 



7! 


r ui 2 1 

V 


la I 2 

Kji 


” exp 

b 1 

n ) (l + <n>)j 



J 

„ 1 + <n)_ 


Therefore, the receiver estimates that the transmitted symbol is rm if the quantity 


f. 

J 


n. 

J 

l 

y=0 



r | a . | 2 1 

V 

r 1 1 2 1 

a. 

1 

1 J 1 

ex P 

1 J 1 

7 ! 

_<n)(l + <n)) - 



1 + <n)_ 


is maximum over all j. For simplicity, we shall further assume that the magnitudes 
of the complex quantities are such that 


-jl 2 = M 2 


PT. 


(12) 


That is, the average number of photons in the signal field is independent of transmitted 

input symbols. Then, f. is the largest among all j's if and only if n. is the largest. 

1 3 


RECEIVED 

FIELD 


INPUT 

PUPIL 



Fig. C-2. Optimum, receiver for equal -strength unmodulated signals 
with random phase. 


Therefore, an optimum receiver is the system of photon counters shown in Fig. C-2. 
We shall express the bounds to the error probability P(e) in the form 

K j exp[— rCE(R)] « p(e) < K 2 exp[--rCE(R)], 

where E(R) is the system reliability function, and C is the channel capacity. We shall 
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be concerned only with the exponential term in the bounds of P(e), since the coefficients 
of K. and K do not depend exponentially on t. 

It has been found that the channel capacity C is given by the expression 


1 + <n> 

C = p In . (14a) 

< n ) 


For rates R in the range 0 < R < R^,, where 

R c (n)(l + <n» 

C (l+2<n» 2 ’ 


(14b) 


the system reliability function E(R) is bounded from below by 


E(R) =* 


( 1+2 ( n)) In 


1 +<n) 

<n) 


R_ 

c • 


For rates R within the range R c « R « C, the function E(R) is given by 


E(R) 


1 + 2<n)-g-- y/l + 4-g-<n>(l + <n» 


<n) In- 


1 + 1 " ^/l + 4-|-<n)(l + <n» 


<n) 


_R 

C 


1 +<n) 


In 


r -\ 

1 + 2<n> (l + <n»-§- + + 4-^<n> (l+(n» 

\ 5 o r • 


In- 


2§ (l + <n» 2 


<n) 


(14c) 


(14d) 


The general behavior of the exponential factor CE(R) is shown in Fig. C-3 for several 
values of (n) and p/(n) =1.0. Also shown in Fig. C-3 is the exponential factor in 
the expression of P(e) when a heterodyne receiver is used. It can be shown that 
when a heterodyne receiver is used for the reception of the signal specified by Eq. (12), 
the system capacity is given by 



(15a) 


For rates R in the range 0 $ 



E h< R ) 




the reliability function is lower-bounded by 


(15b) 
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For rates R in the range -3- < 7?- < 1, the system reliability function is given by 

4 u h 

E h' R > ■ (■ -yy ) 2 ,i5c > 

It has been shown 3 that for arbitrary narrow-band frequency orthogonal signals given 
by the set of density operators Eq. (8), the optimum system reliability function and chan- 
nel capacity are also given by Eq. (14) when 

1 I CT jkl 2 = b| 2 =P T - 

k 

The receiver shown in Fig. C-2, however, is no longer an optimum receiver. According 
to the result in Section I, an optimum receiver should measure dynamical variables 



Fig. C-3. System’ performance of the direct-detection scheme 
illustrated in Fig. C-2. 


represented by operators bjbj, b^, , . . b^.b M simultaneously, where the operators b. 
are given by 

V X ff jk a k : j = 1 , 2 ,... M . 

k 

Let us again denote the outcome of the measurement of bjbj, . . . b^b M by n^, n 2 , . . . 
n M> Then the receiver estimates that rm is transmitted if n. is the largest among 
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all m. (The can be considered as annihilation operators of a set of new normal modes 
with mode functions. 


V j( r,t) = 


°jk ex P 




where the to. are chosen so that the V.(r, t) satisfy the orthonormality condition.) For 
J 1 

this set of signals, the performance of the heterodyne receiver is again given by Eq. (15) 
and Fig. C-3. But the system performance of the direct-detection scheme shown in 
Fig. C-2 is not known. 
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Some Quantum- Mechanical Aspects of Optical Communications 
Carl W. Helstrom 

I. QUANTUM DETECTION THEORY 

A. Ideal Receiver 

A typical optical communication system will transmit periodically one of M signals, 
which, after passing through the intervening medium, appears at the receiver mixed with 
background radiation. The receiver must decide which of the M signals is present in 
its input during each reception interval of duration T; that is, it must choose every 
T seconds one of M hypotheses H^, H^, . . . , H^. Hypothesis H^ states that the input 

to the receiver consists of the signal and background radiation. 

In determining the fundamental limitation on the performance of a receiver in such 
a communication system, it is expedient to imagine the signals as picked up by an ideal 
receiver consisting of a lossless box or cavity with an aperture facing the transmitter. 
Initially empty, the cavity is exposed to the input field by opening the aperture at the 
beginning of the reception interval. After the signal field can be expected to have 
entirely reached the aperture, the aperture is closed. An observer then proceeds to 
make the most effective observations of the field within the cavity, and on the basis of 
their outcomes he decides which of the signals arrived. The same procedure is carried 
out for each interval during which a signal might appear. By virtue of its independence 
of particular choices of detector and amplifier structures, this ideal model brings out 
the limitations on signal detectability imposed by the stochastic and quantum natures 
of the background radiation and of the signal itself. 

The field in the cavity will be described quantum- mechanically by the density opera- 
tor p k when hypothesis H^, that the k th signal was sent, is true. The observer must 
perform whatever measurements he can on the field in order most efficiently to decide 
among the M density operators pj , p^, . . ., p^.. For greatest generality one admits a 
randomized decision procedure whereby the measurements yield M numbers it ^ 

.... ir M , 

0 < ir k < 1, k = 1, 2, . . ., M 

M 

1 "k" 1 ' (D 

k=l 

of which iT k is the probability with which the observer should choose hypothesis H^ . He 
presumably then uses some chance device to make his decision in accordance with the 
M probabilities tt^, tt^, . . . , -ir^. 


Ill 
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The number can be regarded as an eigenvalue of an operator 11-^. The observer 
measures the M operators II 1 , II 2 , .... 11^, and the outcomes of these measurements 
are the probabilities ir 1> ir 2 , . . ., Because of Eq. (1), the M operators must sum to 

the identity operator 1^ 

nj +n 2 + ... +n M = ( 2 ) 

and it is actually necessary to measure only M - 1 of them. 

From the standpoint of the axioms of quantum mechanics as established, for example, 
by von Neumann, the M operators 11^, necessarily to be measured on one and the same 
system, must commute. Hence they can be measured simultaneously or one after 
another without mutual interference. 

B. Broader Formulation 

A less restrictive specification of the M operators is based on the following rea- 
soning. Suppose that an overcomplete set of states |p) characterizes the system, p 
being a continuous n-dimensional parameter. If these are to constitute a resolution of 
the identity, the completeness relation 

J |p><p| d n p = i (3) 

must hold, where d n p is the volume element of the n-dimensional parameter space and 
the integral is taken over all values of p. An example of such a set is furnished by the 
coherent states of the harmonic oscillator, for which n = 2, and p is conveniently a 
complex number. 

The decision procedure is then based on determination of the parameter p by some 
means, followed by the choice of hypothesis with probability ir^tp), k = 1, 2, . . . , M, 
where the w^Cp) are M functions of p summing to 1, 

M 

E ff k (|3) = 1 - (4) 

k=l 

If hypothesis H. is true and the system has the density operator p, , the probability 
h * K n 

that the determination of p yields an outcome in the element of volume d p about 
the point p is 

P k (P) d n p =<p|p k |p) d n p. (5) 

Because of Eq. (3), the n-dimensional probability density function P k (P) is properly 
normalized, 

*We speak of "determining" p, using the word "measure" only in the conventional 
quantum-mechanical sense. 
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( 6 ) 


where 'Tr' stands for the trace of an operator, and the integral is carried out over 
the entire range of (3 . 

The probability that hypothesis H. will be chosen when hypothesis H, is really true 

1 k 


now is 


PrfejlHj = y ir.(p) p k (p) d n p 


= J ir.((3)<p|p k |p> d n p = Tr (p k n.) (7) 

if the operator IF is defined by 

n j = j |p> TT.(p)<p| d n p. (8) 

Because of Eqs. (4), and (3), the M operators II j, II^, . . . , I1 M so defined sum to the 
identity operator as in Eq. (2). If the states |p) are orthonormal, 

<p[p'> = 6 n (p-p'), (9) 

where 6 is the n-dimensional delta function, and the operators II, commute. The pres- 
n k 

ent formulation, therefore, includes the possibility that the operators Il k commute. 

C. Optimization Problem 

The average probability of correctly identifying the transmitted signal is 


M M 

Q= ^ l k K-{H k |H k } = ^ t k Tr(p k n k ) (10) 

k=l k=l 

where l, k is the prior probability of hypothesis H k> 



k=l 

The best receiver can be defined as the one for which Q is maximum. If the M density 
operators p k commute, the most effective measurements to make on the system are of 
a complete set of commuting observables in terms of whose eigenstates the operators 
P k are simultaneously diagonal. The receiver can then process the outcomes of its mea- 
surements as prescribed by conventional detection theory. 

When the M density operators p k do not commute, the determination of the M oper- 
ators Il k that maximize Q is, for M>2, an unsolved problem. For M = 2, measurement 
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of only a single operator II is necessary, and it has been shown 1 that II is a projection 
operator on the eigenstates of the operator ^ z p z ~ ^ 1 P 1 havin g positive eigenvalues. Even 
for M = 2, however, there are few pairs of noncommuting density operators p^, for 
which the eigenvalue problem can be solved and the optimum detection operator II deter- 
mined. 

One problem that can be solved for M = 2 has to do with the optimum choice between 

two pure quantum-mechanical states |l) and |2), which might, for example, be two 

2 

coherent states of the harmonic oscillator. The density operators are 
p ! = h>ol. P 2 = I 2 ) < 2 !> 
and the optimum operator II is 

n = (pj. 

where |t|j) is the eigenstate of ?1 having a positive eigenvalue. The two eigen- 

states of this operator have the form | q) = x^ | 1) + x^ |2), with x^ and x^ coefficients that 
can be determined algebraically. 

In the choice between a coherent state |a) and the ground state |o) of the harmonic 
oscillator, the probability of error attained by measuring II is smaller than that 
incurred by measuring the number of photons in the system or by measuring the oscil- 
lator coordinate parallel to the complex number a, these last two methods corresponding 

to incoherent and coherent detection, respectively. The analysis and the resulting 

2 

error probabilities have been presented elsewhere. 

D. Choice among M Pure States 

If the receiver can be in one of M pure states, represented by the vectors 1 1), |2), 

. . . , |M), and the observer is to decide which one, it is necessary to consider only oper- 
ators that can be expressed in terms of eigenvectors in the M-dimensional space spanned 
by the M state vectors |k), k= 1, 2, .... M. The states ||3) in Eqs. (3)-(8) can now 
be limited to M discrete vectors |Pj), |p 2 ), • • • |P M ) s P annin g the same M-dimensional 
space. The integrals in Eqs. ( 3 )- (8 ) are replaced by sums over this discrete set. 

The "resolution of the identity" 

M 

X |P k > < Pk I = i (1Z) 

k=l 

now requires the vectors |p^) to be orthonormal k = 1, 2, .... M. Indeed, Eq. (12) 
implies that the scalar product of any two of these vectors is 



APPENDIX D 


115 


Let B be the MX M matrix of the scalar products (P m |P n ) : Eq. (13) is then 



Hence B must be the identity matrix, the solution B = 0 being untenable, and the vectors 
|p k ) must be orthonormal. 

As the |p k ) are now orthonormal, the operators 

M 

V 1 Ki’jWkXPtt <»> 

k=l 


defined as in Eq. (8), commute, for they are simultaneously diagonal in the |p k )- 
representation. The average probability of correct decision is 

M M 

Q= £ £ ^(P k )|<P k |j>| 2 . (15) 

3=1 k=l 

where is the prior probability of hypothesis H^, j = 1, 2, . . ., M. It is apparent that 
the orthonormal set |p, ) can be so oriented in the M- dimensional space that one vector, 
say Jp.), is closest to each state vector |j), and 

Upjlj) I 2 > I ( P k 1 0 1 2 > k *3- 

The average probability Q is then maximum when ir.(p.) = 1, it. (p.) = 0, k + j. Randomi- 

1 J K 1 

zation is unnecessary, and 
M 

Q= £ ^l<P j |3)| 2 . (15) 

3=1 

How to orient the orthonormal set |pp to maximize the average probability Q of cor- 
rect decision is an unsolved problem for M > 2. An analysis showing the difficulties 
involved is presented in Section II. 

E . Afterthoughts 

The procedure outlined in part (B) involved determining — somehow— the n-dimensional 
parameter p and then choosing hypothesis H k with probability -rr^p), k = 1, 2, . . . , M, 
where the ^ k (P) are such functions that the average probability Q of correct decision is 
maximum. This probability was expressed in Eq. (10) in terms of M operators 

n k = J IP) * k <PKPl dV 


(17) 
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The integral is taken over the entire range of (3. The decision procedure is not 
necessarily equivalent to measuring these M operators II II ..,,11 

I Z JVl * 

To see this, consider an operator of similar structure, namely the density opera- 
tor p of a harmonic oscillator in thermal equilibrium at temperature t. This opera- 
tor can be expressed in terms of the coherent states Jp) of the harmonic oscillator by 
the expression 

P = J |P> P(P) <P| d 2 p, 

P(p) = (Nf 1 exp(-|p| 2 /N), (18) 

where 

N=l-e w , w = hv/KT, 

3 

hv being the energy per quantum and K Boltzmann's constant. This way of writing p 
resembles (17), but another way of writing it is in terms of the energy eigenstates |n) 
of the oscillator, 

00 

P = 2 l n > q n< n l 
n=0 


% ■ < N+ » _1 ( raj )” <■’> 

A measurement of the operator p can yield only one of the eigenvalues q .. What 
corresponds to the procedure in part (B), however, is determination first — somehow — 
of the complex number p, and then evaluation of the function P(p) in Eq. (18). This can 
yield any of a continuum of positive values between 0 and (N) \ It is possible, there- 
fore, that once the function ir^tp) has been specified, the operator 11^ may acquire a 
spectrum of eigenvalues that does not coincide with the domain of the function ht^P). Mea- 
surement of 11^ may not be equivalent to first determining p and then calculating ir^(p) . 
Furthermore, if the II ^ do not commute and one tries to measure them simultaneously, 
one cannot be sure that the average probability of correct decision is still given by 
Eq. (10). 

Two possibilities then arise: either the set of states |p) to be used in the procedure 
is given in advance, or it is not. If it is given, the probability densities Pj^(P) are also 
known, and the average probability Q of correct decision can be written in the usual way 
as 

M 

« = i y Ok"" "V 

k=l 


( 20 ) 
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This is maximized in the same manner as in conventional detection theory, and the best 
receiver determines the M likelihoods t,^p^((3), k = 1, 2, . M, and picks up the 
hypothesis corresponding to the largest. 

If the set of states |p) is not given in advance, there must be some way of determining 
it from the operators 11^ that maximize Q of Eq. (10), provided the 11^ can be found, for we 
have seen that it may not be possible to attain the same maximum Q by measuring the 11^ 
themselves if these do not commute. How one can specify a set |p) permitting the deci- 
sion scheme to be realized in terms of the probabilities tt^P) is unclear, and the formu- 
lation of part (B) may be incomplete. 

II. OPTIMUM CHOICE AMONG M PURE STATES 

Continuing with the problem of choosing with maximum probability of correctness 
which of M state vectors |l), \z), ..., |m) describes a quantum- mechanical system, 
we write these states as linear combinations of the M orthonormal vectors |(3.), using 
( 12 ), 


M 


l k > = I IV <p j |k> 


M 


j=l 






( 21 ) 


We wish to find how the orthonormal set |pp should be oriented to maximize the proba- 
bility Q of Eq. (16). For simplicity, we take the states as equally likely, = l/M. 

The probability of choosing hypothesis H^, that the system is in state j j), when it 
really is in state | j), is 


Pr{n.~l|H.} = Ujlppl 2 = \^\ Z . (22) 

Since the vector |k) can be multiplied by an arbitrary phase factor without changing any 
probabilities, we are free to assume that = ( p^ j j) is real, and we do so henceforth. 

The average probability of correct decision when the states |k) are equally probable 
is 


Q = 


_ 1 _ 

M 


M 


3=1 



( 23 ) 


We are to maximize this by varying the under the constraint that the cosines 
(j|k) of the angles between the state vectors |j) and |k) be fixed for all j 
and k: 
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Y jt =<j|k>= £ S*iW P i l|S m> 
i, m 

= I 6 ji 6 ki* 


( 24 ) 


since ( (3. 1 6 \ = 6. 

Nr i ir m / im 

-f- / sfe # \ 

Denote by Y k the row vector ^ lk t , 2k . . . 1, and by Y k the complex conjugate 

column vector of the £,^ k , j = 1, 2, . . . , M. Let J k be a matrix that has a 1 in the k** 1 diag- 
onal element and zeros elsewhere. 


( J, ).. = 6. .6.. . 
~k i] i] jk 


(25) 


Then J k is an M X M matrix with £ kk as its k** 1 diagonal element and zeros else- 

where, and 

M 


D 


"I 


J, Y , Y, J, 
~k ~k ~k~k 


(26) 


k=l 


is a diagonal matrix with diagonal elements £ kk . Furthermore, 

M 

MQ = Tr D = Tr ) J, Y, Y, + J. 

~ l_j ~k ~k ~k ~k 


= Tr 


M 

I 

k=l 


k= 1 


J, Y, Y, . 
~k ~k~k 


(27) 


The matrix T whose elements are jf^ = y. k is given by 


M 

I 

k=l 


Y, Y. + 
~k ~k 


(28) 


by virtue of (24). We want to maximize Tr D for fixed T, and we can do so by intro- 
ducing the matrix A of Lagrange multipliers . and maximizing the function 


kj 


jk 


= Tr (D-A D 

<v ~ ~ 

= X T --(J k -A)X k xJ 


of the elements L . of the vectors Y 
kJ ~J' 


( 29 ) 
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If we make arbitrary increments e, in the vectors Y. , we find, to first order in the 
+ / + * * \ ~k 

V where £k = V € lk e 2k ••• e Mk)> 

F({Y.+€.}) -F({Y.}) = 

^ Tr ^k"-* (~ k ~k + ~k ~ k ) ' 
k 

•f. 

At the maximum this must be zero for arbitrary and e^, since the components 
may be complex. Hence we obtain 

<ik-A>Xk “ 


Xil'i-i 1 ■ »' 

which leads to the equations 

ay, = j, y. =l. g. 

k ~k ~k *kk ~k 


Y, + A = Y. + J, = L . G. + , 
~k~ ~k~k *kk ~k 


(30) 


th. 

where G, is a column vector with 1 in the k xn place, and zeros elsewhere. Taking the 
conjugate transpose of the second of these, we see that A must be Hermitian, A = A . 
Writing out the first equation in (30), we find 


M 

I 


^j^jk^kk* 


6., . 

lk 


If we define the matrix U = || IL j || as the M X M array of elements 
U.. = 

i] i] 13 


(31) 


this equation can be written A U = I, where I is the M X M identity matrix. Hence U = 

~ ~ ~ ~ 

A is Hermitian, and its diagonal elements are 1 . 

The "Hermiticity" of U imposes a condition on the coefficients namely 



( 32 ) 


The matrix 
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Y = 


IV = ( ^X 2 


. Y. 


M 


is related to the known matrix T by (28), or by 


(33) 


r = YY + . (34) 

Because of Eq. (32), Y has only —■ M(M+1) independent elements, and by writing out 
Eq. (34) we obtain y M(M+1) distinct equations for them. Unfortunately, the equations 
are nonlinear in a most complicated way. 

From Eqs. (26), (28), and (30) we find 


A r A 


M 

-l 


AY, Y, + A 
~ ~k ~k 


k=l 


M 

= ZiXkXki=S- 05) 

k= 1 

Hence the unknown matrix A is a Hermitian matrix that diagonalizes T, if such can be 
found. Equivalently, 

T = UDU (36) 

1 2 
which also yields a set of-r-M(M+l) equations for the M diagonal elements of D and 
2 2 J] ~ 

the yMIM-l) elements above the diagonal of the Hermitian matrix U. (U^ = 1, k = 

1, . . . , M.) These equations also are nonlinear. 
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Estimation of Modulation Frequency of a Light Beam 
Carl W. Helstrom 

In a system to measure the velocity of a target, the intensity of a beam of light is 
modulated sinusoidally and directed toward it. The reflected light has all its frequencies 
shifted by the Doppler effect, and the frequency of the intensity modulation is also 
changed, by an amount proportional to the modulation frequency and to the radial velocity 
of the target. The reflected light falls on a photoelectric surface, and the emitted photo- 
electrons excite a tuned circuit. The average rate of emission is sinusoidally modu- 
lated, and by finding the resonant frequency at which the response of the tuned circuit 
is maximum, the modulation frequency of the light and hence the radial component of 
the target velocity can be determined. A device for accomplishing this was proposed 
by Nelson McEvoy. 

The question arises how accurate an estimate of the modulation frequency can be 
attained by such a method, and it can be answered at least roughly by an analysis 
involving the Cramer-Rao inequality of statistical estimation theory. 

The photoelectrons are observed over an interval of duration T, which can be taken 
as the integration time of the tuned circuit. Divide this interval into subintervals of dura- 
tion At that are much shorter than the period of the modulation frequency. Let t. be a 
time within the 1 of these subintervals, and let m be the number of photoelectrons 
emitted in that subinterval. The m can be taken as independent random variables with a 

Poisson distribution; the mean value of n. is 

1 

E(n.) = Mtp At, (1) 

where \(t) is the rate of photoelectron emission and is sinusoidally modulated, 

k(t) = \ o (l + cos’m). (2) 

It is the angular frequency 12 that is to be estimated. 

th. 

The joint probability of observing n. photoelectrons in the i interval, i = 1, 2, . . . , 
is 


n. 

p({n.}) = n (n.l)' i [X(t i )At] 1 exp[-X(t i )At]. 
i 


(3) 


The Cram^r-Rao inequality tells us that the variance of an unbiased estimate 12 of the 
parameter 12 of this distribution is bounded below by 


-1 


Var 12 


3 

312 


In p({n.}) 


(4) 
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We proceed to evaluate the right-hand side of this expression. 
The logarithm of the probability in (3 ) is 


In p({n.}) = ^ {n. ln [X(t. )At] - X(t. )At - ln(n.!)} 
i 


and its derivative is 


ln P<{ n J> = I [ n i K_1 ft “ ft At _* 

i 

(5) 

By virtue of (1), the expected value of this expression is zero, 
value is, therefore, 

and its mean- square 

m[wu ln P<t n i}>] = £ (Var n i> x " 2 (fa) ■ 

i 

(M 

Since 


Var n. = E(n.) = Xft^ At, 


we find 


ln p| W> J ■ I W 1 ^) “ 

i 

[ 

. f T dt / axf 
= J 0 x(t)^ 

(7) 

when we pass to the limit At « l/il. 
From (2), 


d\/dU = -X t sin S2t, 


and 


f T dt /0\\ 2 . r T t 2 sin 2 f2t dt 

^0 X(t) X9 ^ " 0 J 0 1 + cosSn 


r T 2 

= X \ t ( 1 - cos «t) dt 
° 


= |x o T 3 [HO(l/OT)]. 

(8) 
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Under the assumption that to T » 1 (many cycles of the modulation are encompassed), we 
obtain from (4), (7), and (8) the lower bound 

Var £! > 3/X. T 3 (9) 

~ ' o 

The Cramdr-Rao inequality provides only a lower bound, and there is no assurance 
that it can be approached. At high signal-to-noise ratio a maximum-likelihood estimate 
of to will have an error variance given approximately by the right-hand side of Eq. (9). 
Hence for such an estimate 


Var to ~ 3/X. T 3 , X. T » 1 

' o o * 


( 10 ) 


since with "shot noise" such as we are measuring here the mean number X. q T of photo- 
electrons serves as a signal-to-noise ratio. Thus 


Var 6 ~ 



1 

X. T’ 
o 


X. T » 1 
o 


(ID 


has the usual form of the error variance of a frequency estimate, being inversely pro- 
portional to the square of the observation time T and to the signal-to-noise ratio. 

If the frequency to be measured is changing at a rate , the observation time T 
should be short enough that 


to T « \/ Var T = 



3 


or 


T « 


( 3 A 0 “ 2 ) 


1/5 


Otherwise an analysis taking both to and to as unknown parameters to be estimated must 
be carried out. For this the two-dimensional form of the Cramdr-Rao inequality can be 
employed. 
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Comments on the Quantal Theory of Communication 
James P. Gordon 
I. INTRODUCTION 

Quantum-mechanical effects, or simply quantal effects, are important to communi- 
cations using optical- frequency radiation. In communication theory, one must recognize 
that information is carried by fields or particles that obey the quantal laws. The 
rewriting of information theory to take this fact into account poses a number of inter- 
esting problems, only a few of which have been solved thus far. On the other hand, most 
systems envisaged now involve some linear processing of an incident field (amplification, 
addition of locally generated fields, etc.), followed by a detector characterized by photo- 
production of electrons at a mean rate proportional to the short-time average of the opti- 
cal power incident on it. By short time is meant a time short compared with the 
reciprocal of the signal bandwidth or intermediate frequency, if applicable, and yet long 
compared with the reciprocal of the optical carrier frequency. In these cases it appears 
that the quantal nature of the field can be rigorously ignored if the basic photoproduction 
process is assumed to have Poisson statistics. 

In Section II we shall briefly discuss the application of quantal laws to communication 
theory, and in Section III we shall give a simple example of how classical fields and 
Poisson photoproduction statistics give correct quantal results. 

II. QUANTAL DESCRIPTION OF COMMUNICATION 

The communication system that we envisage comprises a modulator-transmitter that 
implants information onto a field (or other quantal objects), a transmission path through 
which the field is sent, and finally a receiver that measures some property or properties 
of the field in order to extract the implanted information. We consider that the field 
strength is much attenuated by the transmission path, and perhaps some noise (such as 
thermal background) is added to it. In order to avoid undue complications, we shall 
consider that only one transverse field mode and polarization is important, so the path 
is analogous to a single-mode waveguide. In some reasonably long time interval T, each 
of the independent Fourier components satisfying periodic boundary conditions over 
time T, and hence separated by 6v = 1/T from its neighbors, behaves like a harmonic 
oscillator. The number of such independent harmonic oscillators is equal to the time- 
bandwidth product, and we may liken the model system to an equivalent one in which the 
information is implanted by suitable preparation of harmonic oscillators, which are then 
sent to the receiver at a rate equal to the optical-frequency bandwidth. Now consider 
an imaginary plane cutting the system at the point where the transmitted power 
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reaches its lowest ebb, namely just at the receiver aperture. Everything to the left 
(transmitter side) of this plane may be associated with the preparation of each harmonic 
oscillator as it arrives at the plane, while everything to the right may be associated with 
its measurement. 

A harmonic oscillator is described classically by a coordinate q and a canonically 
conjugate momentum p. The proportionality between q, p and the electromagnetic 
fields E, H is well documented. The energy W of the oscillator is 

W =y(p 2 +w 2 q 2 ). (1) 


In quantal theory q and p become noncommuting operators satisfying Heisenberg's 
uncertainty relation, and W becomes the Hamiltonian operator. 

A pure state of preparation of an oscillator is described no longer by assigning values 
of q and p, but rather by assigning a complex wave function ^(q) after Schrodinger, or 
a ket vector |t) in Hilbert space after Dirac. The index t is a (possibly multidimen- 
sional) index ranging over all possible states that the transmitter may prepare. The 
index may be thought of as identifying the various possible knob settings of the modula- 
tor. As a pertinent example we have the coherent states championed by Glauber, for 
which the normalized wave functions are given by 


CO 

1/4 

u(q-q) 2 ipq] 

ttK 

exp" 

zn + R 



L J 


+5, pW = 


Following Klauder, we can represent these states by kets labeled 

I q» p) • 


(2a) 


(2b) 


These are the states which would be prepared by a highly coherent laser followed by a 
modulator that could adjust amplitude and phase, followed by a noise-free (0°K) path. 
Another example might be the energy states, whose wave functions are the Hermite- 
Gaussian functions, and which we can represent by the kets 

| n) . (3) 

Here the index n represents the number of quanta of energy hv in the state. Because 
of the attenuating path, these pure states cannot be prepared in our system; and this is 
unfortunate because they would give us a communication channel with the maximum pos- 
sible information rate. An exception is the ground state |o), which is identical with the 
coherent state q = p = 0. 

A rather general description of the measurement process involves two stages: first, 
a reversible coupling of the original quantal system to an "auxiliary" quantal system 
whose initial state is set by the receiver, followed after a suitable interval by the irre- 
versible exact "observation" of a complete set of commuting observables for the expanded 
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(original plus auxiliary) system. This description allows for minimum uncertainty com- 
bined measurements of noncommuting observables of the original system. For example, 
if we wish to measure both position and momentum of a harmonic oscillator, we might 
couple half its energy out into a second identical harmonic oscillator which was initially 
in its ground state, whereupon exact "observation" of the momenta of the two separate 
oscillators at suitable later times would give us measures of both the position and 
momentum of the original oscillator. It turns out that the results of all such measure- 
ments if done as well as quantal theory allows can be expressed in the form 



4* r (q) 4* t (q) dq 


l<r|t)| 2 . 


(4) 


Here r is an index that covers all possible measurement results (meter readings), and 
P(r/t) is the conditional probability that the measurement will give a reading r (again 
possibly multidimensional), given that the transmitter prepared the systems in the state 
|t). The state |r), or 4< (q), represents a state of the original system (harmonic oscil- 
lator) which is associated through the measurement process with the result r. If the 
set of prepared states |t) are unit-normalized so that 




|t> = J 


4* t (q) + t (q) dq = 


and we write 

P(r/t)= (t|r) (r|t), 

then we see that p(r/t) is a probability function over r if for any t we have 

2 P(r/t)= < 1 1 [ 2 |r)<r|]|t> = 1. 
r r 


Thus we require only that the set of measurement states satisfy a completeness relation 

2 |r> ( r | = I, (5) 

r 

where I is the unit operator. The notation |r) (r| represents the outer (tensor) product 
of the ket |r) with its corresponding bra (r | in the conjugate space. An important point 
here is that the set of measurement states |r) need not be orthogonal. For pertinent 
examples, the measurement states corresponding to energy measurement are the energy 
states |n), which satisfy the completeness relation 

2 |n) ( n | = I. 
n 

Correspondingly, for energy measurements, we find 
P( n /t) = |( n 1 1) | 2 . 
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Likewise the position eigenkets |q) satisfy the completeness relation 

J I q) <q| d q = i. 

and correspondingly 

P(q/t) = |< q 1 1) | 2 = |4/ t (q)| 2 . 


Here, of course, the probability function p(q/t) is continuous rather than discrete as was 
the case for the energy measurement, and we achieve the standard result that the abso- 
lute square of the wave function is the probability function for position measurement. 

But now our formulation covers minimum uncertainty measurements of both position and 
momentum. The coherent states also satisfy a completeness relation 

J I ' V’ p m> < V’ P m I dq m dp m = < 6 > 

where the individual vectors are unit-normalized, and q m P m replace q, ip in (2) and 
represent the measurement results. Correspondingly, a minimum uncertainty measure- 
ment of both p and q can result in a conditional probability function 

P(q m’ p m /pq) = zk l< p m’ p ml pp >l 2 


= 2 ^f exp 


“ ( q m ~ q)2 _ (p m~ p)2 \ 

2fi 2Kio J 


(7) 


where in the last line we have worked out the answer, using the wave functions of (2a) 
and the indicated integration of (4). Note that the variance of the q part of the measure- 
ment, namely R/oo, is twice as large as we could have achieved by an exact measurement 
of q, for which the result could be 


p ( q m / q > P) = 


l‘ |j q,p (q m ) l = 


, 1/2 


L ttK_ 


exp 


u(q m -q ) 2 


( 8 ) 


where we have used (2a). Thus we cannot obtain information about p without giving up 
some ground in the measurement of q, which is also in accord with the uncertainty prin- 
ciple. We may note that the uncertainty principle has two faces, one saying that we can- 
not prepare a system with exactly defined p and q, the other saying that in addition we 
cannot measure both p and q with arbitrary precision. Note the similarity between (7) 
and (8) and the quantal noise in heterodyne and homodyne receivers. In the former case, 
both in phase and quadrature components of the field can be used to carry information, 
and we find a minimum effective input noise energy per mode of hv, just as in an ideal 
laser amplifier. 
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In the latter, while we give up all notion of using the quadrature field component, we find 
the minimum effective quantal input noise reduced by exactly a factor of two, to hv/2. 

With one further step, we shall have completed the formalities. Equation (4) gives 
the conditional probability result for the case in which the transmitter can prepare pure 
states and the receiver can make what we may call ideal measurements, namely the best 
measurements of a particular type that quantal theory will allow. But in real systems 
neither of these is the case. For example, the transmitter may attempt to prepare a 
coherent state I'p, q), but what appears is a state |p+P n > q+q n ), where p n and q n repre- 
sent the random-noise field added in the path. Given the statistical probability distribu- 
tion P(P n » q n ), Eq. (4) must now be modified to read 

P(r/q, p) = S P(p n , q n ) | < r | P+P n> q+q n > 1 2 (9 ) 

p n’% 

which may be rearranged to read 

P(r/q, p) = (r|p(p,q)|r), (10) 

where 

p(p,q)= s P(P n . q n ) |p+p n . q+q n ) (p+p n . q+q n l 

p n’ q n 

is the statistical density operator representing the noisily prepared system. For thermal 
noise the distribution P(p , q ) is Gaussian, just as in the classical theory, and, of 
course, then the double sum indicated in (9) and (10) becomes a double integral. Finally, 
the receiver may not be ideal. For example, its meters may be perturbed by unessen- 
tial noise, or perhaps the observer simply does not care to distinguish each possible 
value or r. He might then obtain a more coarse-grained conditional probability distri- 
bution such as 


P(Vq.P) = 2 < r I p(p> q) I r ) (11) 

rei 

where the set of A's represent a complete set of mutually exclusive ranges of r. 
Equation (11) can be re-expressed in the form 

P(A/q,p) = Trace [p(p, q) cr( A)], (12) 

where 

cr(A) = Z | r) ( r | . 
re A 

Equation (12) is the generalization of (4) to include noisy preparation and measurement 
of the quantal systems. Without ambiguity, we can rewrite (12) in the general 
notation of (4), namely 
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P(r/t)= Trace [p(t)cr(r)]. (13) 

To recapitulate, p(t) is a density operator for the information-carrying quantal system 
normalized to unit trace, which is associated with the particular transmitter setting t, 
and includes the effects of the lossy and noisy transmission path. And cr(r) is a density 
operator for the information- carrying quantal system, but is associated with the particu- 
lar receiver reading r. The set cr(r) must be complete, in the sense that 

2<r(r) = I, (14) 

r 

so that the conditional probability sums to unity over r. The possible set of measure- 
ment operators, however, need not be restricted to orthogonal sets. 

If (13) encompasses all possible measurements of quantal observables, a backwards 
question may be asked; namely, Does every complete set of density operators o-(r) corre- 
spond to some measurement? The answer is not obvious, and deserves some thought. 

A simple example may serve to illustrate the problem. Helstrom has considered a 
binary system in which the transmitter prepares each oscillator either in the ground 
state 1 0, 0 ) or in a coherent state |"p, q) with some a priori probability X Q for the 
ground state, and = 1 - X. Q for the excited state. For this problem the two states by 
themselves form a complete set. Helstrom has shown that a certain pair of orthogonal 
measurement states, which are the eigenstates of the Hermitian operator 

\ 0 | 0 , 0 > < 0, 0 | -Xjp.q) (p,q|. (15) 

minimize the average error probability for the system. These measurement states have 
the form 

|r) = a r |0, 0} + b r |p, q), r=l, 2. (16) 

The results above can be shown solely on the basis of Eq. (4), without any reference to 
how one might make such a measurement, even in a gedanken world. In this case, how- 
ever, the measurement turns out to be possible in principle if one admits the possibility 
only of energy measurements . Without loss of generality we can set p = 0. We must 
wait only a suitable fraction of a cycle for this to be true, since p and q obey the clas- 
sical harmonic equations of motion, while the oscillator is left undisturbed. When "p = 0, 

then a and b turn out to be real, and hence 
r r 

4 J r (q) = ( q I r > = a r + 0 j 0 (q) + b r Vq(q) (17) 

is a real wave function in which there is a weighted sum of two coherent- state wave func- 
tions. Thinking of our oscillator as a mechanical oscillator, we can alter suddenly the 
potential energy function so that one of the two measurement states (16) becomes 
an energy eigenstate of the new Hamiltonian. Energy eigenstates are found from 
Schrodinger 's equation (for unit mass) 
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R z a 2 qj 

^ y + (W-V(q)H = o. (18) 

2 aq 


Hence, given 4< r (q)> we have for the new potential 


V(q) = W + 


r 2 a2 Vq ) 

2+ r (q) aq 2 


It does not matter whether the other measurement state is an energy state for the new 
potential or not. We calculate all of the energy states for this new potential, and then 
having made the sudden change, make an energy measurement at leisure. Either the 
oscillator was in the singled-out ^(q) state or not, and we find the answer we want. This 
idea of making measurements via a sudden alteration of the potential function in the 
Hamiltonian of the system involved is due to W. E. Lamb. 

Based on (13) as the fundamental premise of quantal communication, a number of 
general things can be said. First, an upper bound to the information rate of the channel 
per quantal system (per harmonic oscillator here) exists which is based solely on the 
set of prepared density operators p(t) and the a priori probability \(t) that each is used. 
In quantal theory, the entropy function S(p) is given by 


S(p) = -Trace [plogp] 

= - S p k log p k , 
k 

where the p k are the eigenvalues of p. The upper bound is given by 

B = S(p) - 2 k(t) S( P (t)), 
t 


where 


p = 2 X(t) p(t). 
t 

The properties of the transmitter determine the set p(t). Thus we can optimize the 
transmitter in a rough way by choosing a set of p(t) and corresponding \(t) so as 
to maximize B. Given any particular variety of receiver, the set of measurement 
operators <r( r) can be determined, and then the information capacity of that channel 
can be determined from the set of conditional probabilities P(r/t) and the \(t), and then 
we can see how good the channel is by comparison with B. 

Helstrom has considered M-nary channels in which the set p(t) is a set of 
M pure quantal states 
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In this case any set of M measurement states |r) always forms an orthogonal set, 
and in the case of the binary channel he has been able to find the pair that minimizes 
the error probability. We have discussed how the binary measurement might in prin- 
ciple be accomplished, but the question of what the optimum set of measurement states 
is for the M-nary channel, and whether the Hermitian operator that generates them is 
in any sense an "observable, 11 remain unanswered. Also, the question of how additive 
noise affects Helstrom's results should be looked into. 

III. BASIS FOR CLASSICAL FIELDS AND POISSON STATISTICS 

To return to a more practical level, we may consider the class of prepared density 
operators for a harmonic oscillator which can be described in the so-called P repre- 
sentation of Glauber, namely 

p = jj p(p, q) | p. q> < p. q)| dpdq, 


with a real, non-negative probability function P(p, q). The probability distribution for 
the number of quanta in such an oscillator is given by 


<n|p|n) = JJ P(p,q)S_ e m 


dpdq. 


where 


m = 


1 

2Ru 


, 2 , 22 . 

(p +co q ) 


is the average energy of the state |p, q) expressed in photon units. Thus the photon 
number distribution is simply a Poisson convolution of the "classical" distribution 
P(p, q). An oscillator weakly coupled to a heat bath and in equilibrium with it is sta- 
tistically described by a Gaussian P function, 


P(P> q) « exp ^ 


2 2 2 ] 

p +io q I 

2nfico 


where 



The field emitted by an ideal oscillator is very nearly described by a P function that 
is a 6 function. It is just this P distribution that we can consider as describing "the" 
classical field distribution. It obeys the classical laws of motion of the harmonic 
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oscillator throughout all of the linear operations that are performed at the receiver. The 
Poisson convolution makes its appearance only at the detecting element, where in essence 
we ask the question, How many photons are there ? While there are many conceivable 
fields that are not described by such non-negative P functions, it is not likely that they 
will ever be encountered in any practical problem in communication. 
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Equivalence of Transmitte r and Receiver Antenna Gain 
in a Random Medium 

David L. Fried 

We are concerned here with the relationship between the performance of a ground- 
based optical heterodyne receiver and a ground-based laser transmitter operating 
through the atmosphere in a ground-space link. In particular, we are interested in the 
effect of the random refractive index on the system performance. We shall show that 
the atmosphere has an identical influence on both systems. 

For the purpose of analysis, rather than consider the actual configuration of an opti- 
cal heterodyne receiver we shall consider a much simpler system composed of a 
diffraction-limited lens and a pinhole, at a distance one focal length away from the lens. 
The pinhole diameter is small enough compared with the Airy disk, which the lens might 
form, that the pinhole diameter may be considered to be effectively of zero extent. In 
back of the pinhole a photodetector is placed to measure the optical power coming through 
the pinhole. The significance of this very simple configuration is that it corresponds 
exactly to the heterodyne receiver. It is easy to show, either formally or heuristically, 
that the wavefront distortion effect of atmospheric refractive index inhomogeneities is 
identical in reducing heterodyne signal power and in reducing the optical power that can 
pass through the pinhole. The simplest proof of this follows from recognition that (i) the 
optical power through the pinhole measures the Strehl definition of the composite lens- 
atmosphere optical system, (ii) the Strehl definition is equal to the integrated MTF of 

the composite system, 1 and (iii) the expression for the integrated MTF as given else- 
2 

where is identical in its form to the expression for the heterodyne receiver relative to 

3 

signal-to-noise ratio as previously given. Although these expressions, as given pre- 
2 3 

viously, ’ apply to values averaged over an ensemble of atmospheric distortions, equiv- 
alent results can be obtained for the instantaneous, or unaveraged performance. 

In place of the laser transmitter, we consider the same lens/pinhole arrangement 
but with a light source placed in back of the pinhole. The equivalence to a diffraction- 
limited laser transmitter is obvious. In this case we seek to determine how atmospheric 
turbulence affects the on-axis far-field intensity. 

For the purpose of analysis we define the three coordinate systems Xj, x 2> and 0, 
which are (for x^ ) linear position in the plane of the pinhole, denoted Xj , (for x 2 ) linear 
position in the plane Xg, a plane just above the atmosphere, and (for 0) angular position 
in the far field away from the atmosphere. (0 is a vector in the direction space just as 
x is a vector in the ordinary space X.) The pinhole and lens define a center line for the 
system. Xj = 0 (or x 2 ) occurs where the center line passes through Xj (or X 2 ), and 0 = 

0 is the direction vector parallel to the center line. 

Our analysis is based on the well-known fact that if Gfx^ x 2 ) is the Green's function 
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that propagates a wave function u^Xj) or into the wave function u^fx^) or X 2 ; 
that is 


u 2 (x 2 ) 


I 


dx. 


G(x 1 ,x 2 ) u x (Xj ) 


X, 


( 1 ) 


then G (x^.x.,) propagates a wave function from X 2 to X^; that is. 


u l (x l 


- I 


dx 2 G (Xj, x 2 ) u 2 (x 2 ). 


X_ 


( 2 ) 


The Green's function G reflects all of the boundary conditions associated with the region 
between Xj and X 2 , including both the presence of the lens and the irregular refractive- 
index pattern of the atmosphere. 

The wave function u^(0 ) associated with an infinite plane wave beyond X 2 traveling 
in the direction 0 may be computed from Eq. (3), 


u 3 (0 ) 


■S 


dx„ exp(ik0 


x 2 ) u 2 (x 2 ). 


X, 


(3) 


where we have suppressed a constant of proportionality as being irrelevant to the details 
of our calculation. The inverse relationship permitting us to calculate u 2 (x 2 ) for a plane 
wave traveling from the direction 0 2 , incident on X 2 is 

u 2 (x 2 ) = J d0 exp(-ik0 ■ x 2 ) u 3 (0). (4) 

© 

To analyze the problem of transmitter performance, we assume that Uj(Xj) is a 
delta function centered at = 0; that is, 

UjtXj) = 6(Xj) (5a) 


u 3 (0) = (0-0 o ), 


(5b) 


. _ jj( . . 

and solve for the power density in the far-field direction 0 Q , that is, for u 3 (0 q ) u 3 (0 o ). 
This is a measure of the efficiency of a ground-based transmitter in directing optical 
power to a receiver in space, an angular distance 0 q away from the optical center line. 
This is directly proportional to the transmitter antenna gain, G T> Combining Eqs. (5a), 
(1), and (3), we get 


S 

X„ 


dx 2 exp(ik0 Q 


x 2 ) G(0 ,x 2 ), 


u_(0 ) = 
3' o' 


(6a) 
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so that G,p is proportional to 


u~ (0 ) u, (0 ) 
3' o' 3' o' 


J dx 2 dx^, exp[ik0 o • (x. 


2 2 


)] G(0,x z ) 


G*(0, xi,). 


X„ X- 


(7a) 


To analyze the problem of receiver antenna gain, we assume that u 2 (0) is a delta 
function centered at 0 , corresponding to a transmitter in space an angular distance 0 
away from the ground-based receiver center line, and ask about the power density at the 
pinhole; that is, we solve for UjfO) UjfO). This quantity determines the received signal 
power and is proportional to the heterodyne receiver antenna gain, G^, . By combining 
Eqs. (5b), (2), and (4), we get 


Uj(0) 


■f 


dx^ exp(-ik0 Q 


X 2> 


(0,x 2 ). 


X„ 


(6b) 


Hence the receiver antenna gain, Gj^, is proportional to 


u*(0) u, 


( 0 ) 


■n 

x„ X* 


dx 2 dx' 2 exp[ik0 Q 


(x. 


-x' z )] G(0,x 2 ) 


G*(0,xy. 


(7b) 


Considering the fact that all of the atmospheric refractive -index irregularities are con- 
tained in the Green' s function G, and noting the exact equivalence of the right-hand sides 
of Eqs. (7a) and (7b), we can conclude that a given optical system, whether used as a 
diffraction-limited transmitter or a heterodyne receiver will be equally affected by the 
atmosphere in terms of instantaneous and average optical signal power transfer. This 
means, for example, that if a large transmitter and heterodyne receiver share common 
optics, then if, because of a particular atmospheric refractive index configuration, the 
heterodyne signal increases in value (or fades) then the fraction of transmitter optical 
power reaching the space platform will increase in value (or fade). Such a situation 
offers, for instance, the possibility of burst communication. It also means that if the 
optical axis is rapidly adjusted to maximize the heterodyne signal strength, thereby 
tracking out some of the atmospherically induced wavefront distortion, this pointing will 
automatically maximize transmitter performance. 


A particularly useful consequency of the result that we have obtained is that we can 

3 

use the results of the previous paper, which describes optical heterodyne receiver per- 

2 

formance with atmospheric effects considered, or even of the other paper, which, ana- 
lyzes the performance of an imaging system to predict the average effect of the 
atmosphere on a laser transmitter. We could, for example, conclude from this that a 

laser transmitter would not be limited by atmospheric effects if its diameter were ~r /z, 

3 ° 

where r Q is a coherence length previously defined. For a transmitter of diameter r Q 
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the average antenna gain will be down approximately a factor of 2 from the ideal 

achievable if there were not atmospheric refractive -index inhomogeneities. 

The reciprocity between transmitter and receiver that we have demonstrated also 

4 

permits us to use other results to calculate the variance of transmitter antenna gain 
as affected by the changing pattern of atmospheric inhomogeneities. 

It should be noted that although this analysis has been couched in terms of a ground- 
space link, it should be possible to modify the proof so that it will apply equally well for 
the two ends of the link both in the atmosphere and on the ground. 
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Atmospheric Turbulence Considerations for a 
Ground-Based Laser Transmitter in a Ground-Space Link 

D. L. Fried 

The kinds of questions that we wish to be able to answer about atmospheric turbulence 
effects on a ground-based laser transmitter concern (i) average antenna gain, (ii) the 
variance and other statistical aspects of antenna gain, and (iii) to what extent can antenna 
gain be adjusted by matching the transmitted wavefront to the distorted wavefront of a 
beacon arriving from the receiver position, and what is the speed with which these 
effects fluctuate. 

There are two ways in which we can approach the problem. The first method involves 
direct calculation of laser beam statistics starting from Schmeltzer's 1 results or per- 
forming an equivalent more extensive treatment. It is then necessary to carry out numer- 
ical evaluations from this. Published results are available for laser beam intensity 
2 3 

variance ’ for both ground-to-ground and ground-to-space links. The published results 

are available in the form of log-amplitude variance data, C.(0). To obtain intensity vari- 

2 4 * 

ance data v ^ , we use the assumption that the fluctuations are log-normally distributed 

so that 


o-j. = exp[4 C £ (0)] - 1. 

Predicted values for transmitter antenna gain, G a ^. mos » relative to the ideal gain, 
Ordeal’ obtainable in the absence of atmospheric turbulence can be determined from the 
log-amplitude mean, £, and variance, C^{0) (again under the assumption of a log-normal 
distribution), from the formula 


G 


G = 


atmos 

^ideal 


exp{2[C.(0)+£]}, 


where gains are directly proportional to power (rather than in the dB form, propor- 
tional to the logarithm of power). The value of £ has recently been calculated, from 
Schmeltzer,* by F. Gebhart (O. S. U. Electro Sciences Laboratory) and the results 
should soon be published. The calculations, however, only apply to a horizontal path 
and so are not directly applicable to our ground-space link problem. 

Because a Gaussian tapered apertured is assumed in all of the published work on laser 
beam statistics, and it appears that a uniformly illuminated finite-diameter source may 

have significantly different scintillation statistics, we do not believe the results of 
3 5 

Fried are directly applicable to our problem. As a consequence of this and the lack 
of calculated data for £ on the ground-to-space path^ there is little that we can determine 
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with certainty from a direct study of laser beam propagation statistics. Consequently, 
most of what we can say about transmitter performance must be deduced from our sec- 
ond approach to the problem. This approach is based on the exact reciprocity, as dem- 
onstrated in the Appendix, between optics acting as a transmitter or as a receiver over 
a given randomly inhomogeneous path. This very powerful result permits us to consider 
the receiver in space as replaced by a laser source, and the transmitter on the ground, 
whose performance we wish to study, as replaced by a heterodyne receiver, and equate 
instantaneous, as well as average atmospheric effects, on the heterodyne receiver 

antenna gain with what the ground-based optics could achieve if it were functioning as 

7 8 

a transmitter. This immediately permits us to apply the results of Fried ’ to the 
evaluation of the transmitter problem. Moreover, since it can be shown that imagery 

Q 

resolution as defined by Fried is equivalent to antenna gain as applied to heterodyne 

o 

reception, we can also use these results to study transmitter performance. Fried's 
7-9 

results together permit us to talk about transmitter antenna gain, gain variance, and 
the effects of rapid angle tracking on gain, for a uniformly illuminated transmitter aper- 
ture operating on a ground -to -space link. 

Using this approach, we draw the following conclusions regarding the performance 

of a continuous (as distinct from burst operation) static (not angle-tracking) transmitter. 

We restrict attention here and throughout the rest of this appendix to wavelengths of 

approximately 0.5 p and 10.6 p. The maximum average achievable antenna gain is that 

associated with a diameter r . For X = 0.5 p, r is of the order of a few inches. Between 

o o 

one and four inches may be taken as a good nominal range of values, although the exact 
value will depend on local atmospheric conditions, the zenith angle of the path, and pos- 
sibly nearby geographic features. For X = 10. 6 p, we expect r to be larger by a factor 

f /r O 

of (10. 6/0. 5) ' , with a nominal range from just under 1 m to possibly as much as 3 m. 
(It is interesting to note that since ideal antenna gain is proportional to (x/D) , the lim- 
iting antenna gain at 10.6 p will be better than that at 0.5 p by a factor of (10.6/.5) 2 ^ ~ 3 ) .) 
The variance of antenna gain will show, however, a strong dependence on the ratio of 
diameter to r . For D = r Q /2, the rms fluctuation will be approximately 30%, but for 
D = r Q it will be approximately 70% . In view of the fact that we expect the fluctuations 
to be log-normal, the 70% rms fluctuation is probably intolerable, and even the 30% 
represents a problem, although probably a tolerable one. In forming these estimates 
for scintillation, we have assumed that there was no gross intensity fluctuation of 
the kind that in the reciprocal problem would make the total power entering the aperture 
vary. A consideration of stellar scintillation data and our estimate of its wavelength 
dependence (that is, C^(0) oc X. -7 /^) indicates that for X = 0, 5 p we may have an addi- 
tional 30% rms fluctuation to convolve with the effect discussed above (giving a total of 
approximately 40% for D = r /2), but for X = 10. 6 p there should be little additional 
effect. We conclude that for a static laser transmitter to minimize fluctuations of gain, 
the optics diameter should be ~r Q /2, in which case we shall have 30-40% rms fluctua- 
tion. It may be desirable to build a transmitter with a variable diameter. 



APPENDIX H 


141 


If a transmitter is designed to operate on a burst communications basis, then it would 
be desirable to have a diameter much larger than t q and maximize the antenna gain fluc- 
tuations. Transmission would then be initiated by monitoring reception from a beacon 
on the space receiver. Whenever the beacon signal strength rose above a critical value 
a burst of transmission would be initiated. It is difficult to estimate the mean time 
between the occurrence of burst conditions, but it is reasonable to estimate that it is of 
the order of D/v or r Q /v when v is the velocity with which the distorting medium moves 
through the optical beam. We may estimate v as being comparable to the wind speed. 

It is questionable whether, at the low rate at which bursts occur, adequate message 
memory is available for this kind of operation. 

It is possible that if the transmitter pointing is slaved to a receiver that tracks a 
beacon on the space station, the transmitter pointing will track out that part of the atmo- 

Q 

spheric disturbance which appears as wavefront tilt. Using the results of Fried and 
noting that we expect very much less than 1 Np rms log-amplitude fluctuation in all cases 
(so that we use what are referred to as the near-field results^), we may expect average 
performance with a diameter D » 1,75 r Q as close to the diffraction limit as the r o /2 
diameter provides in the static case, that is, a gain down ~30% from the diffraction limit 
for that diameter. We believe, but this is only conjecture, that this means a tracking 
system with D = 1.75 r will have the same rms fluctuation of gain that a static system 
with D = r Q /2 has, namely 30%. Experimental data for the spectrum of image motion 
indicates that most of the motion is below 10 Hz, and all of it below 100 Hz in the vis- 
ible range. We believe the frequency will be lower by a factor of approximately 20 for 
1 0-p observation. 

Perhaps the most serious questions to be raised concerning beacon angle- 
tracking transmitters involves the question of "point ahead." To reach the beacon, 
the transmitted signal must be pointed ahead of the beacon an angle 0 = 2V/e, 
where V is the relative velocity perpendicular to the line of sight between the 
transmitter and the receiver, and c is the speed of light. For a nominal veloc- 
ity V = 10^ miles/hour, 0 = 3 X 10 ^ rad. We may take the scale height of 
the atmosphere as H = 6300 m, the distance in which the density of the atmo- 
sphere falls to e * of its sea-level value. We argue that most of the wave- 
front distortion is produced where the atmosphere is most dense, and if the 
received and transmitted beams are nearly congruent over that range, that they 
will see the same atmospheric effects. Hence H0 » 20 cm. This is larger than 

1.75 r from most cases when X = 0. 5 p, and is much smaller than 1.75 r for 
o r o 

X = 10.6 p. Thus it appears likely that point-ahead problems will prevent any 
significant improvement in transmitter antenna gain by beacon angle -tracking in 
the visible range, but will definitely be usable at 10 p. We may comfortably 
expect to be able to use a 1. 5-m or 2. 0-m diameter, beacon-tracking, 10-p 
laser transmitter at 70% of diffraction-limited antenna gain, with an estimated 
30% rms fluctuation of gain. 
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Appendix I 


Spatially Distributed Coherent Heterodyne Reception through a 
Fluctuating Atmosphere 

P. E. Green 

In order to develop maximum SNR at the receiver output, it is most desirable to 
sample more of the received signal wavefront than simply that over which phase and 
amplitude remain more or less constant. The means to be described are available in 
principle to extend the usual notion of coherent heterodyne optical reception so as to 
embody the optimum use of energy over a large area of wavefront. 

We shall use as our starting point receivers (discussed elsewhere*) that receive a 
signal at each of D spatially separate detectors, pass each detector output through 
appropriate matched filter-envelope detector elements (or the equivalent correlators), 
and then, after a nonlinear weighting of each signal, combine outputs, as shown in 
Fig. I- lb. (For convenience, we shall assume throughout that binary frequency-shift- 
keyed signals are transmitted, as in Fig. I- la.) Each input element is a detector 
whose aperture is a portion of the received wavefront roughly equal to the area over 
which phase is constant. The nonlinear operations serve the function of emphasizing in 
the output the contribution of those elements whose signal strength is greatest while 
de- emphasizing the weaker ones. 

We are motivated here by the notion that a considerable improvement in performance 
against noise may be expected from a receiver which, instead of using nonlinear ele- 
ments to do the weighting as in Fig. I— 1 , multiplies each channel before combining by 
the long-time average strength of that channel. To use all of the available data in the 
received wavefront, it is desirable not only to do such a weighting of individual spatial 
regions by the most noiseless estimates possible of their current amplitudes but also 

to estimate the (slowly varying) received phase at each element and compensate for 
* 

this. 

We discuss two procedures for accomplishing the phase compensation and amplitude- 
weighting functions. First, we describe a receiver structure that is discrete in space; 
that is, it is a mosaic of elements, as is the structure of Fig. 1-1. The second struc- 
ture is identical in function to the first but treats the spatial coordinate as a con- 
tinuum. It appears on cursory examination that the spatially discrete version could 
be built now with state-of-the-art components, whereas the spatially continuous 
receiver, although visualizable in detail, requires components that do not exist at 
present. 


Devices that combine signals in such a way may be viewed as "optimum multipath 
receivers," "adaptive arrays," "optimum diversity combiners," "phase-front correcting 
antennas," etc. , depending on one's point of view. 
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(a) 



ffo) <2<*> 

REFERENCE COPIES OF 
TRANSMITTED SIGNALS, 
BUT FREQUENCY-SHIFTED 
BY Uj + Uj 


FILTERS OF BANDWIDTH 
ROUGHLY 


(b) 


Fig. 1-1. (a) Transmitter, (b) Optimum receiver when averaging of slow 

channel fluctuations is ignored ("min-max phase" receiver - 

Kennedy and Hoversten 1 ). Bandpass correlators are shown 
instead of matched filters. 


A. .Spatially Discrete Version 

A receiver 2 that weights each of the D diversity element outputs by a smoothed esti- 
mate of its current amplitude, and at the same time corrects phase differences from 
element to element, can be realized by replacing the operations within the dashed line 
of Fig. 1-1 with the operations shown in Fig. 1-2. Figure 1-2 is labeled to show the 
amplitude, phase, and frequency at various points; the bar indicates a time averaging 
that is assumed to take place by virtue of the narrowness of the "very narrow bandpass 
filter" or channel-measuring filter shown in the center. The filter is shown driven all 
the time, no matter which signal "1" or "2" is transmitted. For optimum performance 
the squared-amplitude passband characteristic of this filter should match the fluctuation 
power spectrum of the path. In each of the two final bandpass filters all channel 
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contributions add in phase with the right amplitude weighting, as follows. Note that the 
input to each bandpass filter has the desired amplitude (product of current value A and 
average A) and that the phase is constant at <t>^ (that of the oscillator at v> A ) plus the 
small fluctuation (<t>-<j>), By driving all D elements that are like that in Fig. 1-4 from the 
same oscillator at w A> all channels can be added together before the bandpass filtering, 
thereby achieving the desired phase realignment. 



DIVERSITY CHANNELS 

Fig. 1-2. A modification of one element of the structure of Fig. I- lb to include 
channel amplitude weighting and phase correction. (The final narrow 
BPF's and linear envelope detectors appear only once. All else is 
iterated D times to form the complete receiver.) 

At least two questions have to be investigated further to determine whether this 
receiver is optimum for the optical channel. Note that in going from Fig. 1-1 to Fig. 1-2 
some sort of "nonlinear operation" has been retained, and the first question concerns 
whether this is a requirement in the optical transmission case. If the optical channel 
had been nonfluctuating, the optimum receiver (under the assumption of only additive 
white Gaussian noise) would have been the single elements within the dotted line of 
Fig. I— 1 , with the nonlinearity replaced by a linearity. If the channel fading statistics 
were Gaussian, the nonlinearity in Fig. 1-2 becomes a linearity and so all diversity 
channels may be combined before bandpass filtering (as just described) with the nonlinear 
gain element taken out. One may conjecture that when the fading statistics are log- 
normal, as is presumed to be the case for the optical atmospheric channel, some form 
of nonlinear gain operation may be needed. This is tentatively shown by the dashed-line 
"nonlinear gain" elements in Fig. 1-2. A second point to be mentioned is that, as shown, 
an equal smoothing time for amplitude and phase is implied. Equality of fluctuation rates 
of amplitude and phase, a property of Gaussian fading, is not a property of the optical 
channel. These points and presumably others, concerning the required modification of 
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the Fig. 1-2 receiver, optimum for Gaussian fading, to convert it to a receiver that is 
optimal for log-normal fading ought to be investigated further. 

The biggest problem in actually building such a system would appear to be the rela- 
tive frequency stability of the transmitter and receiver. Assuming that the transmitter 
frequency is constant, we note that the various conversions ahead of the "very narrow- 
band filter" must translate the signal spectrum to an accuracy of, say, one-tenth of the 
bandwidth l/T^, Under the assumption of atmospheric fluctuation rates of roughly 1 kHz, 
this implies frequency stability to 100 Hz. This order of raw laser frequency stability 
might be achievable, but it appears that the use of a locked loop (completed via path A-B 
in Fig. 1-2) would be more convenient. In space applications this approach will be 
required because of the large Doppler shifts. 

If frequency stability problems become too great, one can always fall back on the 
transmitted-reference technique. Here, we receive on each detector element not only 
the modulated signal but also a reference sinusoid transmitted at the desired nearby fre- 
quency (close enough so that fluctuations are highly correlated). The receiver reference 
tone is then injected at point C of each diversity element (Fig. 1-2) in place of the signal 
shown. This wastes part of the transmitter power, but might be interesting in certain 
cases. In particular, this artifice appears attractive for the spatially continuous system 
to be discussed next. 

B. Spatially Continuous Version 

In Fig. 1-2 the key element was the feed-forward loop shown isolated in Fig. I-3a. 
This operation had to be iterated as many times as there were diversity elements, in the 
structure proposed in the preceding section. A more elegant approach and one that 
might become quite necessary if the number of diversity elements becomes truly large 
is to set up such a feed-forward loop of the entire spatially distributed field. To do this 
in a manner exactly analogous to Fig. I-3a requires that the received field (amplitude 
and phase as a function of time and two spatial coordinates) be both frequency- shifted 
and time-averaged before being injected as a "local oscillator" in a final conversion. 
Of these two operations, only the latter appears possible. Therefore, we must fall 
back on the transmitted-reference technique just mentioned. For this situation, 

Fig. I-3a becomes Fig. I-3b. 

The corresponding two-dimensional receiver is shown in Fig. I-3c. The pilot signal 
excites a laser amplifier whose bandwidth (at a frequency shown as co^ + to,,) is the 
required 1 /T M and whose apertures at the two ends are adjusted so that it can amplify 
many modes simultaneously (each mode corresponding to a different spatial Fourier 
component of the received wavefront). We require an amplifier, rather than simply 
a resonant cavity, so that the local-oscillator signal amplitude will be very large 


* Systems of the Rake 2 type, based on Fig. 1-2, have been built at megacycle frequencies 
embodying 50 or more diversity elements. 
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Fig. 1-3. 


(a) Feed-forward loop. 

(b) Alternative approach using transmitted 
pilot sinusoid (unmodulated at Wj + w 2 ). 

(c) Spatially continuous version of Fig. 1-2, 
using the method of deriving the conver- 
sion signal given in Fig. I-3b. 


compared with that of the arriving signal, thereby minimizing the effect of noise in the 
final detecting element. 

If these operations can be carried out, then, at all points across the aperture 
ahead of the final lens that focuses energy onto the difference -frequency detector, the 
difference -frequency component has almost constant phase, and an amplitude distribution 
in which those points receiving most strongly are accentuated, and those receiving most 
weakly are suppressed. 
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